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Abstract 

Charge Symmetry Breaking and 
Nuclear Pion Production Reactions 

Daniel R. Bohon 

Chair of the Supervisory Committee: 
Professor Gerald A. Miller 
Department of Physics 

Large momentum transfer reactions such as pion production represent the frontier of Chi- 
ral Perturbation Theory and must be understood before more complex reactions can be 
considered. Pion production is also interesting in its own right, one application being the 
hadronic extraction of a charge symmetry breaking parameter: the contribution of the down- 
up quark mass difference to the neutron-proton mass difference. This dissertation reports on 
two primary projects: (1) a calculation of the charge symmetry breaking forward-backward 
asymmetry of the differential cross section of the np — > d7r° reaction, and (2) the devel- 
opment of a new theoretical framework addressing the issue of reducibility in the impulse 
approximation's contribution to pion production. It is shown that the traditional one-body 
impulse approximation must be replaced by a two-body operator which makes a larger 
contribution to s-wave pion production. 
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GLOSSARY 



CHPT: Chiral Perturbation Theory, an EFT for low energy hadronic iteractions. 

CONTACT TERM: An interaction vertex which can be thought of as originating from the 
integrating out of an intermediate heavy field. 

CS(B): Charge Symmetry (Breaking), the approximate invariance (small symmetry break- 
ing) of nature under the interchange of up and down quarks. 

EFT: Effective Field Theory, a quantum field theory formulated on a larger length scale 
than the fundamental theory. 

ERA: Free Recoil Approximation, the use of the approximations = mj^miy and « 
in the pion production kernel when including the distorted wave initial states. 

(H)BxPT: (Heavy) Baryon Chiral Perturbation Theory, inclusion of (non-)relativistic 
baryons (neutron, proton, and sometimes Delta) to xPT. 

lA: Impulse Approximation, the pion production diagram in which one nucleon is a 
spectator and the other produces a pion directly. 

ISOSPIN: A quantum number, similar to spin, which is used to describe systems in which 
the light quarks are treated as identical particles. 

LEC: Low Energy Constant, a parameter/coupling constant for a term of an EFT, often 
a non-renormalizable contact term. 
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MCS: Momentum Counting Scheme, a reorganized power counting for pion production 
reactions accounting for the requisite large three-momentum. 

(N)LO: (Next-to-)Leading Order, a diagram's location within the expansion of an EFT. 

OPE: One-Pion-Exchange, self-explanatory nucleon-nucleon interaction. 

RS: Rescattering, the pion production diagram in which one nucleon produces a pion 
which is then scattered on-shell by the other nucleon. 

SEAGULL: A vertex with two nucleons (incoming and outgoing) and two pions. 

THRESHOLD: In the context of pion production, the kinematical situation where a pion 
can just barely be produced; that is, the reaction in which the final state pion has 
zero three-momentum. 

WT: Weinberg- Tomozawa, a ttuNN vertex which is fixed by chiral symmetry and scales 
as the sum of the incoming and outgoing pion momenta. 
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Chapter 1 
INTRODUCTION 

Pion production reactions occur when two nucleons (protons or neutrons) collide at 
sufficiently high energy to produce a pion. When the collision is that of a single neutron 
with a single proton, it is possible for them to collapse into a bound ffiial-state deuteron after 
producing the pion. A variety of different pion production reactions have been realized in 
laboratories; this thesis focuses on experiments with proton or neutron beams and liquid or 
gas-jet hydrogen targets. Physicists study pion production for many reasons. For example, 
this thesis will describe how a concert of experimental and theoretical effort can result in 
the extraction of the fundamental parameter 5mj\f, the portion of the neutron-proton mass 
difference attributable to the down-up quark mass difference. 

This thesis begins in Chapter [2] with a review of the strong interaction, which is the 
primary force determining the properties of pion production reactions (cross sections, an- 
gular dependencies, energy dependencies, etc.). To achieve a perfect description of these 
observables, one should perform scattering theory calculations making use of the poten- 
tial energy as determined from the fundamental theory of the strong interaction. Quantum 
Chromodynamics (QCD). This fundamental theory is written in terms of quark and gluon 
degrees of freedom. At the same time, it is well-known that QCD exhibits "color confine- 
ment" [TJ. Among other things, this means that objects with non-zero color charge are 
never observed as well-defined asymptotic states. This is in contrast with Quantum Elec- 
trodynamics (QED) where a proton is a clear example of a well-defined asymptotic state 
with electric charge. 

Another property of QCD that affects its applicability is its non-perturbative nature at 
low energiesj^ This property is connected with confinement in an incompletely-understood 

^ "Low energies" is meant to exclude the high energy processes which can be found in stars and some 
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manner. Again, in contrast with QED, which has a couphng constant oem ~ 1/137 in terms 
of which one is able to perform perturbation theory, QCD has a coupling constant as = 0{1) 
preventing such an expansion. Thus, one must either perform nonperturbative calculations 
such as lattice-regularized QCD (LQCD) or else describe the interactions of color neutral 
hadrons, where a separate perturbation theory based on an effective field theory approach 
is possible. 

In this thesis the latter approach is taken and applied to the calculation of pion pro- 
duction observables which are then compared with experiment. After decades of model- 
dependent perturbative theories of the strong interaction, a theory based on the chiral 
symmetry of QCD was developed, called Chiral Perturbation Theory (^PT). Chiral Per- 
turbation Theory is guaranteed to reproduce the physics of QCD in a model-independent 
way, provided one can identify a convergent power counting scheme. Chapter [2] provides a 
development of the aspects of relevant for pion production. 

The next topic, found in Chapter [3] is a review of pion production. Much effort (both 
experimental and theoretical) has been devoted to its study over the past half-century, and 
Chapter [3] provides motivation for this effort as well as an overview of the reaction's features. 
The salient feature is the large collision momentum p = ^Jrf^^m^ that is required in order 
to produce a pion. A variety of reaction channels (particle and angular momentum) are 
possible, ranging from relatively-well understood to quite poorly understood. This work 
will focus on the near threshold up — t- dix^ reaction. 

Chapters |4] to [6j each based on a published paper, apply xPT to the calculation of 
several pion production observables. Chapter |4] [2j is an update of a theoretical calculation 
[3j from 2000 which relates bm^ to the asymmetry in the differential cross section of a 
2003 experiment [3]. The calculation depends heavily on a particular amplitude which was 
updated [5j in 2006. In Chapter [4| the updated calculation of that amplitude is extended 
off-threshold and, using a model estimate for dniN, a new post-diction is made for the 
asymmetry. The theoretical calculation of the asymmetry is found to be ~ 2— 3 experimental 



colliders, and instead refer to nuclei in, say, this piece of paper. 



3 



uncertaintie^ above the experimental value. 

While performing this calculation, the authors became interested in a formalism issue 
which is similar in nature to Ref. [5], but concerns a different amplitude, the impulse 
approximation (lA). In the lA, pions are produced directly from a single nucleon without 
interacting with the spectator nucleon. A naive attempt at a solution is introduced in 
Chapter |4j but it is clear that further investigation is required. 

This investigation is performed in Chapter [5] [B]. This chapter begins with a clearer 
identification of the problem with the traditional formalism: strikingly different results 
are obtained for the lA amplitude before and after application of the Lippmann-Schwinger 
equation to the external nucleon-nucleon wave functions. The same solution as in Chapter |4] 
is pursued but more subtle modifications are made in order to correctly account for the 
large scale p. Unfortunately, the large number of theoretical complications prohibit a clear 
resolution of the issue in this nai've manner. 

A definitive conclusion is reached in Chapter |6] [7], which attacks the problem from its 
root in the reduction from a covariant 4D problem to a traditional 3D one. It is shown 
that the traditional formalism for the lA is itself an approximation that is only valid in 
the absence of initial-state interactions. For the full distorted-wave calculation, the lA is 
correctly included by using the Lippmann-Schwinger equation to turn it into a two-body 
operator. The reduction to such a simple interpretation involves several approximations 
which are shown to be controlled in the Appendices. The final conclusion is that the lA 
amplitude is enhanced by a factor of approximately two. 

Finally, Chapter [7] gives a review of the topics discussed in this thesis. Problems encoun- 
tered in the calculations of Chapters |4] to [6] are recalled and possible solutions listed. Having 
successfully shown in Chapter [6] that the I A is larger than was previously thought, at least 
two calculations need to be updated: the p-wave A^A^ — t- dir amplitudes (and therefore, the 
np — )• d-K^ asymmetry) and the pp — )■ ppir^ total cross section. 



^The theory is not yet well-defined enough to assign an uncertainty. 
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Chapter 2 

CHIRAL PERTURBATION THEORY 



Because of the non-perturbative nature of QCD, the theory cannot be used directljj^to 
calculate pion production observables. One way to overcome this limitation is to use an 
effective field theory (EFT) formulated in terms of the hadronic bound states of quarks and 
gluons, i.e. pions and nucleon^ [8j. This theory, called chiral perturbation theory (xPT), 
is guaranteed to include all of QCD's features through a perturbative expansion containing 
all terms consistent with the symmetries of QCD [13 ITTl [T^ I13j . The truth behind the 
usefulness of EFT is a separation of scales; the long distance (low energy) physics of pions 
and nucleons does not depend greatly on the short distance (high energy) physics of quarks 
and gluons. Taking advantage of this separation of scales, the expansion of xPT is organized 
in terms of the small external momenta divided by the large hadronic scale which happens 
to also be the scale associated with the spontaneous breaking of the chiral symmetry of 
QCD, A^. The usefulness of xPT is primarily the fact that (once convergence is shown), it 
provides model-independent predictions for strong interaction phenomena. 



In Sec. 2.1 , an overview of QCD and its symmetries is given. Mesonic xPT is developed 



and the "power counting" scheme that provides predictive power is explained in Sec. |2.2| 



Section 2.3 expands the theory to include nucleons and a difficulty therein is overcome. 



Finally, in Section 2.4, Charge Symmetry Breaking operators are expressed within the 
formalism of xPT. 



^Lattice QCD does provide means to accomplish this, but computational limitations make pion production 
calculations presently impractical. 

■^"Nucleons" refers to protons and/or neutrons which are put into the same multiplet in xPT. 
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2.1 Quantum Chromodynamics 

2.1.1 The Lagrangian of QCD 

In this section, a brief description is given for the experimentally observed spontaneous 
breakdown of the SU(2)i(8)SU(2)/j symmetry of QCD in the limit of massless u and d 
quarks to SU(2)v. This material has been summarized often, and this summary borrows 
heavily from Ref. |14j . 

To begin, consider a quark field qf{x) with flavor index / (color and Dirac indices are 
suppressed) that transforms under a local color gauge transformation as, 

g/(x) ^exp (^-iJ2@\x)^^ qf{x)^U{x)qf{x), (2.1) 

where A* represents the eight Gell-Mann matrices which generate SU(3) and 0*(x) are func- 
tions of spacetime representing a particular gauge transformation. The Lagrange density 
for the matter fields is then, 

>Cmattor = ^^Qfi^) (0 " qf{x), (2.2) 

/ 

where the flavor sum runs over the six different quarks (u, d, s, c, b, and t) and the covariant 
derivative D is required by gauge symmetry to transform as, 

D^qf^UD^qj. (2.3) 

This requirement is met by introducing a gauge field Af^ (via = + igs^^^i) that 
transforms as, 

UAfJA^ + -d^UU\ (2.4) 

where g^^ is the running strong coupling constant. Also note that the gauge fields are 
matrices in color-space. 



A^ — A^X . 

i 



(2.5) 
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The gauge field is tlien made dynamical by including the term, 

>Cgauge = -^E^;-^^^^ (2.6) 

i 

where the field-strength tensor = d^A]^ — d^A\, — gsP-'^AjtA^ has a third term that 
is unique to non-Abelian gauge theories. This term is written in terms of the structure 
constants of SU(3), 

fi'^ = 1tv([A\A-'] A'^) . (2.7) 

Finally, note that another term, quadratic in the field-strength tensor and totally antisym- 
metric in spacetime indices, is allowed by gauge invariance. This term is important because 
it violates CP invariance; that is, it is not invariant under simultaneous charge conjugation 
and parity transformations. If its coupling constant were non-zero, then the strong force 
would break CP symmetry; however, experimentally, no evidence for this term has been 
found [T3] . 

2.1.2 The Chiral Symmetry of QCD 

Apart from the coupling gs, the only parameters of QCD are the quark masses. The 
"current" quark masse^are listed in Table 2.1 It is clear that the masses of the up and 



down (together, "light") quarks are negligible compared with the hadronic scale of ~ 1 GeV. 
For this reason, it is beneficial to consider the massless, or chiral, limit of QCD as a starting 
point for calculation. The asymptotic states of interest in pion production are composed 
entirely (in the valence sense) of light quarks. Thus the effects of the four heavier quarks, 
with masses mq, are seen through either electroweak or virtual pair production processes. 
In moving to an effective field theory, these processes are "integrated out" and replaced 
with contact interactions that are suppressed by factors of l/mq. 

The strange quark mass is such that it could feasibly be included as a third light quark, 
and indeed this is sometimes done. However, the review article |T3] discusses the fact that 



■^Because quarks are never observed as asymptotic states, any statement of their masses must come with 
a description of what that mass means. Referred to here are the "current" quark masses (the parameters 
in the Lagrangian) in contrast with the model-inspired "constituent" quark masses. 
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Table 2.1: Current quark masses in the MS renormalization scheme taken from the PDG 
|16j . The values shown are the masses mq{fj,) as evaluated at a particular scale /x. For the 
n, d, and s quarks, the value shown corresponds to fi = 2 GeV. For the c, b, and t quarks 
the pole mass [mq{fi = niq)] is used. 



Quark Mass (MeV) 

up 1.7 to 3.3 

down 4.1 to 5.8 

strange 80 to 130 

charm 1,180 to 1,340 

bottom 4,130 to 4,370 

top 172, 000 ±1,600 



the resulting EFT is relatively slowly convergent. For example, Ref. |17] showed that 
the one-loop expansion of the octet and decouplet baryon masses in SU(3) HB^PT does 
not converge. This conclusion was reached by using LQCD to investigate the quark mass 
dependence of the baryon masses. On the other hand, a relatively new, covariant version of 
SU(3) BxPT was used in Ref. [18] to calculate baryon magnetic moments with reasonable 
convergence. The application of EFT to pion production is already expected to be slowly 
convergent due to the large external momentum p. Therefore, to avoid further convergence 
issues, strangeness is not considered in the pion production calculations of this thesis. 

When the masses of the light quarks are neglected, there is nothing in the Lagrangian 
to distinguish one from the other. For this reason they are often grouped into an "isospin" 

(A 

doublet (7=1 I on which one may perform isospin rotations with the Pauli matrices. The 
terminology used is that the light quarks are "isospin- 1/2" particles. As will be discussed 



in Section 2.3 the proton and neutron are also isospin-1/2 particles due to their similar 
masses. The symbol used to represent the quantum number of isospin is T (many authors 
use / instead). 
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As shown in this and the next paragraph, the Lagrangian for chiral QCD manifests not 
only Lorentz and SU(3)c gauge symmetries, but also a global "chiral" symmetry. Consider 
the projection operators, 

Pl = \{1-1^), Pij=^(l + 75) (2.8) 

which satisfy + Pr = 1, Pi = P^, = Pr, and PlPr = PrPl = 0. These projectors 
separate the helicity (p • S) eigenstates of the massless, free Dirac Equation: qi = PlQ, 

I UL 

qpt = Prq. Utilizing these definitions and adopting vector flavor notation, 

^Scr = {Ql i0 qL + Qr 1$ Qr) ■ (2.9) 

Note that the left- and right-handed sectors have decoupled and that the neglected mass 
term provides the coupling, 

Cm^s = -{qiMqR + qRMqL) (2.10) 

with M = d[ag{mu,md)- 

Chiral symmetry can now be formalized with the definition of the global chiral transfor- 
mation, 

3 



qi exp 



QL, (2.11) 



a=l 

and similar for the right-handed fields. Note the appearance of the Pauli matrices that 
generate the SU(2) symmetry group. The U(l) transformation is parametrized by d^. The 



Lagrangian of Eq. (2.9) is invariant under each of SU(2)l, SU(2)r, U(1)l, and U(1)r 
separately. Due to the symmetry breaking pattern observed in nature (discussed below), it 
is beneficial to combine the left and right transformations to form vector and axial vector 
transformations. To accomplish this, consider the Noether [TDJ symmetry currents, 

dSC _ Ta 

= dim = ''''''' 
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and similarly for the right-handed currents. In these equations, bL represents the variation 
in the Lagrange density due to a locals infinitesimal chiral transformation. Now form linear 
combinations, 



(2.13) 



Except for A^^, each of these currents are conserved in the chiral limit of the strong force. 
Quantum fluctuations destroy the classical conservation of A^, and instead one finds [20j . 



d„.A'' 



3ffi 
327r2 



'-fii/pa 



(2.14) 



Another way to understand this divergence is to notice that although the Lagrangian itself 
is invariant under the iso-singlet axial transformation, the measure of the path integral is 
not. 

Because later chapters will discuss the symmetry- breaking role of the quark masses, they 
should be considered briefly now. Again, using Noether's Theorem, one finds that the terms 



of Eq. (2.10) cause three of the currents to develop a divergence, 



d^V^' = iq [M, 
d^A^,=iq{M,^]j,q 
^f,A^' = 2iqMj5q. 



(2.15) 



Thus the extent to which the quark masses are non-zero causes explicit breaking of the 
axial symmetries, and the extent to which — niu ^ causes explicit breaking of the 
vector, iso-vector symmetry. The fact that the divergence of A^ is proportional to the 
pseudoscalar density is referred to as the Partially Conserved Axial Current relation, and has 
been understood (at least macroscopically) since Gell-Mann first described the "Eightfold 
Way" in 1964 [21]. 
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2.1.3 Spontaneous Symmetry Breaking of QCD 

Upon first thought, one should expect to see a spectrum of hadrons corresponding to the 
irreducible representations of the symmetry group G, |22j 

G = SU(2)l®SU(2)/j®U(1)v. (2.16) 

The multiplets are classified as mesons or baryons based on how they transform under U(l)y 
and there should be small splittings within each multiplet resulting from the non-zero quark 
masses. If the spectrum is invariant under separate left- and right-handed transformations, 
the multiplets should contain approximately degenerate positive and negative parity par- 
ticles. The experimentally observed spectrum does not appear to satisfy this rule. For 
example, the negative-parity proton has a mass of 1535 MeV, almost 600 MeV above the 
proton. 

The fact that the full symmetry group is not manifest in the hadronic spectrum implies 
spontaneous symmetry breaking. It can be shown that a scalar quark condensate provides 
such a mechanism, 

(^,,L(a^)gi,i?(x)) = A^^i,-, (2.17) 

where (•) denotes the expectation value of a quantity in the interacting QCD vacuum, and 
A is a quantity with dimension of mass. Under a chiral transformation, one finds 

{qj,L{x)qi,R{^)) ^ K^RiAkLl^ = A^^^j. (2.18) 

Now, it is clear that Sjj = 6ij if and only if L = R, which implies that the scalar quark 
condensate is only invariant under the V = L + R (called the "isospin" ) symmetry trans- 
formation. Another way to state this is that a non-zero value for the condensate implies 
that the symmetry group of the ground stat^is smaller than the symmetry group of the 
Lagrangian, 

{q{x)q{x)) / =^ SU(2)i » S\J{2)r ^ SU(2)y. (2.19) 



*The scalar quark condensate has the same quantum numbers as the vacuum 
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Thus, the ground state of the strong force is only invariant under S\J{2)y^\J{l)y. Gold- 
stone's theorem |23l [2^ then states that the broken generators manifest as pseudoscalar 
bosons with small (non-zero due to the explicit breaking of Eq. ( 2.15| )) mass. These cri- 
terion are of course met by the pions. Note that physically, a scalar quark condensate 
means that it is more energetically favorable for the vacuum to be constantly fluctuating 
into qq pairs. As this document is not meant to be be a complete review of QCD, the 
interested reader is encouraged to see more group-theoretical proofs of these statements in 
the already- mentioned Ref. [ TS] . 



2.2 xPT for Mesons 

2.2.1 Leading Order Kinetic Term 

We have seen that the pions arise as the Goldstone bosons of the spontaneous symmetry 
breaking of SU(2)i (g) SU(2)^. To be more specific, the pion fields are fluctuations of the qq 
operator about the vacuum. Therefore, when including them in the EFT Lagrangian, one 
must ensure that the pions transform in the same way as qq under the full symmetry group 



G of Eq. (2.16). The vacuum is defined by the (non-zero) chiral condensate which remains 
invariant under the subgroup SU(2)v. The physical pion states fill out the baryon-number- 
zero, isospin-one representation of the SU(2)y ® U(l)i/ symmetry group. 

The pion fields can be written as a 2 x 2 special unitary matrix U with the correct 
transformation properties, 

U{x) RU{x)L^. (2.20) 
Although there is freedom in how U is parametrized, the common convention is chosen here, 

U{x) = exp , (2.21) 



where 7r(x) = ^a''^a{x)Ta contains the pion fields, and /tt = 93 MeV is the pion decay 
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constant. Note that in terms of the charge states, the pion matrix is given b; 

TT = I I , (2.23) 




so that in Cartesian coordinates vr^ = (vri it 7r2)/\/2- The value of is experimentaUy 
determined from the leptonic decay tt — )• //z^, 

{fi~{kp^{q)\H\7T-{p)) = (01^X7^(1 -75)d|vr-) x M'^^ 

(2.24) 

= ifnPti X -^EW 

where A^^^v represents the perturbatively calculable electroweak contribution to the pro- 
cess. 

In order to construct Chiral Perturbation Theory, one follows the principle that any 



term consistent with the symmetry group G of Eq. (2.16) must appear in the Lagrangian. 
All such terms will contain derivatives 9^C/ which will scale like the typical (also referred to 
as external) momentum, Pext- Recall that an effective theory is only an approximation to a 
more complete high energy theory. Here this means that is only valid for low energy 

processes with pext ^ where ~ 1 GeV is the energy scale at which chiral symmetry 
is restored. Therefore, one is able to order the terms in the Lagrangian according to powers 
of X = J'ext / A;,^- This organization scheme is called power counting, and is central to xPT. 

Although terms with mass dimension greater than four most certainly appear in the 
Lagrangian, xPT remains renormalizable order-by-order in the expansion. Li other words, 
at next-to-leading order (NLO), loops built from leading-order (LO) interactions are renor- 
malized by NLO tree-level diagrams. This procedure was worked out in great detail by 
Gasser &: Leutweyler 0. A well-thought-out formalism for power counting particular dia- 
grams exists [8j, but will not be presented at this time because a modified power counting 
scheme is necessary for pion production reactions. 



"■The charge states are formed to satisfy 

(2.22) 



+2/3 \ / 7r° %/27r-' 
-l/2j'lV27r- -tt" 



X tt" +1 X \/27r-' 
-1 X V^TV' X -7r° 



Practically, one expands in powers of Pcxt/47r/^, since factors of 47r/,r ~ A^^ commonly appear in loop 
calculations 
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The leading order Lagrangian consists of a single term with two derivatives to which all 
other possible forms can be reduced (up to unphysical total derivatives), 

£ = :^TV (d^.Ud^'uA + (2.25) 



4 

where the prefactor was included to ensure consistency of the kinetic terms with the La- 



grangian for free scalars: ^d^(f>d^<p. To see this, expand out Eq. (2.25) and look at the 
quadratic terms. 



C ^ fU, ( X (-z)^^^ + ...)+ ... 

= ^d^TTad^^TTbTriTan) + ... (2-26) 
= ^d^TTad'^TTa + ... 

Nevertheless, this is not the complete leading order Lagrangian because the quark masses, 
which in reality are non-zero, have not yet been included. 

2.2.2 Inclusion of Quark Masses 

One is not able to use the same symmetry principle to include the quark mass term because 
the quark masses explicitly break chiral symmetry: q^Mqi — )• qj^R^ MLq^. However, if one 
pretended that the mass matrix transformed like M — )■ RML\ then the mass term would 
be invariant. Thus one constructs a corresponding chrially invariant term in xPT, 

C D ^Tr(MC/t + UM^), (2.27) 

where D is understood to mean "includes", the was included for convenience (see 

below) and i? is a mass scale used to make the term have correct mass dimension. Note that 
the antisymmetric combination of M and U is excluded by paritjj^ Again, by expanding 



Eq. (2.27), one is able to connect B with the mass of the pion (here the isospin limit 



m„ = ma = ruq is employed; isospin breaking will be considered in later chapters), 

^Tr(MC/1' + UM'^) = const - Bm^vr^ + ... (2.28) 



'^Under parity 7r(t, x) — >■ — 7r(t, — x), so ?7(t,x) — )■ L'^^(t, — x). 
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Identifying the quadratic term with the mass, it is clear that ?ti-^/2 = Bniq. As might 
be expected, B can also be shown to be proportional to the aforementioned scalar quark 
condensate [Hj. 

As a brief aside, although flavor SU(3) has not been considered in this chapter as a good 
symmetry, one is free to do this by replacing the 2x2 pion matrix with the 3x3 meson 
matrix, 



a=l 



V3 



V2' 



IT ' 



V2K+\ 



(2.29) 



.21 , 
V3 I 



Performing an expansion and looking at the quadratic terms, one finds the following ex- 
pressions for the meson masses, 

,2 _ , _ ^ ^2_^ ^ B{ms + mu) 



mz = B{mu + md) 



B 



m„ 



(4ms +mu + md) 



B{ms + md). 



(2.30) 



3 

These equations can be combined to obtain a rough estimate for the light quark mass 
m = [mu + md) /2 in terms of meson masses and the strange quark masa 



m 



2m|. 



771: 



2 ■ 



(2.31) 



Finally, it is interesting to look at the mass splitting between the charged and neutral 
mesons. To first order, there is no quark mass contribution to this splitting for the pions: 



m 



m 



<^EM + 0. However, this electromagnetic splitting should be the same for the 



kaons, which do have a quark mass contribution. 



- mj^o = Sem + B{mu - md). 



(2.32) 



Thus, using Eq. (2.30) one can derive the following relations, 

,2 



mu 
md 

ms 
rrid 



m 



K+ 



m 



Im 



2ml,) 



m 



m 



RO 
m 



R+ 



m\+ + m1_ 

2 

m . 



0.56 



(2.33) 



m 



RO 



20.1. 



The strange quark mass can, in turn, be determined by a variety of methods; perhaps the most common 
modern method is lattice hadron spectroscopy with non-perturbative renormahzation. For example, see 
Ref. [25]. 
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Note that these relations are not exact for two reasons: they were derived at leading order 
in the EFT, and the EFT itself is based on an inexact symmetry. 

Let us now return to the construction of the xPT Lagrangian. It is important to note 
that the square of the pion mass is proportional to the quark mass. For low energy processes, 
one has m.j^ ~ pext and thus m^r/A^ is considered perturbatively small. This means that a 
single insertion of the quark mass matrix is counted as being "of the same order" as Pext- 



Therefore the mass term of Eq. (2.27) is included in the leading order Lagrangian, 



(2.34) 



where the superscript (2) refers to the fact that each of the terms in this expression scale 
like 

Pgxf Once again, note that this expression is not invariant under chiral symmetry (due 
to the mass terms), but that the breaking of chiral symmetry is motivated by the way in 
which QCD breaks chiral symmetry. 



2.2.3 Beyond Leading Order: Low Energy Constants 

The building blocks of the mesonic xPT Lagrangian are [/, d^U , and M, and all terms 
consistent with the symmetry group G must appear at some order. It is interesting to note 
that there are no terms of order p^^^ because all the Lorentz indices must contract and 
M ~ Pexf further illustrate the EFT principle, the NLO mesonic Lagrangian [TU] is 
displayed, 

2 



£(^) = Li |Tr D^U{D''U)^ } + LsTr D^U{D^U)^ Tr D^'U{D''U) 



+ L3Tr 



D^U{D^'U)^DyU{D''U) 



+ 2BLiTT 



D^UiDf'U)^ 



Tr 



+ 2BL5TT 



D^.UiD^'U)^ {aIU"^ + C/M"^) + AB'^Lq |Tr Mf/"^ + C/M^ | 



+ 452^7 |Tr MU^ - UM^ I + 4.B'^LsTv UM^UM^ + MU^MU^ 



(2.35) 



where terms without U have been ignored because they are unphysical. In Eq. (2.35), the 



Li are referred to as Low Energy Constants (LECs). It is in these LECs that xPT rises 
or falls as a predictive theory. The idea is that one can calculate a quantity (say, the tttt 
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scattering length) to one-loop order using the vertices of £^^^ and absorbing the divergences 
into the operable LECs from C'^^\ The renormalized values of the LECs are then either 
obtained from experiment or calculated with LQCD. 

External fields (such as the electromagnetic field) can be included in the QCD Lagrangian 
by adding the terms ^ = ^iV^f/j'j^YQ — '^^^t'^Tst^' When working in 

terms of right- and left-handed quark fields, the currents are rewritten as t;^ = + l^, 
cipL = ffi — lf_i', where — t- Rr^B) — id^RR^ , — t- Ll^V^ — id^LL'^ under a local chiral 



symmetry transformation. The derivatives in Eq. (2.35) have been written with a capital D 
to indicate that these external fields can be included in xPT by using a covariant derivative 
dfJJ — D^U = dfj_U — ir^U + iUl^ which transforms exactly as d^U does. 

In the sense that is a low energy theory which does not consider heavy quarks (or 
even baryons!) explicitly, the LECs are sometime referred to as "contact terms." As the 
language implies, effects of heavy particles occur on short time scales and when one probes a 
system with a low energy, such time scales cannot be resolved but instead appear shrunk to 
a (contact) point. By analogy, in Fermi's electroweak EFT, an interaction of two electrons 
via a W boson is expressed as a single 4-pt vertex with Gp sts the LEG. 



Of course, the procedure described following Eq. (2.35) does nothing to validate the 
theory; further predictions must be made. Indeed, distinct combinations of LECs contribute 
to several experimentally independent observables including the pion mass, polarizability, 
and decay constant; the 77 — )• tt^tt^ reaction; and low energy pion-pion scattering. As long 
as the renormalized LECs not so large as to prohibit the use of perturbation theory, the 
theory is generally accepted as useful. This is clearly found to be the case for two-flavor 
mesonic ^PT. A modern survey of the statements in this paragraph can be found in Ref. 

ESI. 



^Iso-scalar vector, scalar, and pseudoscalar terms can be included similarly. 
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2.3 xPT for Baryons 

2.3.1 Naive Leading Order Lagrangian 

The next task of this chapter is to introduce baryons to (creating BxPT). Because 

their masses are nearly identical, the nucleons are placed into a isospin doublet, 

(2.36) 

It is known phenomenologically that pions mediate the long-range part of the nucleon- 
nucleon interaction [27]. It is useful to include this behavior when adding nucleons to 
xPT. This can be accomplished by gauging the nucleon kinetic term, promoting the chiral 
transformation from global to local. The local nature of the transformation comes from 
the pion fields by including U{x) in the manner described below. Because one is free to 
perform field redefinitions, it is useful to choose a basis such that the nucleon doublet 
transforms like tp — t- K(L, R, U {x))^l^, where K is some undefined function of the L and R 
chiral transformations, in addition to U{x). The correct pion-nucleon interaction will then 
make the covariant derivative, D^ip, transform like the nucleon itself: D^ip — )• KD^^l^. The 
expressions which satisfy these constraints are, 



DA 



(2.37) 



u{x) Ru{x)K-^{x) (2.38) 
u{x) = ^/U(x) = e^"(^)/2/- . (2.39) 



From Eqs. ( |2.38| ) and ( |2.39| ) one finds that K = VrULJ ^R\/U. To reproduce this 
parametrization, use the field redefinition xp'^ = utpi, ip'j^ = u^iJjr the free form (0^. 
Though this definition may seem esoteric, it permits a simple construction of chirally in- 
variant interactions. Additionally, note that this parametrization preserves the necessary 
transformation under the 811(2);/ subgroup, ■0 — ^ ^V'l since for L = ii = F one finds that 
K = V . One final observation regarding the kinetic term is that by parity the pion-nucleon 
interaction only contains terms with even numbers of pion fields. 
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Figure 2.1: One of the loop contributions to the nucleon mass. 



Next, one constructs another chirally invariant term built from the u{x). Define, 

Uf, = i {y}d^u - udfj,u^^ , (2.40) 

which transforms like — )• Ku^K^ . This term is odd under parity so one must include a 
75 in addition to contracting with 7^ to make it a scalar. Putting all the pieces together, 
the leading order Lagrangian is written, 

£(1) =^(^iiP-mN + y7''75^^m) V', (2.41) 

where uin- ~ 1 GeV is the nucleon mass and gA ~ 1-27 is the axial-vector coupling. 

The key difference between this construction and that of mesonic xPT is the introduction 
of a new mass scale 771^ which does not vanish in the chiral limit. Schematically, this 
statement means that ~ mNtp. This introduces a subtlety to the power counting 
scheme in which loops with nucleons contain terms that are effectively promoted to lower 
order. As an example, it can be shown that in the MS scheme, the (subtracted) diagram 



of Fig. 2.1 (which contributes to the renormalized nucleon mass) contains both a term 
proportional to g\m^/ f"^ and one proportional to g\m'^mN / f"^ [14J. This is a problem if 
one is to assign the diagram a particular order in the power counting scheme. 

2.3.2 Heavy Baryon Formalism 

One way to solve the power counting difficulty brought about by the nucleon mass is to 
remove it by performing the field redefinition if) — )• e^'^'^^ip and then integrating out the 
heavy components of the nucleon field. The goal of this procedure, introduced by Ref. |28j . 
is to produce nucleons that obey d^N ~ Pext-^ at the cost of new effective interactions 
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that are suppressed by powers of l/rriAr. As shown below, the remaining dynamical field 
propagates (in the absence of interactions) with four-momentum = (£^(p) — m^v, p)- 
Since BxPT is only to be used for low momentum transfer processes, it should then be a 
good approximation to assume that <^ mj^i for all The end result is a new EFT, 
Heavy Baryon Chiral Perturbation Theory (HBxPT), that is an expansion in both p^-^t/ 
and Pe^t/mN. 

To develop HB^PT, first define projection operators, 

n± = (2.42) 

which satisfy P^j^ + P^- = 1, P^j^ = Pv±, and Py±Pyzp = with the requirements that 
v'^ = 1 and vq > 1 satisfied by the choice v = (1,0) that is used in this work. This choice 
makes it clear that one is separating out the heavy and light components the free nucleon 
spinor, 

(' VeT^ ( ^ ] e---. (2.43) 

Now one makes the field redefinition/split, 

N{x) = e*™^'^-^P^+V'(x), H{x) = e*™^'^-^P^_V'(x). (2.44) 

Note that fN{x) = N{x) and fH{x) = —H{x). The next step is to use the equations of 
motion to eliminate the heavy field H{x). The theory will then consist of solely the light 
field N{x) which propagates with exp(— i(|3 — rriMv) ■ x). This residual momentum will be 
referred to as fc^, 

Pfi = ruNVfj. + k^ (2.45) 
v-k = -^, (2.46) 

where the latter is found by squaring the former. 

The Euler-Lagrange equations give the equation of motion for the leading order La- 



grangian in Eq. (2.41), 



{lU) -mN + y7^75^Xm) V'(^) = 0. (2.47) 
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Inverting Eq. (2.44 ) , inserting the resulting expression for ip into Eq. (2.47 ) , and multiplying 



by — exp{imj^v • x), one finds 



gA 



rriNi) + iU)-mN + "^Yl^u^ ) {N{x) + H{x)) 



0. 



(2.48) 



Next, left multiply Eq. (2.48) first by and second by to create two differential 



equations. After several steps of algebra, one obtains 

9A 



ivD + ^|i^75 ) N{x) + Ulp^ + ^v u-f5 ) H{x 



(2.49) 



= 

0± - • ^^5) + (-^^ ■D-2mN + Y^±15) H{x) = 0, 
where V_\_^ = — v ■ Vv^ and various relations involving v have been used; for example 



V ■ DH{x) = - 
iv D + 



V ■ DfH{x). Now the heavy field H{x) is eliminated to find 

9A 



75 



QA 



iv ■ D + 2mN 



9A 



75 



;0j 



gA 



-V ■ ti75 



N{x) = 0. 

(2.50) 

Though the proof is technical (see Ref. [2] ) , it can be shown that these nucleon spinors are 
normalized differently than the relativistic ones, such that one needs to multiply correlation 
functions by the normalization factor \J E /rrii^jVQ for each external nucleon. 



All that remains is to organize Eq. (2.50) into an expansion in (2miy) ^ and to manipu- 
late the expression into a more commonly used form. It is clear now that the derivatives in 



the denominator of Eq. (2.50) will produce small quantities and thus one is free to perform 
the expansion. It is important to note here that this equation of motion can be generated 



from a Lagrangian identical to the LHS of Eq. ( 2.50 ) left multiplied by A^. Nevertheless, this 
Lagrangian does not contain all the necessary terms for the EFT because it was obtained 
from simply the leading order relativistic equation of motion. Other terms consistent with 
the symmetries must be included as well. With that in mind, one proceeds by writing 

9A , , V- On 



£d N 



iv ■ D + 



_75 + E 



N. 



2 . ^ (2mjv)' 

n=l 

The first (2mjv)~^ correction (called a recoi/ correction) is given by. 



(2.51) 



(2.52) 
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Note that there are two derivatives in each of these terms and regardless of whether they act 
on nucleons or pions, they will bring down factors of Pcxt making these terms ©(pp^t)- These 
expressions can be dramatically simplified using the fact that N = P^^N. Because these 
NLO interactions will be heavily used throughout this thesis, a bit more time is devoted 
now to working on them. First, consider the lp'j_ term. 



D^-{v Df - ia^"'{D^ - v ■ Dv^){D, - v D,) 



(2.53) 



-{v D) 



2 ' 

where a^'^ = ^ [7'^,7^]. The commutator of the covariant derivatives can be expressed in 
terms of the objects by use of d^{uu^) = and deletion of terms symmetric in /i -H- 
(due to its contraction with cr^'^), 

1 



(2.54) 



Next, use v = (1,0) to obtain, 

9A 



1 ■ U75 



+ 



1 



2mN 



la 



—u^Uu + i Y It • D75, uq] - — ) + 



(2.55) 



N. 



Finally, as the lower components of the nucleon spinors have been projected out (and em- 



ploying the Dirac basit ) , the leading order HBxPT Lagrangian is given by. 



£(1) = Art 



iDn a ■ u 

.2 



N. 



(2.56) 



The NLO expression includes the (2mj\[) ^ terms of Eq. (2.55) along with LEG terms 
originating from the projection of the NLO relativistic Lagrangian, 



1 



tin 



+25ciTr (M+) + ca-ug + cgu • n + Ci[a\ a^jmuj + 2Bc^ M+ - -Tr{M^ 



(2.57) 



where = Mu^ + uM'^u. The C5 term will later be shown to provide a charge symmetry 
breaking contribution to pion production. 



■"The Dirac basis is defined by: 7" 



1 
-1 





-a' 



75 



1 

1 
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2.3.3 The A Resonance 

As outlined in Ref. [13], there are two reasons why it is important to include the A(1232) 
resonance as an explicit degree of freedom in HB^PT. Firstly, the resonance lies a mere 
^ 300 MeV above the nucleon and, even for low momentum processes, should be expected 
to play an appreciable role in pion-nucleon dynamics. Secondly, it is known that the vriVA 
coupling is quite strong ^AAf-Tr ~ 2gAf7p^ In this section a brief outline is given of the rather 
technical formalism for including the S = 3/2, T = 3/2 fields in HBxPT. 

Spin-3/2 fields can be described in quantum field theory with the Rarita-Schwinger 
formalism |30j that considers the tensor product of spin-1 and spin- 1/2 fields, ipn- The 
equation of motion is given by, 

{i0 - tua) i^f,{x) = 0, (2.58) 

where the spin-1/2 component is removed with the constraint j^ip^ = 0. Additionally, 
because the A is an isospin-3/2 particle, one considers the isospin doublets ipl^ix) for i = 
1, 2, 3 with the constraint T*^/)^(a;) = projecting out the two unwanted degrees of freedom. 
In this thesis the representation advocated by [31] is used. 




(2.59) 



The general idea is the same as it was for the nucleons, one splits the ipjj^ into heavy and 
light fields and integrates out the light ones. The additional difficulty is keeping track of 



^^In p-wave ttN scattering, the inclusion of the A as an intermediate state, in addition to the naiVe nucleon 
pole, dramatically improves the calculation of the T — J — 3/2 scattering volume |29j . 
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the extra spin and isospin degrees of freedom. To this end, define the projection operators, 



p3/2 

^3/2 



1 

' 3 
1^ 

3' 



1 



(2.60) 



Because one wants to describe A^A transitions, it is necessary to project out the only the 
nucleon mass (not the fuh Delta mass), 



(2.61) 



Another important consideration is the "point invariance" required of spin-3/2 fields. With- 
out going into full detail, this is the statement that the Lagrangian must be invariant under 
a transformation -0^ — t- + a7^7;y'0'^ introducing an arbitrary admixture of spurious spin- 
1/2 fields. This invariance is accomplished by the addition of terms into the Lagrangian 
which depend on the unphysical gauge parameter A. This dependence can be absorbed into 
a field redefinition with the resulting generalized kinetic term, 

1 



A 



(2.62) 



Introducing the pions, one may write down a general expression with a covariant deriva- 
tive and axial coupling, both of which are rank-2 tensors in spin and isospin |31| . However, 
this would unnecessarily complicate matters for the purposes of this thesis. All reactions 
considered here will have nucleons as asymptotic states and will not proceed beyond NLO. 
Accordingly, diagrams containing ttAA, ttttAA, vertices will not appear up to the order 
considered. The purpose in describing the Delta lies simply in its propagator and the ttNA 
vertex. Therefore, it is sufficient to display these relevant terms only. 



After performing the steps described in Sec. 2.3.2 to eliminate the spin-1/2, isospin-1/2, 
and heavy spin-3/2 isospin-3/2 components of the full "0^, the resulting Lagrangian is once 
again expressed as in expansion in (2mAr)~^. For the choice v = (1, 0) one finds. 



C 



C 



(1) 

AfA 

(2) 
rAfA 



-T} 

hA 

2 
-1 

2mN 



2mN 



+ .. 



T^^hAid^w^N + h.c 



+ .. 



(2.63) 

(2.64) 
(2.65) 
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where 6 = — itin = 293 MeV, = ^Tr(rjM^), is the analog of for the NA 
transition, and h.c. represents the Hermitian conjugate of the preceding term. In the cal- 
culations that follow, the value = 2.1gA is adopted based on p-wave scattering data 



summarized in Ref. [29]. Finally, note that in Eq. (2.65), terms with even numbers of pions 



and couterterms were left out, again due to their absence in the present calculations. 
2.4 Charge Symmetry Breaking 

In this section. Charge Symmetry (CS) and its breaking (CSB) are defined. Examples are 
given of how CS(B) manifests in nature, as is a discussion of its inclusion in xPT. A modern 
review of this topic can be found in Ref. 



2.4-1 Definition of Charge Symmetry 

Charge symmetry is defined as the approximate invariance of a hadronic system or reaction 



under an isospin rotation of vr about the y-axis Formally, this means that the equation 



[^5,^^031=0, (2.66) 

is satisfied, where Hs is the strong Hamiltonian, and Pes = e*'^^^/^ is the charge symmetry 
operator. A simple manifestation of CSB is the degeneracy of the nucleon doublet: Pcs\p) = 
— |n), and Pcs\n) = \p) with rUn ~ mp. 

The meaning of the name "charge" symmetry becomes more clear when considering nu- 
clei with large nucleon number, A. The charge of a nucleus is proportional to the number of 
protons, Z; if n of those protons were turned into neutrons, the charge would become pro- 



portional to Z — n. However, if Eq. (2.66) is satisfied, the nuclear properties (i.e. excitation 
energies, see Fig. 2.2) will remain the same. 

A distinction which needs to be made is that charge symmetry is not the same as 
isospin invariance. This confusion is understandable since the charge and isospin operators 
are related, Q = e(l/2 + T3) (for nucleons). However, charge symmetry refers to the specific 

^■^The j/-axis is used because the z-axis of isospin is related to the electric charge, and rotating by 180° 
about the y-axis modifies the charge. 
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Figure 2.2: Example of CS: spectra of two members of the A = 11 multiplet. The notation 
is SpinP^"*y. Taken from Ref. [33]. 



situation described by Eq. (2.66) while isospin invariance is more restrictive, 

[Hs,t/2]=0, (2.67) 

where r is the total isospin vector of the system. In the context of the nucleon-nucleon in- 
teraction (where isospin invariance is referred to as "charge independence" ) , this distinction 
has been formalized with the definition of different "classes" of forces [34]. A class II force, 
for example, maintains charge symmetry but breaks charge independence, 

Vn = c[T3{i)Ts{j) - ^r(i) • t(j)], (2.68) 

where i and j are particle labels. The Coulomb interaction leads to a class II force due to 
the fact that the electric charge operator is Q oc (1 + T3). 

2.4.2 Manifestations of CSB 

Charge Symmetry Breaking has two sources, the light quark mass difference and electromag- 
netism |35]. One example where the light quark mass difference dominates is the splitting 
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of the neutron and proton masses, 

m„ - m„ 1.29 MeV 



{mn + mp)/2 938.92 MeV 



0.1%. (2.69) 



If the source of this CSB were electromagnetism, one should expect the proton to be more 
massive. Therefore, one attributes the 0.1%, plus the electromagnetic breaking, to the fact 
that d quarks are more massive than u quarks. Another example of CSB is the T = 1 nn 
(~ —19 MeV) and pp (~ —17 MeV) scattering lengths (after models are used to remove 
electromagnetic effects). Reference [36] argues that the primary source of this difference is 
the vector meson p — uj mixing due to the non-zero value of — niu- 

Isospin violation is typically dominated by electromagnetism [37]. For example, consider 
the T = 1 nn (~ —19 fm) and np (~ —24 fm) scattering lengths (again, after models are 
used to remove any residual electromagnetic effects), which at first glance would not seem to 
be affected by electromagnetism. As shown in Ref. [38j, the primary source of this difference 



is the the mass difference between the charged and neutral mesom^^ which mediate the 
nucleon-nucleon (NN) interaction. 

Although there are a multitude of examples of CSB in hadronic physics (Chapter |4] will 
discuss the angular distribution of the differential cross section of np — )• dir^ in great detail) , 
the task of incorporating the light quark mass difference and electromagnetic CSB in xPT 
will be the focus of the remainder of this chapter. 



2.4.3 CSB in xPT 

The inclusion of CSB in xPT was worked out by Ref. |39j . Let us begin by discussing CSB 
resulting from the light quark mass difference. Reference [39] begins by noting that in the 
absence of strangeness, chiral symmetry is described by the group SU(2) (g) SU(2) ~ S0(4). 
The quark mass difference term in the QCD Lagrangian is then identified as the third 

^^This mass difference, — m^o = [139.6 — 135.0] MeV = 4.6 MeV, is clearly due primarily to electro- 
magnetism since each particle has the same number of u and d valence quarks. 
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component of the S0(4) vector P = {—qrq, qi^^q), 
£QCD _ 




(2.70) 



A term with the same transformation properties must appear in xPT, and the corresponding 
LEG must be proportional to m„ — m^, 

^xPT 3 ^ J^^VraiV - ^iVvrar • vrivj , (2.71) 
where 5m at oc (m„ — rud)- To derive this form, one need not discuss the group theory of 



SO (4) since the LEG terms have already been written down in Eq. (2.57). In that equation 



one encounters the object M_|_ which contains a term proportional to m„ — m^, 

M+ = u^Mu^ + h.c. 

m,, — rrij ( i \ \ ( i 1 , , 

= • • • + 1 - TT^-^ •7r-^779+---r"3 l - wr-^ • ^ - 7779 + • • • + ^-C- 



2 V 2/^ 8/2 y 2/^ 

2/^ 



• • • H 7, 2 ( rs - — 2 T • TTTTs + . . 



(2.72) 

where the first . . . represents a term proportional to [rau^md) /2. The GSB term is traceless, 



so only the C5 term of Eq. (2.57) will contribute. Truncating the expansion at two pions 



(diagrams with four or more pions will be parametrically suppressed) one finds, 

£(2) 3 2Bc5(m„ - md) [t^ - • tt^s) , (2.73) 

which tells us that bra^ = 4:Bc^{mu—md) ■ Since B(mu+md) = , a final useful expression 
is 

5mN = 4:C5m^ ■ . (2.74) 

rriu + rud 

Sometimes the ratio that appears in this expression is referred to as a parameter e = 
{mu — md)/{mu + rud) which can be estimated from meson masses in the SU(3) limit: 
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e 



1/3. Note that the behavior is required as this term is a member of the £^^^ 



Lagrangian. 

There is also an electromagnetic CSB interaction which comes in at this order parametrized 
by the LEG Sm^. Reference [39] describes how hard photon exchanges between quarks 
break the S0(4) symmetry in a manner described by the 34 component (qi'j^Tsq) of a 
rank-2 anti-symmetric tensor, 



Since the electromagnetic terms in xPT were not derived in the above sections, connection 
with the corresponding LEG, /2, is not made at this time. 




(2.75) 



so that the corresponding term is. 




(2.76) 
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Chapter 3 

PION PRODUCTION IN NUCLEON-NUCLEON REACTIONS 



This chapter presents a overview of pion production in nucleon-nucleon reactions. First, 
Section [3.1| motivates the study of pion production and explains some features of the reverse 
reaction, deuteron breakup. In Section [3. 2| the thematic large threshold momentum is intro- 



duced and various kinematical relations are worked out. Next, in Section 3.3 a discussion 
is devoted to the various reaction channels which are relevant for this work. Section 13.41 



presents the hybrid formalism that is used in subsequent chapters. Finally, in Section 3.5 
a brief summary is given of the historical calculations which are being improved upon. 



3.1 General Remarks 

In this section, an outline is given of some general features of pion production reactions, 
which have been studied for many years j40j . Much of this information is drawn from the 
modern reviews found in Refs. [29^ I^Tj I42j . Since pions are the lightest hadrons, pion 
production is the inelasticity of lowest energy of the A^A^ interaction and is interesting for 
several reasons. As explored in Ref. [l3] and more recently in Ref. [M], p-wave pion 
production can be used to study three-body nuclear forces. This connection results from a 
ttNNNN contact term which contributes to a variety of quite different physical phenomena. 
Secondly, pion production is a gateway to understanding production of more exotic mesons. 
One example is the pp — )• pYK (where y is a hyperon and K a kaon) reaction which can be 
used to extract information about the hyperon-nucleon potential. Similarly, A^A^ — ?• NNrj, 
can be used to study the A^A^ — t- A^*A^ transition due to the pronounced role of the A^*(1535) 
as an intermediate state. Other intermediate states can be studied; as will be discussed, the 
A(1232) makes a large contribution to thep-wave A^A^ — )• dir amplitude. Another motivation 
which is explored extensively in Ch. [4] is the determination of Charge Symmetry Breaking 
observables, such as the light quark mass difference. One final motivation that should be 
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mentioned is the usefulness of pion production as a laboratory for the development of xPT 
as the EFT for the strong interaction. The relative momentum of the initial-state nucleons 
required in order to produce a meson is large, and this fact provides a challenge to xPT 
(which is an expansion about small relative momenta) ; pion production could be considered 
the "frontier" of xPT. The continued development of xPT is important in order to better 
understand nuclear phenomena such as nucleosynthesis, fusion, etc. 

Before describing pion production reactions, a brief discussion is given of the reverse 
reaction, deuteron breakup ird — t- A^A^. The principle of detailed balance relates the two 
cross sections 



'^d-K+^pp ^~^^PP^iin+ ■> (3-1) 



where p is the nucleon momentum and q is the pion momentum in the center-of-mass system. 
When considering ird scattering, one sees that deuteron breakup is a primary contribution to 
the imaginary part of the scattering length, aj^d, a quantity which is described theoretically 
most recently by Ref. [46J . Expressing the optical theorem in terms of the scattering length 
a, one finds 



Ima^-rf = lini ( • (3.2) 

g->0 \ 47r / 



By correcting for the other possible final-states (NN^, NNe~^e~ , and NNtt) one can at 
least approximate a relation between Im a^^ and Im flTrd-^TVAf • The most often quoted modern 
experiment [U], extracts the following value from the study of pionic deuterium, 

a^_^m^ = -0.0261(5) + i0.0063(7), (3.3) 

where ttZtt = 139.6 MeV is the charged pion mass. This connection with pionic deuterium 
will be revisited in Ch. |4| where a similar experiments' |48j determination of the total cross 
section for nn — )• dir^ is used. One final important feature of ird scattering is the prominent 
contribution of the A(1232) resonance, clearly seen from the energy dependence of the total 
cross section shown in Fig. |3.1[ This feature emphasizes the necessity of including the 
resonance as an explicit degree of freedom in the EFT. Having gleaned useful information 
from deuteron breakup, let us now return to pion production reactions. 
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Figure 3.1: Average total cross section |(cr^+rf + c^-d) of pion-deuteron scattering (taken 
from Ref. |l9]). The (un)filled circles show experimental data with two different target 
lengths and demonstrate the dominance of A intermediate state. The dashed curve shows 
the naive prediction obtained by summing the contributions of free vrp and ttti scattering, 
while the solid curve shows the multiple scattering calculation of Ref. ^0] . 
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3.2 Kinematics ofNN dvr 

Several of the important features of NN — t- dn can be understood in terms of the large 
relative momentum required in the initial-state. The deuteron binding energy (defined by 
Ed = ran + — Eb) is only E^ = 2.2 MeV, implying that virtually all of the energy of 
the produced pion is coming from the kinetic energy of the colliding nucleons (mjv ~ 940 
MeV). Let the first nucleon have momentum pi and the second p2, such that in the center 
of mass frame pi = —p2- The total momentum of the A^A^ system is P = pi + p2 = 
and the relative momentum is p = (pi — P2)/2 = pi- Working in the non-relativistic 
limit and neglecting the small binding energy, one finds that even at thresholcj^ the relative 
momentum is large, 

= m-TT =^ p = Jm-i^mN = 356 MeV. (3.4) 

rriN 

This momentum sets the scale for pion production reactions. Because this scale lies between 
the "small" pion mass and the "large" chiral symmetry breaking scale A^^, ~ m-Tv, one can 
already see that the convergence of HBxPT is going to be slow. 

In addition to the slow convergence of the EFT, there are at least three other ramifica- 
tions of the large threshold momentum. Firstly, when two particles collide at relatively large 
velocities, they are bound to approach each other closely. This means that pion production 
reactions probe the short-range (r ~ 1/p ~ 0.5 fm) part of the nucleon-nucleon interaction, 
which is completely understood neither in HB^PT nor in nonperturbative QCD. In fact, 
this difficulty requires one to use phenomenology for the initial-state interaction. Of course, 
full coupled-channel calculations exist (for a modern example, see Ref. [51]); however, they 
suffer from separate difficulties, most notably the requirement of chiral symmetry. More 



will be said about the use of phenomenology in Sec. 3.4, where the "hybrid" formalism is 
described. 

Secondly, and quite importantly, note that the wave function of the initial-state A^A^ 
pair is a somewhat distorted plane wave with momentum p while the wave function of 



^ "Threshold" refers to the reaction performed at just the right energy to create a pion with with no 
remaining energy for the pion to have a three-momentum. Clearly, this reaction cannot be measured in 
the laboratory, but is instead studied via extrapolation. 
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Figure 3.2: Momentum mismatch for (a) one-body and (b) two-body production operators, 
illustrated with the real-space wave functions (figure taken from Ref. [41j). Note that this 
figure considers an unbound final-state with real, but small, momentum; the deuteron case 
manifests the same mismatch. 



the deuteron in momentum space contains much lower Fourier components. This leads 
to a prominent momentum mismatch between the initial- and final-states for a one-body 
production operator. In order to obtain a large overlap between these two states, a two-body 
production operator (primarily a meson exchange current, as will be seen) which transfers 
the large three-momentum is required. This concept is illustrated in Fig. |3.2[ 

Thirdly, the fact that the initial-state contains the large scale, p, means that the energy 
E dependence of the cross section will only depend mildly upon the dynamics of the initial- 
state, according to {E/p)^ [52j. The same is true of the production operator, but not of the 
small-momentum final-state nucleons that interact quite strongly and therefore influence 
the cross section energy dependence [53J. Since the focus of this work is near threshold 
kinematics, more will not be said about this point. 

Let us work out two more kinematical relations before proceeding into a discussion of 
selection rules. It is most straightforward for a theorist to work in the center of mass frame, 
but pion production experiments are performed in the lab frame, where a beam of neutron^ 
with energy En = run + 7l {Tl is called the "lab kinetic energy" ) is incident on a stationary 

^The pp — >■ dTT"*" reaction is obviously easier to perform, but since the np — >■ d7r" reaction is related by 
chiral symmetry and avoids Coulombic effects, it is chosen for this works' calculations. 
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proton target. The invariants are written, 

SL,i = (-E„(pi) + nip, pi)^ = (m„ + TTip)^ + 2mpTL 



(3.5) 



and equating them yields a relation between the center of mass relative momentum and the 
lab kinetic energy, 

^2 ^ mlTLi2mn + TL) ^ tunTl 
{ran + "Zp)2 + 2mpTi 2 

In the final state, in the center of mass frame, there is a pion with momentum q and a 
deuteron with momentum — q, such that the invariant energy is ^scj = ET^{q) + i?d(q). 
It is common to refer to the pion momentum in terms of the dimensionless parameter 
r] = g/mTrj^ The algebra does not simplify nicely to allow an elucidating expression relating 
q and p. Nevertheless, one sees that at threshold (q = 0) when neglecting the nucleon mass 
difference {nin = mp = miy), 

2 



2Jmj^+p'^ = 2mN + mT,-Eb =^ p'^ = {m-^ - Eb)mN + ^^—^^-j^^ , (3.7) 



which is the relativistic version of Eq. (3.4). 
3.3 Selection Rules 

Though the arguments can be tricky to follow at times, it is crucial to understand the 
selection rules for pion production as well as the features of the various channels. Using only 
angular momentum algebra and the Pauli Principle, one is able to make striking progress 
in understanding more features of the reaction. The additional quantum number of isospin, 
and its conservation dictated by the chiral symmetry of QCD, already provides us with the 
relation app_^^^+ = 2cr„p_j.^^o due to the fact that an np state can be either T = or T = 1 
with equal probability, \np) = (|10) + |00))/\/2j^ It is simplest to begin the discussion with 

•^There are two potential pitfalls here: 1) the frame must be specified for r; to be meaningful (the CM 
frame is used unless stated otherwise), and 2) this definition is ambiguous for three-body final states like 
in pp — >■ ppyr"; the maximum pion momentum 77max is a commonly used parameter in this situation. 

''charge symmetry breaking in the cross section should not be expected to be visible on top of the much 
larger CS amplitudes. 
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a deuteron in the final-state, although a discussion is also given for the pp — )• ppn'^ reaction 



due to its special place in historical calculations (see Section 3.5). 

The deuteron has the following quantum numbers: total angular momentum J = 1, spin 
S = 1, and isospin T = 0. Furthermore, it is a parity even particle, which would lead one to 
assume that it is a ^5*1 statej^ Interestingly, the nuclear force has a tensor component which 
mixes in the other allowed state, ^-Di. The pion is an 5" = 0, T = 1, odd parity state and it 
can be produced with arbitrary orbital angular momentum with respect to the deuteron 1^^. 

First consider s-wave pion production with = 0. The final state is J = 1 and odd 
in parity and, since the strong force conserves both these quantities, the initial state must 
be as well (so Li = 1,3,5,...). Although isospin breaking operators exist, the largest 
contributions to pion production come from isospin conserving transitions where the initial 
state is T = 1. According to the Pauli Principle, since it is anti-symmetric under 1 -f-)- 2 
in space and symmetric in isospin, the initial state must be symmetric in spin with S = 1. 
To satisfy J = 1 with S = 1, the only possible state is ^Pi. So one sees that only a single 
channel is able to contribute to s-wave pion production. 

Next consider p-wave pion production with 1-,^ = 1- The final state can now be J = 0, 1, 2 
and is even parity. That the parity is the opposite of s-wave production tells us that the 
initial state must be = in order to satisfy the Pauli Principle, thus one has L = J. 
The possible initial states are then ^5*0 and ^D2. Higher will only contribute appreciably 
at higher energies than this work is concerned with. A schematic argument for this point, 
which was formally shown by Ref. [52], is now given. 

The cross section is given by the transition rate per unit flux, 

-=-|T|V (3.8) 

V 

where v is the beam velocity, T is the transition matrix, and p is the density of states. An 



exposition of the details of how the cross section is calculated is given in Section 3.4, but for 



now let us just concern ourselves with the energy dependence. The T-matrix is proportional 



^The spectroscopic notation used is Lj, where L = 0, 1, 2, . . . is denoted S, P, D, 
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to the pion wave function which, to leading order, is just a plane wave, 

e^q-x = ^^Y.i^j,{qr)YU^)YC^{i). (3.9) 
l,m 

which for small q in the channel leads to the dependence, T/^ oc q^^ . Then, using that 
the density of states p oc q, one finds that near threshold, 

ai^ « q^'-+\ (3.10) 

validating the claim that only the low channels will contribute appreciably to the cross 
section. This also means that the deuteron breakup cross section diverges (a ~ l/g) at 



threshold, according to Eq. (3.1) 



Finally, consider the reaction pp — )• ppir^ in the Z,r = channel, which is clearly T = 1 
in both the initial and final-states. Again, since the final-state nucleons have minimal 
kinetic energy and are interacting strongly, Lf = will make the largest contribution. 
Antisymmetry then forces the final-state into ^6*0; the initial-state has opposite parity with 
^Pq. The energy dependence of this reaction is more complicated for two reasons: 1) the 
three-body final-state kinematics are complicated and 2) the strong meson exchange current 
is forbidden by isospin conservation enhancing the effects of Ij^ > channels. Further details 
about these issues can be found in Ref. [iT\. 

To develop an even better feel for these reactions, plots have been included of modern 



experimental data for pp — )• dyr"*" in Figs. 3.3 and 3.4 and for pp — )• ppir^ in Figs. 3.5 and 3.6 
3.4 The Hybrid Formalism 

Having now described the general features, kinematics, and selection rules of pion production 
reactions, the next topic to discuss is the methodology by which one includes the important 
initial- and final-state interactions. It should be stressed that the ideal method would be to 
perform a non-perturbative calculation of these interactions using an elastic kernel derived 
from HBxPT. Unfortunately, this is not yet possible since chiral potentials are only reliable 
for p < 250 MeV [59]. Instead, one uses the hybrid approach in which one calculates the 
production operator perturbatively using HB^PT and convolves this operator with wave 
functions that are calculated non-perturbatively from phenomenological potentials that only 
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Figure 3.3: pp ^ dn^ cross section up to r/ ~ 7 (figure taken from Ref. [M]). The 
prominence of tlie A peak is clear. Tliree different experiments are compared witli two 
theories. 
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Figure 3.4: COSY (labelled "this work") pp — t- dn^ cross section up to ~ 0.5 (figure taken 
from Ref. |54|). A deviation from the neutral pion production experiment is highlighted 
with a fit to the latter data. The linear dependence on the pion momentum near threshold 



predicted by Eq. (3.10) is visible. 
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Figure 3.5: CELSIUS (filled circles) pp — ?• pp-K^ cross section as a function of maximum 
pion momentum ??max (figure taken from Ref. The data is compared with the lUCF 



experiment in Fig. 3.6 along with theoretical (heavy meson exchange) predictions from Refs. 

(MUST]. 




0.2 0.4 06 

Figure 3.6: lUCF (filled circles) pp — t- ppir^ cross section as a function of maximum pion 
momentum r/max (figure taken from Ref. [SH])- The data is compared with three older 
experiments, fits (dashed and dot-dashed curves), and the theoretical energy dependence 
due to phase space and the final-state pp interaction (dotted curve) . 
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obey chiral symmetry in their long-range pion tails. This hybrid approach may seem ad-hoc, 
but it has been successfully applied to a variety of other processes in the two-nucleon sector 
(see Ref. pOj for a nice review, and Refs. |6H 162] for other applications). Note that in the 
remainder of this thesis, the "operator" or "kernel" refers to the part of pion production 
which is calculated perturbatively, while the procedure of combining this operator with 
the non-perturbative parts of the matrix element to obtain an observable is referred to as 
"convolving" or "folding". 

To summarize, the operator is calculated as the sum of two-particle-irreducible diagrams 
involving four nucleon lines (two incoming and two outgoing) and one pion line. Then, a 
phenomenological potential is used to calculate the A^A^ scattering wave functions and the 



deuteron bound state wave function (see Appendix A.3). Finally, the operator is convolved 
with the wave functions to obtain the matrix element. As pointed out in Ref. [B^, the 
phenomenological potentials always involve some form of high- momentum cutoff. Therefore, 
in order to retain consistency, it is important to implement this cutoff in the convolution 
integrals as well. More will be said about this at various points throughout the remainder 
of this thesis. 

Let us now be more precise about the implementation of the hybrid formalism by deriving 
an expression for the cross section in the center of mass frame where the pion momentum is 
q. First one must define expressions for the initial- and final-states. To form an interacting 
A^A^ state with total momentum P = pi + P2, one uses a superposition of free particle states 

= / ^*(p)|iV(p + P/2),Af(-p + P/2)), 

where spin and isospin have been ignored for now. The wave function 'ip{p) is obtained 
by solving the Schrodinger equation with the appropriate A^A^ potential. For an overview 
of the standard numerical techniques required to calculate such wave functions, see Ap- 



pendix A.3 In the center of mass frame, P = 0, and if the nucleons forming the deuteron 
have momentum ki^2! then ki + k2 = K = — q. The invariant matrix element is then 

M{N{pi),N{p2) ^ 7r{^),diK)) = J ^^^*^{k)M {p, k,q) iPr.p{p). (3.12) 
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Note that one treats the initial-state as two separate particles, but the deuteron as a single 
particle. The sum of the two-particle-irreducible diagrams calculated in momentum space 
with external momenta q, pi,2, ^1,2 (the operator) is denoted Ai, and is convolved with the 
external wave functions. Also note that the wave functions will include spin and isospin 
kets on which the operator acts. 

Let us now perform the calculation in position space by inserting Fourier representations 
of both wave functions 

^NN{r) = I ^e^P-'-V^^(p). (3.13) 
A Fourier representation of the operator with respect to 1 = k — p is also inserted 

— 3e^'-A^(l,q). (3.14) 



As described in Appendix B.2 , 1 is the momentum that appears in pion production reactions: 
X(p,k,q) ^X(l,q). 



Now one can rewrite Eq. (3.12) 

M = [ (frd^r'(fr"^*^{r")M{v)^np{r') 
d^rij2{r)M (r)V'np(r). 



d^/c d^p ik-r" ~t(k-p)-r -tp-r' 

(27r)3(27r)3 



(3.15) 



With these choices of Fourier representations, the momentum integrals evaluate to delta 
functions which allow evaluation of the spatial integrals. The result is an integral over a 
single spatial variable. 

Next, one expresses the invariant S-matrix element in terms of M 

{'K\q)d{K) I S I p{pi),n{p2)) = 1 - i{2T,f6\q + K - pi - p2)M{puP2 ^ q, K). (3.16) 
In the center of mass frame, the spin-averaged differential cross section is, ^ 



where the four A^A^ spin states have been averaged over and the three spin states of the 
deuteron summed over. The vector part of the delta function gives q = — K and the energy 
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part gives Ei + E2 = ujq + E^. This allows one to perform all but the cI^k integral, 
What remains is to derive expressions for the wave functions and the operator appearing in 



Eq. (3.15); this is done in Section 4.2.2 



There is one remaining subtlety: in using phenomenological NN potentials, one assumes 
that the A^A^ pair is in its center of mass frame. This is not true in the reaction's center 
of mass frame because the (final-state) deuteron has a non-zero total momentum, — q. In 
the frame where the deuteron is at rest, the T matrix is written in terms of the stationary 
potential, V, and the two nucleon propagator, G, using the Lippmann-Schwinger Equation 
T = V + VGT. However, since the T matrix is invariant, it can also be expressed it in 
terms of the boosted potential in the center of mass frame, V , and the boosted two nucleon 
propagator, G' as T = V' + V'G'T 7^ V' + V'G'V' . Inserting the boosted expression of T 
into the stationary one yields 

V' + V'G'V' = V + VG{V' + V'G'V'). (3.19) 

Rearranging this equation, one obtains 

V' = V + V{G - G")V' 

G" = {^V - GV')G' (3.20) 

« G' - GV'G'. 
Thus the boosted potential up to next-to-leading order is 

V' = V + V{G - G')V' 
^V + V{G - G')V. 

Corrections to the leading order approximation, V' = V, will come in at orders higher than 
will be considered in this work. 

3.5 History of Pion Production Calculations 



(3.21) 



Let us begin this section with an overview of experimental results for the total cross section 
oi pp — )• dn^ just above threshold. In this regime, according to Eq. (3.10), the total cross 
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Table 3.1: Experimental total cross section parameters for pp — )• dvr^ (Coulomb corrected) 
and the isospin-equivalent np — t- dvr*'. The parametrization used is a{pp — )■ (ivr^) = a'q+firf'. 



Experiment 


a (/ib) 


/3 (^b) 


np d-K^ (TRIUMF, 1989) [65] 


184 ±5 


781 ± 79 


TT+d^pp (1991) [66] 


174 ±3 


982 ± 38 


pp dTT+ (lUCF, 1996) [67] 


208 ±5 


1220 ± 100 


pp dTT+ (COSY, 1998) [51] 


205 ±9 


791 ± 79 


Pionic deuterium (PSI, 2010) [l8] 


252lfi 


N/A 



section takes the form a = Aq + Bq'^. Thus the total cross section is typically reported in 
terms of two parameters, 

cr{pp — (ivr"*") = ar/ + j37f 

1 (3-22) 

a{np — ^ (ivr^) = - [ar] + (irf) 



Table 3.1 displays the results of the most recent experiments. 

Pion production has been studied theoretically since the 1960s when Koltun and Reitan 
published their seminal paper obtaining a = 146 /ib [68]. This calculation came well 
before the rise of xPT and considered just two amplitude^ the impulse approximation 
(lA, sometimes referred to as "direct production") and pion rescattering (RS). Diagrams 



representing these processes are shown in Fig. 3.7 In the lA amplitude a pion is produced 
by a single nucleon and does not interact with the spectator nucleon. The leading order 
ttNN vertex is proportional to the pion's momentum; therefore, at threshold (where the 
momentum is zero) the lA proceeds through a recoil correction which is suppressed by 
p/rriN- In the RS amplitude, the ttuNN vertex is commonly referred to as the Weinberg- 
Tomozawa (WT) vertex (KR used the ttN scattering phase shifts to fix the strength of the 
WT vertex). 

Although Koltun and Reitan's results are as accurate as should be expected for pp — t- 
®Of course these amplitudes also appear in xPT with additional constraints coming from chiral symmetry. 
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fa) 



(b) 



Figure 3.7: Amplitudes included in Koltun and Reitan's work: (a) impulse approximation 
and (b) pion rescattering. Solid lines represent nucleons and dashed lines pions. 



dvr"'", this was not the case for the pp — >• ppn^ reaction; models found threshold cross sec- 
tions to be 5-10 times smaller than experiments [69j. The clear reason for this situation 
is the fact that the large RS diagram which contributes to pp — t- (ivr^ is isospin-forbidden 
for pp —7- pp-K^, leaving smaller and less understood diagrams to dominate. This discrep- 
ancy with experiment prompted several more model-dependent solutions, among which were 
heavy-meson exchange (in the kernel) \56\ [57] and off-shell effects in the RS diagram |70l [7T] . 
Although individually these solutions seemed to resolve the issue, adopting both simulta- 
neously changed the situation from dramatic under-prediction to dramatic over-prediction. 
Chiral Perturbation Theory should provide a final solution to this problem since it is model 
independent (provided convergence is achieved) , but the first such calculations [72\ 175] only 
worsened agreement with the data. 

The next development was the recognition by Ref. [74] that the large momentum scale 
p = ^rriT^mN necessitated a reordering of the diagrams into a modified power counting 
scheme which expands in powers of x = p/n^N ~ 'm-K/p- Sometimes this new scheme is 
referred to as the "momentum counting scheme" (MCS) in contrast with the "Weinberg" 
scheme that simply counts all momenta as ©(m^). In the MCS, one examines a diagram 
and deduces whether the momentum in a vertex factor will carry the small final-state 
momentum or the large initial-state momentum. To be more explicit about this modified 
power counting scheme, let us consider p-wave pion production from the two diagrams (lA 



and RS) shown in Fig. 3.7 
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Figure 3.8: Irreducible diagrams which represent the true power counting for the impulse 
approximation 



However, before doing this, let us take a brief aside to discuss a point that will become a 
major theme of this thesis. Formally, the lA diagram requires a modification: it must also 
include a meson exchange due to the fact that a nucleon cannot emit a pion and remain 
on-shell. The definition of an irreducible diagram in ^ = 2 HBxPT, as laid out by Refs. 
|751lll|. is that no intermediate states are off shell by < m^jra^. At best, the nucleon could 
be off-shell by before and after the pion emission; this led Ref. [74J to claim that the 

power counting of the lA diagram should be determined by considering the diagrams of 



Fig. 3.8 The question of how to correctly include one pion exchange (OPE) in the lA is 
one that will be revisited several times throughout the following chapters. 

Now let us return to the task of power counting and start by considering the first diagram 



on the RHS of Fig. 3.8 For p-wave production, the pion production vertex gives a q while 
the intermediate nucleon propagator gives a l/m-,^. Counting q ~ ttItt, as dictated by the 
Weinberg counting scheme, one sees that the lA contributes to the term in the expansion. 
The RS diagram for p-wave production (smaller than it is for s-wave production) gets a m-,^ 
from the WT vertex, a 1 /p^ from the pion propagator, and a q from the pion emission. Thus 
the RS diagram should be counted as mT^q/p^. In the traditional Weinberg counting this 
ratio becomes ~ ^ while in the MCS counting, it becomes ~ x^. This example illustrates 
the subtleties that pion production theorists have been attempting to work out in the past 
decade. 

Unfortunately, the first attempts at applying the MCS [76l [TTJ [78] seemed to indicate 



45 



a possible divergence of the power counting expansion; the NLO contributions were not 
always suppressed as they should be. It was found that the modified power counting scheme 
becomes still more subtle when loops are involved, and after the work Refs. \4:3\ I79j. most 
questions seemed to be answered at least qualitatively. It was shown that with correct 
counting, the first set of loops sum to zero for pp — t- ppir^. On the other hand, the NLO 
pp — )• dn^ loops remained finite and the form of the kernel required an ad-hoc counterterm 
in order to render the convolution integral finite IHU] . This final issue was solved by Ref . [3] 
which recognized another reducibility issue with the diagram involving the recoil correction 



to the WT vertex. This work will be reviewed in Section 4.3.2, but the summary is that the 
NLO loops modify the WT vertex by simply multiplying it by a factor of 4/3. This factor 
brings the theoretical calculation into greater agreement with experiment. 

This is the point where this thesis joins the task. The theoretical model-dependent 
ambiguities have been largely resolved by the use of with an appropriately modified 
power counting scheme. The pp — t- ppn^ cross section within appears to be free 

of theoretical inconsistencies, but due to its small size and the physics community's still- 
incomplete knowledge of the relevant LECs, agreement with experimental cross section data 
is still a relevant question. The next task which one should consider for pp — t- dvr"^ is the 
calculation of other pion production observables. 
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Chapter 4 

CHARGE SYMMETRY BREAKING IN THE np dn^ REACTION 

The asymmetry in the angular distribution of np — t- dn^ due to Charge Symmetry Break- 
ing is calculated using Heavy Baryon Chiral Perturbation Theory. Recent developments in 
power counting have proven successful in describing total cross sections, and we apply them 
to the asymmetry calculation. Reducibility in one of the leading order diagrams is examined. 
We compare the updated theory with experimental results for a set of physically reason- 
able CSB parameters and find that the theoretical asymmetry is larger than experiment 
by ~ 2 — 3 experimental uncertainties. This chapter is a modified version of our published 
paper |2]; we have reorganized and rewritten parts in effort to make it flow better with the 
rest of this thesis. 

4.1 Introduction 

In the preceding chapters, we have reviewed the development of (Heavy Baryon) Chiral 
Perturbation Theory and its application to pion production in nucleon-nucleon reactions. 
We have motivated the investigation of pion production by noting how its special feature 
(the large threshold momentum p = y^m^J^n/v = 356 MeV) provides a challenge to HB^PT. 
Once theorists arrive at an ordering scheme for the expansion, error estimates will become 
more reliable and calculations of different and more exotic reactions possible. 

One particularly interesting quantity that pion production can help determine is the 
magnitude of the CSB parameters discussed in Section |2.4[ Charge Symmetry Breaking 
in the nuclear force is by no means a new topic (see Refs. |34| I81j for example), but pion 
production provides a unique opportunity to study CSB with an observable that vanishes 
in the CS limit. Additionally, the algebraic combination of Snifyf and 5mjv that appears 
in the np — t- di:^ asymmetry is distinct from the sum which appears in the nucleon mass 
difference. This distinction allows for a clean connection between a hadronic experiment 
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and the down-up mass difference. 

Recall that CS is broken both by the mass difference between up and down quarks and 
also by electromagnetic effects. Both of these effects enter HB^PT through interactions of 
neutral pions with nucleons [82] . Unfortunately, the lack of a neutral pion beam prohibits a 
direct measurement of these operators. The addition of a second nucleon opens new doors; 
the np — )• d-K^ reaction circumvents the need for a neutral pion beam and also minimizes 
electromagnetic effects in the initial- and final-states. To be more specific, the angular 
distribution of this reaction is symmetric about 90° in the center of mass frame when CS is 
respected. A recent experiment at TRIUMF [4] was able to observe that this distribution 
is asymmetric at the ~ 2a level. 

This report advances previous work in several ways. The authors of Ref. [5] showed that 
a vertex which was thought to be higher order in fact contributes at leading order to s-wave 
pion production. This advance lead to an improved understanding of the total cross section 
of pp — )• di:'^ ; we extend this calculation off threshold for np — )• dn^ and make a comparison 
with the experimental data. Another ingredient to the CSB calculation is the CS p-wave 
amplitude. We point out a problem with the standard A(1232) amplitude and implement 
a cutoff in order to match p-wave data. Finally, we calculate the asymmetry of np — )• dir^, 
including both RS and lA CSB diagrams, bringing the calculations of Refs. [83^ [3] up to 
date. 



This chapter also contains a discussion (alluded to in Section 3.5) about reducibility in 
the lA diagram. It has been assumed that the energy transfer required for pion production 
occurs in the external wave function when one considers the lA diagram. At the same time, 
it has been recognized that OPE must be included in order to perform the MCS power 
counting for the lA. As a first attempt to resolve this discrepancy, we tried introducing 
"wave function corrections" and found them to be negligible for most of the amplitudes. 
The following two chapters will continue in this path toward an eventual solution which is 
quite different than the naive wave function corrections. The eventual solution is recognized 
by considering the way in which 3D wave functions are obtained by integrating over the 
relative energy in a 4D formalism. 



In Section 4.2 we review the kinematics and selection rules specific to the np — )• dvr pro- 



48 



cess. We also form expressions for the initial- and final-states in this section. In Section 4.3 



we present the diagrams with vertices from the CS part of the Lagrangian (Appendix A.2). 
In this section we also provide a review of the power counting/reducibility developed by 
Ref. [5j and discuss its impact on neutral pion production. Next, in Section [4. 4[ we present 



diagrams with vertices from the CSB part of the Lagrangian which contribute to the reac- 



tion. Our results for the asymmetry are given in Section 4.5 In Section 4.6, we compare 
our work with another recent calculation of the asymmetry and discuss the problems we 
encountered along with possible solutions. 

4.2 The np — ;> dvr^ Reaction 
4-2.1 Kinematics and selection rules 

In this section, we build upon Chapter |3] and discuss the specific kinematics used in the 
calculation of this chapter. There are two frames to keep in mind: the center of mass frame 
(C), and the lab frame where the proton is at rest (E for experiment). The experiment of 
interest [i] is np — t- dir^ performed in the E frame. In this frame the invariant is expressed 
as se = {rup + m„)^ + 2mpTL where Tl is the kinetic energy defined by En = + 
Tl. The asymmetry experiment was performed at Tl = 279.5 MeV, slightly above the 
threshold value of Tl = 275.1 MeV. To simplify the formalism, we use the C frame to do 
the calculation. In terms of the pion momentum q and the deuteron mass rrid we have 



y/sc = Y + + + q2 . In this near-threshold regime, the parameter r/ describes 

the excess energy, ry = |q|/m7r. Equating the invariants we find rye = 0.169 (gc = 22.86 
MeV) at the experimental energy. 



Also discussed in Section 3.3 were the selection rules for the isospin related pp — )• dvr"*" 
and np — )• di:^ reactions. It was shown that to produce s-wave pions, a ^Pi initial-state was 
required, and likewise for p-wave pions, a ^5*0 or state. In the context of this calculation 
we need to distinguish these CS, or "strong" amplitudes with Tj = 1 from CSB amplitudes. 
Charge Symmetry Breaking operators transform as vectors under isospin, thus the initial 
neutron-proton state must have Tj = 0. The same arguments that were used for the strong 
operators then show that s-wave pions are produced from ^Pi np pairs while p-wave pions 
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Table 4.1: Channels for the np wave function in np — )• dvr^ 





Strong 


CSB 




'Pi 


'Pi 


I. = 1 


'So, 'D2 


^Si, 3^1, ^D2 



are produced from the coupled channels ^Si and '^Di in addition to 'D2. These conclusions 
are summarized in Table l4Tl 

The observable of interest in the experiment is the forward/backward asymmetry in the 
differential cross-section given by 

A non-zero asymmetry will only be observed when initial-states of opposite parity interfere. 
However, the interference can only occur for states with the same total spin because the 
spin z-components are summed. 



+1/2 +1/2 

'0 ) = 0. (4.2) 



mi=-l/2 m2=-l/2 



2 2 



1 TTI-i + 771-2 ) ( ~ "T-l) ~ "1-2 



2 '2 



Thus for calculating the asymmetry, we are concerned with two terms: (s-wave strong) -(p- 
wave CSB) and (p-wave strong) -(s-wave CSB). 

4-2.2 Initial and final states 



Recall the expressions derived in Section 3.4 for the differential cross section in the context 
of the hybrid approach, 

ail 54 TT^ s p 4 ^-^ 

M= j (i^rrd{r)M{v)^np{r) (4.4) 
For the np — )• dvr^ reaction, it remains to find expressions for the external wave functions. 
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In the absence of interactions the wave function ■0„p(r) is expressed by performing a 
partial wave expansion on an anti-symmetrized wave function of a free proton and a free 
neutron with relative momentum p. First we consider the strong operators where the np 
pair is in an isospin-1 state, 

(r I V„p> = Vt=i^ (e^P-nmi,m2) ® |T,,i,r,,2> - e-''^-^\m2,mi) ® |r,,2,T,,i)) 

1 ^ 1 (4.5) 

= — (e*P-nmi,m2) - e-^P'^ma, mi)) |r = 1, = 0) , 



where Vt=i is the isospin projector. The bra (r | indicates that we are choosing a basis for 
space, but not for spin or isospin. Implicit in the notation is the dependence of (r | V'np) 
on the momentum p, the spin z-components of the two nucleons, mj, and the isospin z- 
components of the two nucleons, T^^i, with the requirement that T^^i + Tz^2 = 0. 

The exponentials are now expanded and the presence of the strong interaction is added 
by changing the spherical Bessel functions, jiipi^) — ^ &^^'"ULj{r) /pr. The ul^j functions and 
the 8l phase shifts are obtained by solving the Schrodinger equation with a phenomenologi- 



cal NN potential (we use Argonne V18 [H]). This technique is described in Appendix A. 3 



Finally, the spherical harmonics are combined with the spin kets to form states with definite 
total angular momentum. The notation for these states is \{SL)J,mj) \ T,Tz). For the 
allowed quantum numbers, we find 

(r I ^Pnpe^+'Lj, T = l))= 47r(i)^e^^^ ""^^ (1/2 mi, 1/2 m2 | S m,) 

pr 

(4.6) 



X 

m. 



{S ms,L mi-ms\ J rui) Y^^l^Jp) (f | {SL)J,mi) |1,0) 



where m<j = mi + m2 and the second Clebsch Gordan coefficient allows us to make the 
sum over m^ = m/ + m^ rather than m^. For the CSB operators, we have T = np wave 
functions and find 

(r I ^Pnp^'+'Lj, T = 0)) = ±47r(i)^e^^^ ^^^^ (1/2 mi, 1/2 m2 | S m^) 

pr 

(4.7) 



X ^ (5 ms,L mi - m^ I J rm) Y^*i^^{p) (f | {SL)J,mi) (g) |0, 0) 

rrii 

where the ± refers to Tz^i = ±1/2. Similar analysis gives the final-state wave function of 
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the deuteron as a function of its polarization, ruf, 

((10)l,m/ I r) + ^ ((12)l,m/ | f)j . (4.8) 

4.3 Strong Amplitudes 

4-3.1 Diagrammatic Expansion 

Before we calculate the effects of CSB, we need to discuss the CS diagrams which lead 
to a calculation of the total cross section in agreement with experiment. As is always 
the case in an EFT, there are a host of possible diagrams one could draw, and predictive 
power is only obtained by following a power counting scheme. In HBxPT, one orders 
contributions in powers of the external momenta divided by the chiral symmetry breaking 
scale. In the np — )• dir^ reaction, both q and p appear as external momenta and we need 
to keep track of both in the power counting. Recall that because the expansion parameter 
X = p/rriN = \/ rriTr/mN = 0.40 is large, the expansion should not be expected to converge 
quickly. 

The Lagrangian is given in Appendix A.2[ The index of a "type i" vertex is given by 



iyi = di + ^-2, (4.9) 

where di is the sum of the number of derivatives, rriTr's, and (5's (the AN mass difference), 
and fi is the number of fermion fields. In standard power counting at tree level, the sum 
of the Vi for each vertex in a diagram indicates the power of x at which that diagram 
contributes. This rule, however, requires modification due to the relatively large value of p 
as discussed in Section [3. 5[ We will come back to this issue again in Section [4.3.2 



There are three two-particle irreducible diagrams which can be drawn using the vertices 



from C^^\ They will be referred to as the impulse approximation, or lA, (Fig. 4.1 i); pion 
rescattering, or RS, (Fig. |4.l| 3); and Delta (Fig. |4.l| ;) diagrams. 



Let us recall a few facts about these three diagrams. We have seen in Section 3.2 that 
the RS diagram is enhanced relative to the lA diagram due to the "momentum mismatch" 
being provided for in the RS but not the lA diagram. Using the leading order vertices, 
the RS diagram contributes to the channel with s-wave pions. Although the Delta diagram 
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a 



ib) 



ic) 



Figure 4.1: Leading order contributions to np — ?• dir^. Solid lines represent nucleons, the 
double solid line represents a A, and dashed lines represent pious. 



Figure 4.2: Kinematically consistent, but topologically reducible impulse contributions. 
Solid lines represent nucleons and dashed lines represent pions. 



provides the momentum transfer required, the Delta resonance is at 1232 MeV and the vrA'^ 
energy is w 1080 MeV, so the Delta diagram is also somewhat suppressed for our situation 
of interest. Finally, we note that both the lA and the Delta diagrams with leading order 
vertices contribute to the channels in which the pion is in a p-wave. 



We saw in Section 3.5 that there is ambiguity as to what is meant by Fig. 4.1 1, and we 
must now make a decision on how it is to be included. On the one hand, we know that a 
single nucleon cannot emit a pion and remain on-shell. But on the other hand, the diagrams 



in Fig. 4.2 (the same diagrams that were used in Fig. 3.8 to power count the lA) which 



remedy this problem by including OPE appear to be two-particle-reducible, topologically. 



Near threshold the energies of the exchanged pions in Fig. 4.2 are cj w mTr/2. However 
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cz:> 



+ 



+ 



+ ... 



Figure 4.3: "Hybrid" approach. Solid lines represent nucleons, dashed lines represent pions, 
dotted lines represent pions with uj = 0, and filled ovals represent A^A^ strong interactions. 



□ 



Figure 4.4: Complete impulse contribution. Solid lines represent nucleons, dashed lines 
represent pions, dotted lines represent "wave function" pions with w = 0, and the square 
represents the operator which is used for the full impulse approximation. 



the diagram of Fig. 4.1 1 can be thought of as being evaluated as the sum of the diagrams 



shown in Fig. 4.2, This is due to hybrid nature of the calculation; once the operator is 
calculated, it is convolved with NN wave functions. One of the major terms of the strong 
interaction potential at low energy arises from static OPE (cj = 0)j^ The effects of static 



OPE are schematically shown in Fig. 4.3 



Thus we make the following naive choice: to obtain the correct impulse contribution, 



we add up the contributions from Fig. 4.1 1 and Fig. 4.2 and then subtract what is already 



included in the wave functions (the last two diagrams of Fig. 4.3). This calculation is 
schematically shown in Fig. |4.4[ 



^We ignore the effects of the rest of the wave function for the moment. 
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The OPE propagator is of Yukawa form, 



DAu,r) = -^—, (4.10) 



where = y^Tn"^ — ui"^ ~ y/S/lmT^. Thus subtracting off the final two diagrams in Fig. 
amounts to making the replacement 



4.4 



(4.11) 



in the exchanged pion propagator. The final four terms of Fig. 4.4 comprise a correction to 
the impulse diagram, which we will refer to as a "wave function correction." We find that 
this correction is ~ 4% of the total impulse amplitude at the experimental energy, and we 
include it in our calculation. However, it is important to note that a different method is 
developed in Chapter [6] to handle the lA within the hybrid formalism. To give a complete 
picture, we leave this attempt as it was at the time this paper was published. 

4-3.2 Power Counting 

Now we will look more closely at the size of these diagrams in the MCS counting scheme. 



we 



using the counting techniques summarized in the review article |41j . In Section 3.5 
explained how to power count the lA and RS diagrams for p-wave production, and in this 
section, we expand to the rest of the diagrams and channels. The first step is to consider 
how the external momenta and masses appear in the calculation of a particular diagram. 
Then, MCS can be used to assign each diagram's status as LO, NLO, etc. 

Let us examine the Feynman rules which we will use. The propagators are calculated 



from the Lagrangian in Eq. (A. 6) 



Dn{p) 



+ ie 



Da{p) = ■ (4-12) 



p2 + (m^ — (p^y) — ie 

Also from the leading order Lagrangian, we see that both the ttNN and the vrA^A vertices 
are proportional to |q|, the pion momentum at that vertex. Note that this momentum is p in 
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the OPE vertices. The only other vertex appearing thus far is the WT, which is proportional 

to Uq^in + '^<j,out- 

The external particles have the same momenta in each diagram. The produced pion 
has q = {ojg, q) ^ {jn-,^, 0), and the incoming nucleons have pi^2 = (-S'1,2, ip) ~ ("i7r/2, ±p) 
in the Heavy Baryon formalism in which the nucleon mass is subtracted off of the energy 
component. 



We saw in Section 3.5 that for the p-wave lA diagram of Fig. |4.2^ the final emission 
contributes q, the nucleon propagator l/m-7r, and the OPE p ■ ■ p so that the whole 
diagram is ~ q/mT^. These same arguments give that s-wave RS diagram is ~ ^^—j^, and 
the p-wave Delta diagram is ~ m^s - Thus, according to the MCS scheme, the p-wave lA 
and Delta diagrams are order and the s-wave RS diagram is order x^- However, more 
practically, we note that q/rrij^ = V ^ and 6 ~ 2m7r so that the lA and Delta diagrams 
are numerically ~ x^- This ordering comes in agreement with the fact that the s-wave 
diagrams should make a larger contribution to the cross section than p-wave diagrams. 



Naively, these three amplitudes from Fig. 4.1 represent the complete LO calculation, x^ 
for s-wave and x^ for p-wave. However, it is well documented that these three amplitudes 
alone do not correctly reproduce the experimental data for the reactioij^ Near threshold 
[r] ~ 0.05), the most recent np — )• dn^ experiment [65] found a = a/r] ~ 90 /ibj^ while for 
these first three diagrams, we find a ~ 55 /jh. 

The solution to this problem was discovered by Ref. [5], who noticed that the = 1 



"recoil" correction to the WT vertex, which is found in Eq. (A. 7) goes like (qin + qout) • 



(Pin + Pout)/(2mAr). For Fig. 4.5i, this vertex (a filled circle) ~ p^/(2m7v) = ?TT.7r/2 and thus 
this diagram is of order the same order as the u = rescattering diagram. Similarly, 
one finds that the s-wave portion of the recoil correction to the impulse diagram (Fig. |4.5[3) 



is of order x^- In this diagram, the filled square represents the sum analogous to Fig. 4.4 for 



the recoil diagram. We find that the wave function corrections are more important (~ 20%) 

^Admittedly, this is not necessarily a serious problem since we have only considered the leading order 
Lagrangian and have already stated the convergence will be slow. A more serious problem is expounded 



upon at the end of this section 

^Note that in this chapter we 
smaller than in the other chapters 



^Note that in this chapter we use a different definition for a; its values in this chapter are two times 
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ia) 



ib) 



Figure 4.5: Recoil corrections. Solid lines represent nucleons, dashed lines represent pions, 
and the filled circle and square represent u = 1 vertices. 



in this case. Finally, the s-wave portion of the Delta diagram's recoil correction is found to 
be higher order and is therefore ignored. 

The recoil corrections to the propagators have also been included in the calculation where 
applicable. For this reaction, the only such diagram is Fig. |4.2| 3 where the 3-momentum in 
the nucleon propagator is large (~ p). For that propagator, we use the corrected version, 

^n{p) = -0 — « - — • (4.13) 



Using this propagator rather than the u = version doubles the size of Fig. \A.2}p . Never- 
theless, this diagram (minus its = analog) is already very small. Thus the net effect 
of correcting the propagators is small for this reaction at this order. 

Including all these recoil corrections (especially Fig. |4.5^) brings the theoretical cross 



section near the experimental results as shown by the solid curve of Fig. 4.6 Due to the 
relative scatter of the data shown in Fig. |4.6| it is difficult to tell how well the theory is 
reproducing the experiment. Regardless, it is clear that theoretical improvement has been 
made. 

It should also be noted that the subtlety of reducibility and recoil corrections in this 
reaction resolves questions about NLO loop diagrams discovered by Ref. [80j. Namely, the 



sum of all the NLO irreducible loops in Fig. 4.7 is found to be proportional to p. This is a 
problem because such sensitivity of the operator to the A^A^ wave function is not physical. 
The solution to the problem is, once again, to consider including OPE in the operator, 



57 
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0.15 



0.20 



Figure 4.6: Cross section for np — )■ dn^ in terms of a = cr/r/ as a function of the pion center 
of mass momentum, rj = q/m^o. Circles with error bars display the experimental results of 



Ref. [65 1 . The dashed line displays the results of including the diagrams in Fig. 4.1 and the 



solid line displays the results of also including the recoil terms discussed in the text. 



(a) 



(b) 



/ 

/ 



id) 



Figure 4.7: Irreducible loops. Solid lines represent nucleons and dashed lines represent 
pions. 
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ia) 



ib) 



Figure 4.8: Reducible loops. Solid lines represent nucleons and dashed lines represent pions. 



this time for the rescattering diagram. There are two resulting diagrams shown in Fig. 4.8 
It was shown in Ref. [5] that in these topologically reducible loops, the recoil corrections 
to the nucleon propagators need to be included in addition to the WT's recoil correction. 
This key argument is an application of the MCS power counting scheme. Reference [Sj 
then showed that (part of) the energy dependence of the WT vertex "cancels" one of the 



nucleon propagators leaving a reducible diagram similar to Fig. 4.73,. This diagram is equal 
in magnitude and opposite in sign to the aforementioned NLO sum, resolving the issue. 

The other term that remains from the original loop integral after this manipulation 
is still of reducible form but now has an on-shell WT vertex ~ Im-,^. This term would 
already appear upon convolution of the rescattering diagrams discussed above (including 
recoil corrections) with external wave functions, i.e. this term is truly reducible. This result 
can be stated another way: cancellation of the irreducible loops comes from a short range 
piece of including OPE in the rescattering operator. Thus a complete NLO calculation 
must include this OPE. However, its remaining reducible piece would only be included in 
the calculation if the recoil correction to the WT vertex is included in the rescattering 
operator. Therefore it is both consistent and necessary to include the WT recoil correction 
in the pion production operator. 

4.3.3 P-Wave Observables 

Another important test of the theory is how well it describes p-wave pions [44]. This is 
especially important for the asymmetry, which involves strong p-waves at leading order. 
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The differential cross section can be expanded in Legendre polynomials, 

da 



Qo + aiPi(cos(6')) + a2P2(cos(6i)) + 



(4.14) 



where 9 is the angle between p and q. Note that the total cross section plotted in Fig. 4.6 



is a = ATTao/rj. As discussed in Appendix B.3 a2 receives contributions almost exclusively 
from p-wave pions. The ratio 02/0^0 is therefore used as a test for this part of the theory. 



We find that the diagrams of Fig. 4.1 along with their recoil corrections overestimate the 
data by approximately a factor of two. 

Upon closer inspection we find that the ^Sq amplitude (which is supposed to be small) 
is relatively large. This amplitude is coming mainly from the Delta diagram, as can be seen 



in Appendix B.3 where the values of the reduced matrix elements are listed. To remedy the 
situation in the simplest way possible, we implement a cutoff for the Delta diagram, 

-i -i / A2_\ 

(4.15) 



p2 + /i2 



—i —i 



A2 



,2 ■ 



^p2+/i2 p2+A2yA2-/x- 

One can show that doing this essentially softens the OPE potential for r < log(A//x)/A. 
Note that one consequence of such a cutoff is that it modifies both the and the 
channels. Clearly this cutoff is not an acceptable long-term solution for an effective field 
theory, but the fact that such a procedure is necessary is interesting given that the reaction 
occurs at an energy ~ 150 MeV below the Delta resonance. That p-wave pion production 
is highly sensitive to the strength of its contact term, 



ct 



dl 

fn 



nUt -a- Vtt I N 



(4.16) 



was discussed by Ref. [33]. Also note that if we were to instead use a dipole cutoff, we 
would find a larger value is needed, A — )• ^/2K. The vrA^A vertex should fall at least as fast 
as the jNN vertex which, according to "quark counting", goes as 1/Q^. In this chapter, 



we use the cutoff in Eq. (4.15) for simplicity, but recognize that the corresponding A will 
be small. 

In a complete calculation, the cutoff dependence of the offending ^Sq amplitude will 
be absorbed by the contact term, di = di{A). Naive dimensional analysis (NDA) of this 
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0.01 



A (GeV) 



Figure 4.9: Cutoff dependence of the Delta diagram's ^5*0 amplitude. The vertical axis shows 
the change in the dimensionless radial integral / which depends on the cutoff according to 



Eq. (4.15). This change is defined according to SI (A) = I{oo) — /(A). 



term leads one to expect di 1/A^. Other tree- level diagrams will contribute to the ^5*0 
amplitude at this order; however, for the purposes of this thesis, we will assume that the 
Delta diagram is large enough to justify ignoring these other diagrams. In this simplified 
case, we can investigate the validity of the NDA scaling by looking at the A-dependence of 
the Delta amplitude. 



The log-log plot of Fig. 4.9 shows the change in the -M-^g^ amplitude (the corresponding 
radial integral is what appears in the plot) as a function of the cutoff. The true value of di 
is set by requiring that some experimental p-wave observable is exactly reproduced; for this 
simplified analysis we use the condition, 

-^^50 =-^^50+ -^4 =0- (4-17) 
For large values of the cutoff, the amplitude is unmodified and the contact term takes on a 



large magnitude in order to satisfy Eq. (4.17). For small values of the cutoff, the amplitude 



is suppressed and the contact term is smaller in magnitude. The calculated scaling behavior, 
given by the slope of the log-log plot, is ~ —2, indicating that 6Ai^g^{A) 1/A^. 

There are two ways to interpret the 1/A^ scaling. Firstly, higher order contact terms 
are not required since the calculated scaling is not dependent on a larger negative power of 
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Figure 4.10: Legendre coefficients of tlie differential cross section for different values of 
the cutoff. Data is from an np — )• dir^ experiment (Ref. |65j, circles) and an pp ^ dTr~^ 
experiment (Ref. [67j, squares) in which the data have not be corrected for Coulomb effects. 



the cutoff. Secondly, it appears that the contact term contributes at lower order than NDA 
would lead one to believe. This discrepancy is either due to the cutoff dependence of the 



other diagrams which contribute to the ^Sq amplitude (and were ignored in Fig. 4.9), or 
else a modified power counting scheme will be required. In such a modified scheme, the di 
contact term would be promoted to leading order. 

In fact, there already exists a modified power counting scheme for spin-triplet A^A^ scat- 
tering at low energies. The original identification of this breakdown of Weinberg counting 
was given in Ref. [85]. As discussed recently in Ref. [86], the relevant small scale is the 
strength of OPE in the triplet channel, specifically, the tensor interaction parametrized by, 

A. = ^, (4.18) 

which is numerically ~ Im.,^. In pion production, the A A — )• A A transition that partic- 
ipates in A^igjj (of which the final state is spin-triplet) proceeds exclusively through this 
tensor interaction, so it is not surprising that we observe this scaling behavior. Further 
discussion of this issue is beyond the scope of this thesis. 



Figure 4.10 shows the effects of this cutoff on the ratio a2/ao, where A = 10 GeV 
represents the original theory (such a large cutoff has no significant effect) . The amplitudes 
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-0.30 ^- 



Figure 4.11: Analyzing power for different values of the cutoff. Data is from a pp ^ dn'^ 
experiment (Ref. [67J, circles) in which the data have not be corrected for Coulomb effects. 



for np — )• dir^ are related to those for pp dTr~^ (which are bigger by \/2) when charge 
independence is respected. Thus, the ratio plotted should have the same value for both 
reactions. By adjusting the cutoff to fit the data, we find A = 310 MeV. 

Another useful observable for testing p-wave pion production is the analyzing power. 
Ay, which is defined 

^^^'^ = da,{e) + dam ^'-''^ 
'^^^'^^^^"eii^i ^ \M{m,,y = ±1/2, e)\\ (4.20) 

where mi^y = ±1/2 refers to the fact that the beam is polarized perpendicular to the 
scattering plane. In the z-basis, these states are 

= ±1/2) = = = (4.21) 

v2 



As shown in Appendix B.3, Ay is proportional to the product of s-wave and p-wave am- 
plitudes. Figure |4.11| shows the effects of the cutoff on this observable. Again, charge 
independence implies that Ay should be the same for both neutral and charged pion pro- 
duction. We find the best agreement with the analyzing power data for A = 560 MeV. 
Below, we will display our results using both the original theory and a cutoff taken to be 
the geometric mean of two fits, A = 417 MeV. 
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Figure 4.12: Leading CSB contribution. Solid lines represent nucleons, dashed lines repre- 
sent pions, and crosses represent 1^ = 1 CSB vertices. 



4.4 Charge Symmetry Breaking Amplitudes 

4- 4-1 Leading- Order CSB Amplitudes 

The fact that the up and down quarks have different mass is reflected in the Lagrangian 
by including terms which break chiral symmetry [87J. The leading such terms are given in 



Eq. (A.7) and have coupling constants 5mN and 5mn which are constrained by 

5mN + SniN = nin — rup. (4.22) 



Recall from Section 2.4 that the 5mN term has its origins in the quark mass difference and its 
size is ~ [md — mu] = e{md + mu) with e 1/3. Chiral symmetry tells us that (m^ + mu) oc 
m^, and so dimensional analysis along with the hadronic scale misj yields drriN ~ ern^jm^. 
The 5m AT term is of electromagnetic origins, but is of the same order as bm^. These CSB 
operators appear in the rescattering diagram depicted in Fig. |4.12 where the CSB vertex 



is denoted with a cross. The size of this diagram is 5mN/p €m'^/{miyp) = ex^- Note 
that although the full nucleon mass difference appears explicitly in the Lagrangian at this 
order, the corresponding operator [NW-^N) does not change the parity and thus does not 
contribute to an asymmetry. 

The Cottingham formula |88j (a model for the electromagnetic contribution to m„ — nip) 



can be used to obtain dmj^f = —(0.76 it 0.30) MeV. The constraint of Eq. (4.22) then fixes 



5mN = 2.05 ± 0.30 MeV. In this chapter, we use this value of 6mN for our final results, 



64 



though other values have been suggested. Reference |3j discusses models which predict 
values for 6m]\j leading to 1.83 < ^ 2.83. Additionally, a recent LQCD study [89] 

found = 2.26 it 0.72, having used the MILC collaboration's determination of rriu/md 

m. 



Though we do not present the full details of the calculation of this CSB diagram, let us 
briefly consider on the isospin algebra which one encounters. 



M oc [<5m^(T(l) • t(2) + r(l)3T(2)3) - <5m^(T(l) • t(2) - r(l)3r(2); 



(4.23) 



Evaluating this between T = states yields. 



M oc 



drriN 



(4.24) 



This is the distinct algebraic combination we referred to in the introduction. Since this 



combination is different than that appearing in Eq. ( 4.22 ) , one can use Eq. ( 4.22 ) to eliminate 



5mis[ in favor of the well-known nucleon mass difference. 

Let us now discuss another LO source of CSB coming from a more detailed evaluation 



of the strong RS diagram shown in Fig. 4.1 3 which was made by Ref. |91]. We became 
aware of this group's efforts when nearing the end of our study and worked with the authors 
of Ref. [91J to resolve differences between our calculations. For clarity in this thesis, we 
include this new CSB amplitude into the main body of the calculation. 

Charge symmetry breaking occurs in the strong RS diagram as a consequence of the time 
derivative in the WT vertex. The diagram is typically considered by treating the external 
nucleons as identical. The fact that they are not identical means that the energy transferred 
in the pion exchange will not only include the m7r/2 required by the kinematics, but also a 
term ib(m„ — mp)/2 (positive if the WT vertex is on the proton, negative if on the neutron). 
The effect of this contribution is shown in Ref. [91] to be equivalent to a change in the CSB 
rescattering diagram, 



SrriN 



5mN 



5mj\f 



5mj\f 



+ 



5mN + ^rriN 'idrriN 



(4.25) 
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Figure 4.13: ly = 2 CSB contributions. Solid lines represent nucleons, double solid lines 
represent A's, dashed lines represent pions, crosses represent v = 2 CSB vertices, and the 
boxed cross represents the full impulse CSB diagram including OPE. 



4-4-!^ Next-To-Leading-Order CSB Amplitudes 



Another CSB term given in Eq. (A. 8) involves one derivative and one ~ em'^/m'^ 



and is thus a u = 2 vertex with momentum dependence |q|. This vertex appears in the 

er/x^- In Fig. 4.13 1, the boxed cross 



diagrams of Fig. 4.13 whose sizes are f^iq/m-j^ 



represents the sum analogous to Fig. |4.4| for the CSB impulse diagram. Again, the wave 
function corrections are small (2%). 

Little is known about /3i, the impulse CSB coupling. As a starting point, Ref. [9^ notes 
that this CSB operator can be viewed as arising from n — rj mixing. The result shown is 
that 

,2 ■ 



/3i 



e mz 



m 



N 



(4.26) 



with the idea that the value Crj ~ 0(1) is "natural". As discussed in the review [32], 



0.10 < ^ < 0.51. 
Air 



(4.27) 



Also, Ref. [93] gives (7r°|F|r/) = -0.0039 GeV^ and we use = 547.51 MeV. These values 
result in —0.47 < < —0.21. Thus it is at least plausible that the /3i term could originate 
naturally from r] — ir mixing. We note that the rj' could also give such a term, but do not 
consider it here. 
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Using Eqs. ( |4.26[ ) and ( |4.27[ ), we obtain -3.2 x lO"^ < /J^ < _i.4 x lO'^. Note that the 
value used in the original calculation of the asymmetry by Ref. [83] was /3i = —8.7 x 10^^, 
which we refer to as the "extreme value". However, according to the above discussion, 

2 

the "natural" size is f3i ~ —'^^ ~ —6 x 10^^. Thus even though its origins may not lie 
exclusively with the r], the aforementioned "extreme" value for /3i is not extreme at all from 
the effective field theory's point of view. 

4.5 Asymmetry Results 

Contributions to the asymmetry come from interference terms between the p-wave part 
of the strong amplitude and the s-wave part of the CSB amplitude, and vice versa. The 
issue is somewhat complicated because, in contrast to threshold emission, each diagram can 
contribute in both the s-wave and the p-wave. However, contributions to the sub-leading 
parity (s-wave for the impulse and Delta diagrams and p-wave for the rescattering diagrams) 
are formally higher order. An example of this is the strong rescattering diagram for p-wave 



pions considered in Section 3.5 which was found to contribute according to ~ rjx 



The contributions to the asymmetry are depicted in Fig. 4.14 Figure 4.14 a, includes 
strong p-waves and CSB s-waves and has size (r/) x (ex^)- Figure 4.14 3 includes strong s- 
waves and CSB p- waves and has size (x) x Figure 



4.14 



includes strong p- waves and 



CSB s-waves and has size (rjx^) x i^X^)- Thus we find that in these kinematics Fig. |4.14 a. 
(~ e^X'^) is the LO contribution, and Fig. 4.14 a, c (~ e??x^) both come in at NLO|^ Other 
interference terms involving these diagrams are higher order, ~ e??^X^ or smaller. Finally 
we note that this work is not intended to be a complete NLO calculation since loops and 
higher order vertices may contribute at this order. The identification of all NLO diagrams 
is beyond the scope of this work. 

The details of the calculation are spelled out in Appendix |B} The coupling of the WT 
vertex (and its recoil correction) is determined by chiral symmetry, but the LECs /tt, qa, 
and hA appear in the calculation and we use /,r = 91.9 MeV, gA = 1.267, and Ha = 2.1gA- 



We use the following masses: niiy = {nin + mp)/2 = 938.919 MeV, rur, 



rrir. 



1.293 MeV, 



''Here we refer to the ~ ^vx'' terms as NLO even though they are suppressed by because no terms 
suppressed by a appear. 



67 



□ 



X 



{b) 



^ - - X 



X 



X 



Figure 4.14: Interference terms for the asymmetry in np — t- dn^. Solid lines represent 
nucleons, double solid lines represent A's, and dashed lines represent pions. The filled circle 
and square represent v = 1 strong vertices and crosses represent CSB vertices. The boxes 



represent full impulse diagrams in the sense of Fig. 4.4 
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rriA = 1232 MeV, and m.j^ = m^o = 134.977 MeV. In Appendix B.l the amplitudes are 
expressed in terms of reduced matrix elements based on the partial wave decomposition. 



Next, in Appendix B.2, we present an example calculation: the reduced matrix elements 



corresponding to the RS diagram. Finally, in Appendix B.3, we connect the observables 
(cross section and p-wave observables) with the reduced matrix elements. Also in this section 
one can find a table of our results for each individual reduced matrix element. 

In Table [4?2| we display our asymmetry results. Because the asymmetry is linear in the 
CSB amplitudes (and therefore the CSB parameters), we are able to present our results as 
a set of coefficients, {x, y} defined by 

^/.xl0^ = x.(^)+y(/3ixl03). (4.28) 



The primary advance made in this work over the previous calculation [3] is the inclusion 
of the rescattering and impulse recoil corrections; this improvement is shown in moving 
from the top four rows to the next four rows. At LO in the asymmetry calculation this 
simply increases qq, but at NLO it affects both the numerator and the denominator of 
the asymmetry. The final result, using the set of physically reasonable parameters from 



Section 4.4, is shown in the last column. Our most accurate calculation (the third to last 
row) uses a cutoff to suppress the ^5*0 amplitude and includes the NLO interference terms 
to find, 

= 40.2 X 10"^. (4.29) 

The effects of using different values for the cutoff can be seen in the last two rows of the table. 
Due to the large number of theoretical issues which still need to be addressed (primarily the 
cutoff in the strong p- waves and the model dependence of the CSB couplings), we believe 
that providing an uncertainty at this time would be false advertising. 

The experiment of Ref. [4j found Ajb = [17.2 ± 9.7] x 10"'' and thus our calculation 
overestimates the data by approximately 2.5o"^^p*. It is interesting to note that using 



the most "extreme" set of parameters discussed in Section 4.4 (Suin = 1-83 MeV, /3i 



8.7 X 10 '^) in the cutoff NLO calculation yields Afb = 26.5 x 10 Clearly the current 
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Table 4.2: Asymmetry in np — )■ dvr^ as a function of CSB parameters 5mN and "LO" 



and "NLO" represent the sums discussed in Fig. 4.14 The Delta cutoff is of monopole form; 
a dipole cutoff would be \/2 bigger. 



Afb X 10'' = X 




X 10^) 




A/b(2.05 MeV, -0.0032) x 10^ 


Order 


Delta Cutoff 


X 


y 


LO (no recoil) 


None 


33.7 





69.1 


LO (no recoil) 


A = 417 MeV 


27.6 





56.6 


NLO (no recoil) 


None 


37.6 


1.4 


72.6 


NLO (no recoil) 


A = 417 MeV 


32.5 


1.8 


61.0 


LO 


None 


25.0 





51.2 


LO 


A = 417 MeV 


18.9 





38.9 


NLO 


None 


28.1 


1.4 


53.1 


NLO 


A = 417 MeV 


22.1 


1.6 


40.2 


NLO 


A = 310 MeV 


20.1 


1.7 


36.0 


NLO 


A = 560 MeV 


24.0 


1.6 


44.3 
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calculation cannot be fully judged until more accurate values for the CSB couplings are 
known. 

4. 6 Discussion 

4-6.1 Comparison with other calculations 

The first calculation of the np — )• dn^ asymmetry used a A^A coupled channel formalism and 
included the CSB impulse vertex as well as other, smaller effects arising directly from the 
neutron-proton mass difference |83]. This study reported Ajh = —28 x 10~^. The second 
calculation included only the CSB rescattering vertex, and found Afj, = 60 x 10^^ [3j. Both 
these calculations were preformed before the work of Ref. [5] which brought the total cross 
section into agreement with experiment. Our work brings the asymmetry calculation up to 
date. 

We have seen that there is reason for concern regarding the theoretical description of 
p-wave pions, which comprise the entire strong contribution to the LO asymmetry. Because 
the total cross section is dominated by the RS diagram, small changes to the p-wave am- 
plitudes are able to significantly modify the asymmetry while only slightly changing the 
total cross section. As a temporary solution, we implemented a cutoff in the Delta diagram 
and thereby achieved acceptable agreement with the p-wave data. Another solution to this 
problem is to use a coupled-channel A^A potential for the initial-state. This approach was 
taken by Ref. [M] who were able to achieve good fits to these data without a cutoff, since 
the OPE of the Delta diagram is then part of the wave function. 

Another difference between our calculation and that of Ref. [91] is that in order to 
improve on the theoretical uncertainty (~ 30%) of the Legendre coefficient ao, they used 
experimental data (from pionic deuterium) to obtain oq = 1.93 /xb. This is significantly 
larger than the theoretical NLO value we use, ao = 1.49 /xb (oq = 1.28 /ub for A = 417 MeV), 
and leads to a smaller value for the for the asymmetry. Note that experiments for neutral 
[65] and charged [67J pion production found ao = 1-39 /ib and ao = 1-64 fih (Coulomb 
corrected) respectively. The final difference between our calculation and that of Ref. [911 is 
that they do not include any of the NLO asymmetry terms. 
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For the sake of comparison, we used our code to calculate the LO-only asymmetry 
with A = 417 MeV, using qq = 1.93, and using Ref. |9l]'s quoted values for qa and Jt^. 
These choices should reduce any differences between our calculations to those resulting from 
different wave functions. For these choices we obtain Ajb = 14.0^ x lO"'^, which is to be 
compared with their result of Afh = 11. ^ 10~^. Although they did not present it this 
way, one can use the Cottingham sum rule along with their result to obtain Aj^ = 23 x 10~^, 
which only overestimates the data by less than one experimental uncertainty. 



4.6.2 Outlook 

Several issues remain to be understood theoretically. Firstly, it appears that a contact term 
will be required to suppress the ^Sq channel in the strong amplitude if one uses a purely A^A'^ 
initial-state. The interesting physics observation here is that the Delta part of the A^A^ wave 
function seems to be much more active than it is A^A^ scattering. This contact term comes 
in at NNLO in the p-wave calculation, so it is more important at this time to come to a 
definitive conclusion on the lA reducibility issue. The wave function corrections introduced 



in Section 4.3.1 will be investigated in a more complete manner in the next chapter. 

Secondly, the difference between the more recent experimental determination of ao (from 
pionic deuterium) and the older np — )• dir^ data (which agrees with NLO theory) plays a 
large role in our over- prediction of the asymmetry. This situation becomes even worse when 
the cutoff is used to decrease the p-wave amplitudes. For these reasons, we conclude that 
further calculations are necessary. In particular, one should extend the calculation to next 
order while examining both the power counting of recoil terms and the reducibility of loops. 

More generally, the existence of multiple mass scales greatly complicates the power count- 
ing for this reaction, and it is clear that a converging expansion cannot yet be definitively 
claimed. Another interesting aspect of this calculation is the use of a hybrid formalism. 
One can argue that using phenomenological potentials to determine the A^A^ wave func- 
tions is equivalent to working to all orders in the EFT. Thus there is a mismatch when 
the calculation of the operator is truncated at some order. One way to remedy this situa- 
tion (introduced by [94j ) is to use a cutoff when calculating the Fourier transforms of the 
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operators. This approach is taken in the next chapters. 
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Chapter 5 

IMPULSE APPROXIMATION IN THE np dn^ REACTION 

REEXAMINED 



The impulse approximation (one-body operator) in the np — )■ dn^ reaction is reexamined 
with emphasis on the issues of reducibility and recoil corrections. An inconsistency when 
one pion exchange is included in the production operator is demonstrated and then resolved 
via the introduction of "wave function corrections" which nearly vanish for static nucleon 
propagators. Inclusion of the recoil corrections to the nucleon propagators is found to change 
the magnitude and sign of the impulse production amplitude, worsening agreement with the 
experimental cross section by ~ 30%. A cutoff is used to account for the phenomenological 
nature of the external wave functions, and is found to have a significant impact for A < 2.5 
GeV. This chapter is a modified version of our published paper [6]; we have reorganized and 
rewritten parts in effort to make it flow better with the rest of this thesis. 

5.1 Introduction 

As we have seen, there are two significant difficulties involved in pion production calcu- 
lations. Firstly, the large threshold momentum p ~ ^Jrf^^m^ necessitates a reordering, 
sometimes called the MCS, of the xPT expansion. Secondly, the identification of diagrams 



as either reducible or irreducible cannot be determined strictly by topology. In Section 4.3.2 
we reviewed the work of Ref. [5j , which addressed these two difficulties in the NLO pp — t- dvr^ 
calculation. Pulling a pion exchange from the external wave functions into the term of the 
kernel containing the recoil TTirNN vertex resulted in an irreducible diagram. This diagram 
cancelled other NLO loops that had an unphysical sensitivity to the wave functions in the 
convolution integral; the remainder of the new diagram was shown (in Ref. [5]) to effec- 
tively enhance the LO ttitNN vertex by a factor of 4/3. In this paper we perform a similar 
investigation of reducibility in the impulse approximation. 

An attempt to address this issue was put forth in our study (Chapter |4| published 
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version [2j) of charge symmetry breaking in np — )• dir^ where we introduced "wave function 
corrections". These corrections were calculated for the final state, and found to be a small 
fraction of the impulse diagram that they are correcting. However, this calculation suffers 
from a particular approximation which we will describe. Fixing this approximation has 
a very significant effect. Furthermore, we show that the wave function corrections in the 
initial state are larger than one would expect in the MCS scheme^ 



In Section 5.2 we review the np — t- dvr reaction and the impulse approximation's role. 



Then, Section 5.3 examines the inconsistency that is found when one includes static OPE 
with the impulse approximation. Also in this section, wave function corrections are pre- 



sented as a solution to the problem. Section 5.4 discusses the correct implementation of the 
recoil corrections to the nucleon propagators. The effects of including a cutoff are shown 



in Section 5.5, Finally, we discuss the total cross section in Section 5.6 and conclude in 
Section 15.71 



5.2 Pion Production Review 

The pion production operator in momentum space depicted in Fig. |5.1| is a function of the 
pion momentum q and 1 = k — p, where k (p) refers to the final (initial) relative momentum 
of the nucleons. The momentum transfer between the two nucleons in Ai is defined as 
q' = P2 — k2- We use a phenomenological, non-relativistic potential (V) that is static: the 
energy of each individual nucleon is conserved by V. Given this choice, working in the 
center of mass frame requires q'^ = iOq/2. It should be noted that this is an approximation 
(the "fixed kinematics approximation" ) that one needs to adopt in order to work in position 
space. If one works in momentum space and fixes q' ^ via energy conservation at the NNn 
vertex, this is called the "equation of motion approximation." For the RS diagram, Ref. 
|95| found that both of these approximations have problems, particularly when considering 
initial state interactions. Nevertheless, to work within the hybrid formalism one of the 
approximations is required, and it appears that the former is preferable to the latter. 

^The topic of this chapter is (once again) taken up in Chapter [g] where we present a final definitive 
conclusion. For completeness, we are leaving the arguments made in this paper (nearly) as in their 
published form. 
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Figure 5.1: Pion production operator. Solid lines represent nucleons, dashed lines represent 
pions, and ovals represent interactions. 
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Figure 5.2: Feynman rules 



In this work we employ threshold kinematics, where q = (m7r,0), q' = (m,r/2,l), 
Pi,2 = ("171/2, ±p), and ki^2 = (0, ±k) with |p| = 359 MeV. Also, at threshold only s- 
wave pions {In = with respect to the deuteron) are produced and the initial state is purely 
^Pi (see Section 



3.3 



for the details of this partial wave decomposition). We use the hybrid 



methodology where "operators" (two-particle irreducible diagrams) are calculated perturba- 



tively (as in Appendix B.2) and then convolved (as in Section 3.4) with A^A^ wave functions 



which are obtained using phenomenological potentials (as in Appendix A.3). Appendix A. 2 



discusses the chiral Lagrangian that defines the theory. Using Eqs. (A.6) and (A.7), we 



obtain the Feynman rules shown in Fig. 5.2 
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Figure 5.3: Impulse approximation operator alone (a), with OPE in the final state (b), and 
with OPE in the initial state (c). Solid lines represent nucleons, dashed lines represent 
pions, and ovals represent the full NN T-matrix. 



At leading order, ©(x^jj^the s-wave amplitude is dominated by the rescattering dia- 
gram, where a single pion is emitted from one nucleon and inelastically scattered by the 
other nucleon into an on-shell-produced pion. The only other leading order s-wave diagram 
(recall that the Delta diagram is p-wave at threshold) is the recoil impulse approximation. 



shown in Fig. |5.3[ a). In Appendix C.l we present the details of how the lA diagram is tra- 
ditionally calculated; that is, ignoring the reducibility issues. For the initial and final states, 
we make use of three different potentials: Argonne vl8 |84j . Nijmegen II [96], and Reid '93 



96] . The results for the s-wave reduced matrix element Aq (see Appendix B.l) are shown 



in Table 5.1 along with the LO rescattering results. These numbers have been calculated 
many times before; nothing new has been considered yet. Note that the rescattering results 



in Table 5.1 include the factor of 4/3 (which truly belongs at LO in the expansion) that 



we discussed in Section 5.1 Finally, we point out that the experimental data (Section 5.6) 
imply a reduced matrix element of 80 < Aq < 94. 

Recall our previous discussions of the fact that the impulse approximation operator of 
Fig. |5.3[ a) cannot formally be convolved with the initial and final states as described above 
because the nucleon emitting the pion cannot remain on-shell. To put it another way, q' ^ 

^Recall that the expansion parameter is defined, x = = ^ ~ ^/rn^irrriN. 
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Table 5.1: Reduced matrix elements of the rescattering and impulse production operators 
for three different potentials. 



Diagram 


Avl8 


Nijm II 


Reid '93 


Ares 
^0 


76.9 
4.9 


83.4 
1.3 


80.3 
3.5 



vanishes in this diagram because there is no way for energy to be transferred. The common 
approximation made in pion production calculations is to ignore this formal difficulty. Also 
recall the prescription of Ref. [74] , which pulled an OPE from the final state wave function 
(Fig. |5.3[ b)) in order to argue that the impulse approximation was leading order. This can 
be seen as an application of the Lippmann-Schwinger equation for the final state deuteron, 

li'd) =GV\iPd), (5.1) 

where G represents the two-nucleon propagator and V represents the full potential. In Figs. 
|5.3[ b) and |5.3| [^c), we have replaced V with one pion exchange. This replacement is known 
to be a good approximation for the deuteron [29], but is not as valid for the initial state. 
Including contributions to V other one pion exchange is beyond the scope of this work. 



Equation (5.1) begs the question: are the diagrams in Figs. |5.3| [^a) and |5.3| [^b) the same 
size? If so, we will be able to conclude that neglecting the required energy transfer in Fig. 
|5.3[a) does not significantly alter the true magnitude of the lA amplitude. 



5.3 Including One Pion Exchange 



Calculation of Fig. |5.3[ b) using Eq. (5.1) is detailed in Appendix C.2, For lack of a 
better name, we will call this the "OPE reducible" diagram. In this calculation, we take 
G = {E — Ho)~^ = {—Ed — p^/m^v)""'^. The energy of the exchanged pion in this case is 
taken to be g"' = 0. This choice is consistent with the fact that the OPE is the ffi'st term 



in the V of Eq. (5.1), which should be the same V that is used to generate the initial and 



final wave functions. The results are shown in the first row of Table 15.21 
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Table 5.2: Reduced matrix elements of the impulse approximation with OPE in the final 
state [see Fig. |5.3[ b)]. 



Diagram 


Avl8 


Nijm II 


Reid '93 


xOPE,red,f 
^0 

4 0PE,irr,f 
^0 

/.OPE,irr,f .OPE,red,f 
^0 ^0 


75.2 
75.6 
0.5 


64.6 
64.7 
0.1 


79.3 
79.8 
0.5 



We find an inconsistency between the impulse approximation (Fig. 5.3 'a)) and OPE 
reducible (Fig. |5.3| ^b)) diagrams: although they are equivalent according to the Lippmann- 
Schwinger equation, they are of very different size numerically. Using Avl8, they are 4.9 
and 75.2, respectively. Of course this inconsistency is not surprising when one notes that 
three-momentum transfer is provided for in the latter diagram but not the former. 

To resolve this problem, we reconsider the diagram in Fig. |5.3[ b) as a fully irreducible 
operator, that is a member of the kernel prior to convolution with external wave functions. 
This is justifiable in that the left intermediate nucleon is off-shell by m7r/2, more than the 
m^/mjv typical of reducible diagrams. If we view it in this way, we are free to chose the 



energy of the exchanged pion to be g"^ = 171^^/2 as mentioned in Section 5.2 The single- 
nucleon propagator for the left intermediate nucleon is taken from the rules shown in Fig. 



5.2 The calculation of the reduced matrix element for this "OPE irreducible" operator is 



detailed in Appendix C.2 and the results are shown in the second row of Table 5.2 We find 
that this diagram, which correctly accounts for energy transfer, is approximately equal to 
the OPE reducible diagram. 

The question remains: should Fig. |5.3| ^b) be included, and if so, how? Until a clear 
procedure is defined for going from the full four-dimensional ttNN coupled-channels for- 
malism to the more common three-dimensional uncoupled formalism, this question is open 
to interpretation. We continue to take the view proposed in Chapter |4j which is that the 
OPE irreducible diagram should be included with the OPE reducible diagram subtracted 



off to prevent double counting. This difference is shown in the third row of Table 5.2, and 
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Table 5.3: Reduced matrix elements of the impulse approximation with OPE in the initial 
state (see Fig. |5.3| |^c)). 



Diagram 


Avl8 


Nijm II 


Reid '93 


/,OPE,red,i 
^0 


-11.2 


-23.7 


-15.0 


<OPE,irr,i 
^0 


-11.2 


-23.7 


-15.0 


<OPE,irr,i .OPE,red,i 
^0 ^0 


~0 


~0 


~0 



is referred to as the (nearly vanishing) "wave function correction" . Note that the numerical 
cancellation is not as trivial as it appears. Schematically, OPE is two derivatives acting 
on a Yukawa e~™''/mr that has different ranges for the reducible and irreducible cases. 
In going from reducible to irreducible, the radial integral gets bigger because the range 
increases. However, the derivatives bring down inverse powers of the range such that the 
overall amplitudes are similar in size. 

Let us now discuss a new issue. For s-wave pions at threshold we have a cr ■ (pj + p/) 
at the vertex where the pion is produced. For this reason Chapter [4] only considered OPE 
in the final state, assuming initial state OPE (Fig. |5.3[ c)) to be suppressed by the small 
final state momentum k. However, one needs to be careful when applying power counting 
to calculations that involve external NN wave functions. At small distances the momenta 
of the nucleons (derivatives in position space) are distorted away from their constant values 
at asymptotically large distances. As an example of this difficulty, it can be shown that 
at short enough distances the k^/2m7v operator becomes larger than the p^/2m7v operator 
in the context of the RS diagram. For this reason, we also calculate Fig. |5.3[c) (for the 



details, see Appendix C.2). The results of this calculation are shown in Table |5.3| where, 
again, the full wave function correction is nearly zero. 

There is another formal point to discuss with regard to the above calculations that was 
left out of Chapter [4j Expressions for nucleon propagators in irreducible diagrams differ 



based on the power counting scheme used. Consider the situation shown in Fig. 5.4 , Starting 
from the full, relativistic propagator, the authors of Ref. |9^ showed that, if ~ m-,^ and 
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P-q 



P = mj^v + p 



Figure 5.4: Momenta of the nucleon propagator where v = (1,0,0,0) is used 



p ~ yJmT^niN, the correct propagator after emitting the pion is 

-gO + (2p-q-q2)/2mAr + ie' ^ ' ' 

where the second term in the denominator is the recoil correction, suppressed by one power 
of x'li ^ ^JrrQm^. Note that we have ignored the resulting vertex corrections and anti- 



nucleon effects, both of which are suppressed for any choice of cP . Equation (5.2), which we 
will refer to as the "new" method, comes in opposition to the "old" method which derives 
the propagator from the non-relativistic chiral Lagrangian, 

(5.3) 



pO _ gO _ (•p _ q)2y'277),^ _|_ if' 



where the third term comes from the NLO Lagrangian and is a candidate for promotion in 
MCS counting. Although we will not present the argument of Ref. |97) in full, /2mj\j ~ 



m7r/2 is clearly the connection between Eq. (5.2) and Eq. (5.3) 



Let us now consider the external pion vertex of Fig. |5.3[b), where in terms of the 



momenta shown in Fig. 5.4 we have p = (m-,^ /2,'p) and q = (m^,0). Since p — g ~ m,, 



we promote the recoil corrections in the old propagator and find 

ZLt„ 



■Yirr,f ^ 

'new _ 



^^irr,f ^ 

m7r/2 — rriTT — p^/2mAr 



(5.4) 



Next, consider the left-side OPE vertex of Fig. |5.3[c) where we have p = (m,r/2,p) and 
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q = {—uj, p — k), and thus 

iG 



new 



w + (p2 - k2)/2mAr + ze 

(5-5) 

^ ~ + - k2 /2mAr + ie ' 

Note that in the absence of distortions (|p| ^/rf^/rnN, |k| ~ 0), iG^cw = ^Goid for both 
the initial and final state propagators. For the sake of clarity we will define as the "free 
recoil approximation" (FRA) the use of these free particle values for the nucleon momenta. 

5.4 Nucleon Propagator Recoil 



In Section 5.3, the FRA was used for the recoil corrections to the nucleon propagators. In 
this section we calculate the diagrams again, treating the momenta properly as operators 
instead of numbers. In this section we use the old nucleon propagators for the irreducible 



diagrams. In doing so, we avoid the G^^^ of Eq. (5.5), which would be difficult to evaluate 
exactly in position space. It was pointed out in Ref. |98] that the old nucleon propagators 
have formal convergence problems owing to the large external momenta. Nevertheless, we 
expect to gain insight into the validity of the FRA using these propagators, 

^^irr,f ^ 

- p2/2mAr' 

(5-6) 

iG'"'' = . 

TTiTT — k^/2mAr + ie 

Of course, according to MCS, k should not be counted as ~ y/n^/nTN and the k2/2mjv 



term should therefore not appear as in Eq. (5.6) until higher order. We choose to retain it 



here as an investigation into the effects of the distortions. For the reducible diagrams we 
continue to use G"^^^ = {E — Hq)~^. 

The matrix elements can be calculated exactly in position space with Green function 



methods (see Appendix C.3 ). The results are shown in Table 5.4 We find that the final state 
wave function correction evaluated without the FRA gives a sizeable negative contribution 
of approximately —10. Additionally, we find that the initial state corrections become as 
important as the lower-order final state corrections. To verify the surprising results of this 
calculation, we examine as an example the ttzat — )• oo limit of the radial integral for the 
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Table 5.4: Reduced matrix elements for the wave function corrections with proper treatment 
of the momenta using the old expressions for the nucleon propagators. 





Avl8 


Nijm II 


Reid '93 


4 0PE,irr,f 
^0 

4 0PE,red,f 
^0 

^OPE,irr,f .OPE,red,f 
^0 ^0 


80.8 
92.4 
-11.7 


70.8 
81.2 
-10.4 


89.1 

103.3 

-14.2 


.OPE,irr,i 
^0 


5.1+24.5i 


15.1+34.7i 


7.8+27.8i 


^OPE,red,i 
^0 


16.2+8.2i 


22.9+9.9? 


18.2+8.7i 



^OPE,irr,i _ ^OPE,red,i ^n.i + iG.Si -7.7+24.8i -10.4+19.H 



irreducible initial state OPE in comparison with its analog from the previous section (which 
is independent of mj\[ since k = is used), 

I^Jdrr' (^|:^ + ^) G'°P'^'"^'H2/(m./2,r) +g(mj2,r))i?,,(r), 

(5.7) 



where the functions / and g are defined in Appendix C.2 As shown in Fig. 5.5, in the 
large-mTv limit the recoil term in the propagator vanishes and we recover the leading order 
result. 

5.5 Cutoff Dependence 

It should not come as a surprise that odd things are happening in the short-distance part of 
the wave function, especially when we take derivatives. In the hybrid formalism we are using, 
this domain of the wave function is calculated from a phenomenological potential: Woods- 
Saxon for Avl8, one boson exchange for Nijm II, and Yukawa for Reid '93 (the very short 
range is exponential, exponential, and dipole, respectively). Because these potentials are 
fitted to experimental phase shifts, the wave functions derived from them can be considered 
as infinitely high order in the EFT. Thus one should consider using a cutoff to account for 
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-10 ./ 



Figure 5.5: Irreducible initial state OPE integral as a function of tjin using AvlB. The solid 
line displays the FRA result and the real (imaginary) part of the exact propagator result is 
shown as a dashed (dotted) curve. 



this mismatch between the operator and the wave functions. Use of such a cutoff is referred 
to as EFT*, and was introduced in Ref. |63j . 

In this section we investigate the effects of cutting off the convolution integrals that 



account for the the presence of initial and final states as discussed in Section 5.4 We use 
the procedure of Ref. [94J , which modifies the Fourier transforms with a Gaussian cutoff, 



M{r) = I ^e^'-5i(l2)A^(I) 



{51 



Note that the traditional one-body lA is not affected by such a cutoff scheme. For the OPE 
operators we define g\{uj,r): 



47r ^ J (2^)3 l2 + ^(a;)2- ^ ^ 

The exact evaluation of this integral and of the derivatives required to compute the diagrams 



of this work are shown in Appendix C.4 As desired, the cutoff regulates the behavior of 



g{r) at the origin, as shown in Fig. 5.6 The cutoff dependence of various reduced matrix 



elements is shown in Fig. 5.7 



The fact that we observe significant cutoff dependence of the wave function corrections 
above the typical scale of ~1 GeV is surprising. Indeed, this sensitivity indicates the need 
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Figure 5.6: Comparison of g (solid curve) and (^a (dashed curve) with uj = and 
A = 1 GeV. 
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Figure 5.7: Cutoff dependence of various reduced matrix elements for AvlS (solid curve), 
Nijmll (dashed curve), and Reid '93 (dotted curve). 
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Table 5.5: Reduced matrix elements for three different potentials. 



Diagram 


Avl8 


Nij 


ni II 


Reid '93 


A = oo 


A = 1 GeV 


A = oo 


A = 1 GeV 


A = oo 


A = 1 GeV 


Rescattering (NLO) 


76.9 


76.2 


83.4 


81.5 


80.3 


79.1 


Impulse 


4.9 


4.9 


1.3 


1.3 


3.5 


3.5 


Final wfn cor (FRA) 


0.5 


-0.3 


0.1 


-0.7 


0.5 


-0.5 


Final wfn cor (exact) 


-11.7 


0.4 


-10.4 


-1.5 


-14.2 


-0.7 


Initial wfn cor (FRA) 


^ 


0.1 


^ 


0.1 


^ 


0.1 


Initial wfn cor (exact) 


-11. 1+16. 3i 


-13.0+11. 7i 


-7.7+24.8'i 


-13.0+14.2-i 


-10.4+19. 1'i 


-13.0+13.0i 


Total (FRA) 


82.2 


80.8 


84.8 


82.2 


84.2 


82.2 


Total (exact) 


59.0+16.3'i 


68.5+11.7'i 


66.7+24.8i 


68.4+14.2i 


59.2+19. 1'i 


69.0+13.0i 



for a counterterm because observables must be cutoff independent. As pointed out in Ref. 
|80j . if one considers the difference of terms that comprise the wave function correction, 

1 1 ml/A 1 

^' (5.10) 



q/2 _|_ q/2 _|_ 377j2 j 4 q/2 _|_ ^2 q/2 _|_ 377^2 ' 

it can be argued that wave function correction is N^LO in the MCS scheme. However, if 
this view is to be accepted, the fact that the wave function corrections are much larger in 
magnitude than the LO impulse approximation should be considered surprising. 

5.6 Cross Section Results 



Shown in Table 5.5 is a summary of the findings discussed in this paper at A = 00 and 



A = 1 GeV along with the rescattering diagram. 



We also discuss the total cross section; as we recall from Section 3.2 near threshold the 
total cross section is parametrized, 

a = ^ (ar? + /37?3) , (5.11) 
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a (fib) 




(a) FRA nucleoli propagators 

Figure 5.8: Cutoff dependence of the total 
curve), and Reid '93 (dotted curve). 



a (lib) 




(b) Exact nucleon propagators 



section. AvlS (solid curve), Nijmll (dashed 



where q = rriT^r]. At threshold, one can only calculate a, 

where s = {rriT^ + rrid)'^. Also recall that charged pion production is related to neutral pion 
production by isospin symmetry. This symmetry is the reason for the 1/2 present in the 



definition of a. The most recent experimental data were shown in Table 3.1 



The theoretical total cross section as a function of the cutoff is shown in Fig. 5.8 
The theoretical results include all diagrams up to 0{x^). Thus theory can assign a rough 
uncertainty to the threshold cross section of 2 x ~ 30%. 

5. 7 Discussion 

Before the findings of this work, the total theoretical cross section at A 1 GeV [Fig. |5.8(a)| 



was in agreement with the most recent experiment (fifth row of Table 3.1 ) at approximately 



the Icr level. Fixing the FRA approximation decreases the cross section, and if we stop 



here (Fig. 5.8(b)), the agreement between theory and experiment becomes more tenuous 
(approximately the 2.5cj level). A second conclusion of this work regards the MCS power 
counting, which dictates that |p| ~ ^Jn^^rriN while |k| ~ m-,^. We find that the wave function 
corrections of the initial state are of similar size to those of the final state once the FRA 
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is removed, contradicting the previous sentence. There exists a contact term (an NNNNn 
vertex) at N^LO along with tree-level diagrams proportional to the Cj LECs and two-pion 
exchange loops. Since all the LECs except the contact term are fixed by other data, it will 
be interesting to see if that contact term is of natural size. 

The next chapter investigates further into the concept of reducibility. Specifically, we 
will define a clear procedure for deciding what to include in the impulse diagram. There has 
been a lack of consensus in the literature as to the inclusion of OPE, and if it is included, 
how that should be done. Understanding this issue is important, not only for calculation of 
the total cross section but also for p-wave pion production, since the leading contribution 
to p-wave pion production comes from the impulse diagram. 
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Chapter 6 

IMPULSE APPROXIMATION IN NUCLEAR PION PRODUCTION 
REACTIONS: ABSENCE OF A ONE-BODY OPERATOR 

The impulse approximation of pion production reactions is studied a final time by de- 
veloping a relativistic formalism, consistent with that used to define the nucleon-nucleon 
potential. For plane wave initial states we find that the usual one-body (IB) expression 
OiB is replaced by O2B = —iK{mT^/2)OiB/miT, where K{mT^/2) is the sum of all irre- 
ducible contributions to nucleon-nucleon scattering with energy transfer of 171-^/2. We show 
that O2B ~ OiB for plane wave initial states. For distorted waves, we find that the usual 
operator is replaced with a sum of two-body operators that is well approximated by the op- 
erator O2B ■ Our new formalism solves the (traditionally ignored) problem of energy transfer 
forbidding a one-body impulse operator. Using a purely one pion exchange deuteron, the 
net result is that the impulse amplitude for np — )• (Itt^ at threshold is enhanced by a fac- 
tor of approximately two. This amplitude is added to the larger "rescattering" amplitude 
and, although experimental data remain in disagreement, the theoretical prediction of the 
threshold cross section is brought closer to (and in agreement with) the most recent data. 
This chapter is a modified version of our paper [7] ; we have reorganized and rewritten parts 
in effort to make it flow better with the rest of this thesis. 

6.1 Introduction 

In this chapter we take up once again the calculation of NN — )• dn with the pion in an 
s-wave. We continue to focus on the contribution of the impulse approximation in which 
the produced pion does not interact at all with the spectator nucleon. Recall that pion 
rescattering, not the lA, is known to make the largest contribution to the total cross section 
|68| . The A(1232) resonance is also known to contribute significantly to this observable, 
but in the p-wave channels. Our motivation for the present study is to obtain increased 
precision in the total cross section calculation and to prepare for future application to other 
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Pl P2 



Figure 6.1: Impulse approximation without initial state interactions. Solid lines represent 
nucleons, dashed lines represent pions, and double solid lines represent a deuteron. 



observables to which the lA contributes, such as p-wave pion production. 

As we have seen, a challenge in the calculation of pion production is the presence of 
strongly interacting initial/final states. When considering the lA diagram it is at first 
unclear whether OPE should be included in the kernel or not; historically, calculations have 
not done this, claiming instead that the necessary energy exchange is happening in the 
external wave functions. This method was brought under question in Chapter [5] (see also 
Ref . [80j ) . Ideally, one would like to derive the correct method from a relativistic formalism 
that cleanly separates effects in wave functions from those appearing in the kernel. 

Consider the lA contribution to NN — )• dir in the plane wave (PW) approximation where 



initial state interactions are neglected (see Fig. 6.1). The amplitude for such a process has 
been estimated to go like Ai^"^ ~ ^^cr •pi(/)(p) ~ ^^y^m/vm^0(p), where (f){p) is the bound 
state wave function, evaluated in momentum space. The suppression by m^'^^ was noted 
in Ref. [71], which also included an analysis that a more detailed treatment of the power 
counting based on including initial and final state interactions introduces a power of 
via an energy denominator such that the amplitude varies as ■sjm^. Nevertheless, we see 
directly an explicit times <\){^^frnPm^^ . 

In the physical region where m-^ = 140 MeV, the wave function falls as a power of 
momentum greater than unity. For small values of relative momentum, the deuteron wave 
function also falls more rapidly than an inverse power of its argument. If one takes to 
be small, the deuteron remains weakly bound |99^ llOOj and therefore its momentum wave 
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function will also fall rapidly in the chiral limit. Thus the power counting can only be 
considered a very rough estimate. If we follow |74j . the impulse term is a leading order 
term, but the deuteron wave function is quite small for physical values of p and there 
is also a substantial cancellation between the deuteron s- and d-states. Thus this term's 
contribution to the cross section [68] is small and there is a contradiction between power 
counting expectations and realistic calculations. 

This contradiction was also discussed at length in Chapter [5] where we introduced "wave 
function corrections" as a possible solution. This proposal included OPE with an energy 
transfer of 171-,^/ 2 in the impulse kernel, but then subtracted off a similar diagram with static 
OPE in order to prevent double counting. The result depended strongly on the treatment 
of the intermediate off-shell nucleon propagator and no definitive conclusion was reached. 
This chapter is intended to settle the debate regarding the inclusion of OPE in the impulse 
approximation. We demonstrate, by starting from a consistent relativistic formalism, that 
non-static OPE is to be included with no subtraction necessary; the impulse amplitude that 



should be used is given in Eq. (6.28). Furthermore, we show that the traditional approach 



of using a one-body kernel is correct only in the absence of initial state interactions. 

The key issue in understanding the energy transfer mechanism of the lA is the reduction 



of the hybrid formalism from four to three dimensions. In Section 6.2 we review aspects of the 



four-dimensional Bethe-Salpeter (BS) formalism for the two-nucleon problem. Section 6.3 



presents the N — t- Ntt operator and Section 6.4 shows that for plane wave NN — t- dir, 
the traditional lA is approximately valid. Next, Section 6.5 considers the full distorted- 
wave amplitude by calculating the corresponding loop diagram, including the effects of the 
non-zero time components of the momenta of the exchanged mesons. We are then able to 
interpret the resulting three-dimensional distorted-wave amplitude as a sum of two-body 
operators. We demonstrate the new formalism by explicitly evaluating s-wave A^A^ — t- 
dvr amplitudes at threshold. To aid the flow of the arguments, approximations made in 



this section are verified to be sub- leading in Appendices D.3 to D.5 A comparison with 



experimental cross section data is made in Section 6.6, where we also discuss implications 
and future directions. 
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6.2 Bethe-Salpeter Basics 

Recall the definition of the nucleon-nucleon potential from the Bethe-Salpeter formalism. 
We follow the approach of Partovi and Lomon |101j and also consider the relationship 
between the Bethe-Salpeter wave function and the usual equal time wave function as recently 
discussed in Ref. |102| . 

Partovi and Lomon write the Bethe-Salpeter equation for the nucleon-nucleon scattering 
amplitude 7W as 

M = K + KGM, (6.1) 

where K is the sum of all irreducible diagrams. The quantities Ai and K depend on the 
total four-momentum Ptot and the relative four-momentum k. The two individual momenta 
are pi^2 = -ftot/2 ± k and G is the product of two Feynman propagators: 

^=(-^ '—^] (-^ '—^] =^1^2. (6.2) 

The quantities M and K differ from those of |1U1| by a factor of — i/(27r). Partovi and Lomon 
replace the relativistic G by the Lippmann-Schwinger propagator g for two particles. For 
scalar particles, g is obtained from G by integrating over the zero'th (energy) component of 
one of the two particles |lU2j . For fermions, one must also project onto the positive energy 
sub-space of both particles. This is accomplished in the center of mass frame by taking 

m 



i?(k)(p2,-4m2 -4k2 + 



It] 



where i?(k) = yk^ -|-m^. Note that g contains the important two-nucleon unitary cut. 
The non-relativistic potential U is defined so as to reproduce the correct on-shell NN 
scattering amplitude M using the Lippmann-Schwinger (LS) equation 

M = U + UgM. (6.4) 



The quantity U is obtained by equating the Ai of Eq. (6.1) with that of Eq. (6.4) to find 

m 

U = K + K{G - g)U. (6.5) 



92 









K 









Figure 6.2: Bethe-Salpeter equation near the deuteron pole. 



In solving Eq. (6.4) for the on-energy shell scattering amplitude, U never changes the 



value of the relative energy away from 0. Equations (6.4) and (6.5) are consistent with 
Weinberg power counting in which one calculates the potential using chiral perturbation 
theory and then solves the LS equation to all orders. The term G — g may be thought of 
a purely relativistic effect arising from off-shell (short-lived) intermediate nucleons, and in 
the present context a perturbative effect. 

Consider the deuteron wave function in the final state of a pion production reaction. 



For P near the pole position, the second term of Eq. (6.1) dominates and we replace the 



scattering amplitude with the vertex function F: — )• F, and 

F = KGT. 



(6.6) 



This equation is shown pictorially in Fig. 6.2 The Bethe-Salpeter wave function ^ is 
defined as GT so that 



= GF = GK^. 



(6.7) 



The wave functions of the scattering state and the deuteron are shown in Fig. 6.3, If one 



uses Eq. (6.4), the bound-state wave function (f) is obtained by solving the equation 

= 5f/0 = gUgU(t>. 



(6.8) 



The second equation shows that U also is evaluated at vanishing values of time components 
of the relative momenta. We will treat the amplitudes \I',(/> and the state vectors |\I'), \ (t)) 
(either bras or kets) as interchangeable. 
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(a) A''A'^ scattering (b) Deuteron 

Figure 6.3: Bethe-Salpeter wave functions. 



The next step is to relate ^ with which can be thought of as the usual bound-state 
wave function. This is most easily accomplished by using the projection operator P on the 
product space of two positive-energy on-mass-shell nucleons. We then have 



PG = GP = Gp = g, 



(6.9) 



with the last step resulting from the explicit appearance of two positive-energy projection 



operators for on-mass-shell nucleons in Eq. (6.3). We define Q = I — P and use the notation 



= P'i/, = and PKP = Kpp, PKQ = Kpq , etc. The Q-space includes all terms 
with one or both nucleons off the mass-shell. The amplitude ^p contains the ordinary 
nucleonic degrees of freedom so one expects that it corresponds to (j). This is now shown 



explicitly. Use I = P + Q in Eq. (6.7) and multiply by P and then also by Q to obtain the 



coupled-channel version of the relativistic bound state equation: 

*P = GpKpp^p + GpKpQ^Q 

= GqKqp^p + GqKqq^q. 



Solving Eq. (6.11) for ^'g and using the result in Eq. (6.10) gives 

^Q = [l-GQKQQ]-'GQKQp^fp 

*P = Gp (Kpp + Kpq[G-q^ - KqqY^Kqp) ^p, 



(6.10) 
(6.11) 

(6.12) 
(6.13) 



but one can multiply Eq. (6.5 ) hy P ■ ■ ■ P etc. to obtain the result 



Upp = Kpp + Kpq[Gq - Kqq] ^Kqp, 



(6.14) 
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thus Eq. (6.13) can be re-expressed as 



xjTp = GpUpp^p = gll^p. 



This last equation is identical to Eq. (6.7). Thus we have the result that 



(6.15) 



(6.16) 



^'p is not the complete wave function, but we expect that is a perturbative correction 
because the deuteron is basically a non-relativistic system. 

6.3 The N ^ Ntt Amplitude 

We now turn to the application of the Bethe-Salpeter formalism to the problem of threshold 
pion production. First, we remind the reader of the one-body pion production operator 



in HBxPT, which was derived in Section 2.3.2 and Appendix A. 2 Here we repeat the 
derivation without making (yet) the non-relativistic reduction. The nucleon field is split 
into its heavy (H^) and light (N^) components. 



^-(x) = e-™-" {N,{x) + H,{x)) 



(6.17) 



P_^(x) 



where P± = (1 it ^)/2 and v is the velocity vector satisfying = 1 and chosen in this 
work to be v = (1,0). The heavy component is integrated out of the path integral and the 
resulting free equation of motion for the light component has a solution, 



N{x) = 




—i(E—m,pf)t+ip-x 



D.l 



(6.18) 



we show 



where -E = y p ^ + and x is a two-component Pauli spinor. In Appendix 
that the LO Feynman rule for the s-wave N — t- A^vr amplitude vanishes at threshold and 
that the NLO rule is 



. m-^ QA 5 ^\ 



(6.19) 



where the derivatives act on the nucleon wave functions. 



95 



6.4 The NN dir Reaction: Plane Wave Initial States 



Traditionally, the impulse approximation to pion production is calculated by using the 



operator of Eq. (6.19) as the irreducible kernel to be evaluated between non-relativistic 



nucleon-nucleon wave functions for the initial and final states. Between two-component 
nucleon spinors, 7^747'^ — <7i, so 



-i- znrcr ■ 



\Pl,P2) 



(6.20) 



2mN 2/„ 

where the superscript on Ai indicates that we have neglected initial state interactions. Next, 



we show that Eq. (6.20) is only an approximation to the full impulse amplitude derived 



from the relativistic Bethe-Salpeter formalism. We will see that this approximation is only 
valid in the absence of initial state interactions. 

For the case of plane waves in both the initial and final states, a one-body operator is 
forbidden by energy-momentum conservation. 



{P3,P4\On\pi,P2) = 0, 



(6.21) 



with all the pi on mass shell. The correct formalism must be able to explain the required 



energy transfer. Our primary thesis is that the diagram of Fig. 6.1 must be obtained from 
the Feynman rules as 



{r\GiO^\pi,P2) = \K{mj2) GiO^l pi,p2) , 



(6.22) 



where Gi is the Feynman propagator of the intermediate off-shell nucleon and K(m.,^/2) is 
the sum of all irreducible diagrams with energy transfer of 171^^/2. The second equality of 



Eq. (6.22) results from the relation between F and ^ in Eq. (6.7). This manipulation is 



necessary because {(p\ will be used for evaluation instead of (^'|, meaning that the relative 
energy must remain zero in the final state. Thus the full kernel for pion production via 
the impulse approximation is KGiOj^ rather than just O^. Because KGiO-,^ is a two-body 
operator, the momentum mismatch which suppresses the lA in the traditional treatment is 
removed. 

There are two points to emphasize here. Firstly, this treatment is not equivalent to the 
heavy meson exchange operators of Refs. \103\ [56] which are intended to account for the 
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relativistic initial and final state interactions not present in phenomenological potentials. 



Secondly, although the assertion of Eq. (6.22) greatly changes the way impulse pion produc- 



tion is calculated, one should not perform the same manipulations for the similar impulse 
approximation to photo-disintegration. The reason for this is simply that near threshold 
the nucleon remains essentially on-shell and the diagram is therefore clearly reducible. 

Next, we use ^ = ^ p + = (f) + and focus on the <f> term; the other term contains 
non-nucleonic physics and may be treated as a correction. Thus the impulse approximation 
is given by 



\Kimj2)GiOMP2). 



(6.23) 



Consider the spacetime structure of the product, GiOt^. The relativistic propagator Gi 
is decomposed into three terms: 1,7°, and 7*. Between two-component nucleon spinors 



7Si7° 



7%S./ 

y7Si7° 



and so we can make the replacement 



TUN 



(Pi - qY - + ie 



0, 



. -E(pi) -m^+ niN ^ 
-2E{pi)m.„ + ml + ie 



(6.24) 



(6.25) 



1 



where in the second line we have used that -E'(pi) = rriN + mTr/2 at threshold. Note that 
this propagator agrees with that obtained from the Feynman rules for B^PT at LO. 



In order to make connection with the traditional Eq. (6.19), we use the approximations 
K ^ U [corrections are 0{g — G)] and Gi ^ —i/m-,^ [corrections are 0(m^/mAr)]. Putting 



these substitutions into Eq. (6.23), 



m-n 



-a 



\PlP2)- 



(6.26) 



The quantity U is related to the potential energy hy U = —iV. Ignoring the fact that U 
should be evaluated for non-zero energy transfer, we use the equal-time Schrodinger equation 
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K 



Figure 6.4: Bethe-Salpeter formalism applied to pion production for plane wave initial 
states. 



to replace V — )• —E^ — jm^ and then neglect the binding energy to find ~ -^ib^- 

This means that for a PW initial state, the traditional impulse approximation should be 
roughly adequate. This is borne out in the actual calculation of the reduced matrix elements 



for Eqs. (6.20) and (6.26), 



4PW 
^IB 

^2B 



-24.0 
-25.6, 



(6.27) 



where we have used Appendix B.l 's definition of the reduced matrix element (we suppress 



the subscript on Appendix B.l 's for clarity) and used the same static phenomenological 
potential for V (here, Argonne vl8 |84| ) that is used to calculate the wave functions. See 
Fig. |6.4| for a pictorial description of this section. 

It is important to note that the Bethe-Salpeter equation can also be used for the RS 



diagram as shown in Fig. 6.5 In fact, the diagram on the right in Fig. 6.5 (with K 
approximated by OPE) is one of the two topologically reducible diagrams used by Ref. 
[5j to resolve the problem arising from calculation of NLO loops as we have discussed in 
previous chapters. 



In the next section, we will show that for distorted wave initial states, Eq. (6.26) is 
replaced by 



-'*^'2B 



(6.28) 
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Figure 6.5: Use of the Bethe-Salpeter equation in the rescattering amplitude 



where the first term contributes at leading order in the theory and the second term at 
next-to- leading order. 

6.5 The NN — > dvr Reaction: Distorted Wave Initial States 

6. 5. 1 Definition of distorted wave operator 

There is no reason to expect the result A^2B^ ~ -^iB^ ^'^ carry over for a distorted wave 
(DW) initial state where = m-n-mj\[ no longer holds. Indeed, we will show that the 
traditional expression for the impulse approximation does not hold for DW amplitudes. 
The fully- relativistic initial-state wave function is denoted \^)i, 

mi = \pi,P2) + GK\^)i, (6.29) 



where the first term is exactly the initial state used in the definition of A4^^ of Eqs. (|6.20|) 



and (6.22). The complete DW impulse operator is defined as, 



The second term includes the production operator KGiOj^ from Eq. (6.22) along with 
initial state interactions, 

M'M = f{^\KGiO^GK\^)i 

(6.31) 

^ f{(P\KGiO^GK\<j))i, 

where in the second line we have once again used = (/) + ^'g and neglected the Q-space. 
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(/"-m^ + ^,l-k) 




-i'-k/) 
-i,-k) 

I (-/^* + m^-7,-l) 



fm^r _ 7 _ ^ 



Figure 6.6: Impulse approximation using distorted waves. Solid lines represent nucleons, 
dashed lines represent pions, and ovals represent wave functions. 



As noted by Ref. [104j, the kernel of Eq. (6.31 ) is a loop integral which is shown in Fig. 



6.6 with K being approximated by OPE. Note that four-momenta are conserved at every ver- 



tex. One pion exchange is the first contribution to K in xPT besides a short range operator 



which is irrelevant for the s-wave NN — t- dvr amplitude (see Appendix A.2). Nevertheless, 
one must exercise caution due to the large expansion parameter of pion production. To this 
end, we employ the deuteron of Ref. [105] which is calculated from a purely-OPE potential 
with suitable form factors. As discussed in Appendix |D.2[ this deuteron wave function is 
quite accurate and increases the rescattering amplitude by only 3% over a phenomenological 
deuteron. Having then employed this deuteron wave function in the calculation of the tra- 
ditional DW impulse approximation, we will be able to avoid any complications from higher 
order parts of the potential in our subsequent investigation of the two-body operator of Eq. 



(6.31). In other words, although the full potential must be present in an exact calculation, 
we expect to gain insight into the correct formalism by using an OPE-only deuteron. We 
continue to use the phenomenological potentials for the initial state. To verify that the use 
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of K = OPE in the initial state does not spoil our results too much, Appendix D.5 examines 



heavy meson exchange in the initial state. As will be discussed, this effect is parametrically 
suppressed. 

Note that the relative momenta of the nucleons before and after the loop (p and k) are 
external momenta to the loop integral over / = (Z*^ , 1) , but are eventually integrated over in 
the hybrid formalism. Let us focus solely on the energy part of the loop integral and ignore 
the vertex factors and overall constants. We define the integral /, 

df 1 1 1 1 



Inl'^ -E + ie P 
1 



itItj — E + ie —P + m-7, 
1 



'J - E + ieP - - 7/2 + uji - ie 

1 



P — — 7/2 — LOi + ie P — + 7/2 + ujf — ieP — rriT^ + 7/2 — + ie' 



(6.32) 



where = (1 — p)^ + is the on-shell energy of the initial-state pion, = (1 — k)^ + 
is the on-shell energy of the final state pion, and E = l^/2mjv is the kinetic energy of a 
single intermediate nucleon. Note that pf mT^niN — ^niN and kj —7m at. 

It is straightforward to show that if the energy components of the exchanged pions in 
the above loop are set to zero (violating conservation of four- momentum) , one obtains the 
traditional impulse approximation. In this case, the pion energy denominators are pulled 
out of the integral which is then evaluated by closing the contour in the lower half plane, 

1 1 

-Wj —7 — \^ /rriN 

The quantity in the first set of brackets can be recognized as the product of OPE with 
the final state wave function while the second set is the product of the initial state wave 
function with OPE. This is precisely the operator that the traditional evaluation includes. 



/iB 



1 1 

7 - P/rTT-AT -Ojf 



(6.33) 



6.5.2 Reduction to time ordered perturbation theory 



Our goal is to evaluate the integral in Eq. (6.32), showing that it is a sum of diagrams 



from time ordered perturbation theory (TOPT) terms which can be combined to obtain 



Eq. (6.28). To begin, we rewrite the first two factors as a sum, 

1 



1 



1 



P-E + ieP 



E + ie 



-m-n 



1 



P-E + ie P 



E + ie 



(6.34) 



101 



>( 



X 



/"r term 



/9R term 



Figure 6.7: Definition of the two terms in Eq. (6.35). Crosses represent the propagators 



which are absent due to the partial fractions decomposition. 



This is the key to our method because after making this spht, we see two terms which 
each have the propagator structure of a rescattering box loop. Consider the first term 

- and 



in Eq. (6.34); this loop integral looks like a two-body operator multiplied by 



augmented with an initial-state interaction. The second term looks like the same with 
final-state interaction. We define these two integrals to be IgB ^^'^ -^2B' respectively, 



hB — I2B + Ab- 



(6.35) 



Figure 6.7 illustrates the splitting described in Eq. (6.35). 



Next, we perform partial fraction decomposition on each of the pion propagators, split- 
ting each of the two terms into four terms. Then, we continue the decomposition process 
until each term can be expressed as a single residue. For /ge ^^^^ isolate the poles con- 
taining iOf and then close the contour around them (for IgB' the iOi poles are isolated). By 
isolating the poles in this way, the resulting expression is easily recognized as the sum of 



six TOPT terms. For clarity, we show these terms pictorially for /ge ^i§- 6-8 where we 
have left the overall 3^ implicit. We assume for now that the stretched box diagrams are 
small, as they were in the rescattering toy model investigation [95] and denote the sum of 
the four remaining terms with a /. 

Finally, motivated by the interpretation which is presented in the next section, we alge- 
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Figure 6.8: TOPT terms resulting from the integral. 



braically re-combine these four terms to find 



-ra 



2^ (mj2)2-(a;^ + 5,(l))2 

1 4(1) 



1 + 



2B 



1 



1 



Sail) 



l^i + Sail) 



1 



1 + 



•j -2E - ^ [ ujf + 6b{V)i iUt, I uji 
where we have separated out terms involving 6a and (5^ 

12 



- 2^; - 7 -ijf 

(6.36) 

(6.37) 



(m./2)2-(u;, + <5,(l))2' 



2mN 

12 



7 

2 2 



2mjv 2 ' 



(6.38) 



because (as will be shown in the next section) they are sub-leading and we will neglect them 
in the main body of this work. The only approximation made in the evaluation of the loop 



integral to obtain Eqs. (6.36) and (6.37) is to neglect the stretched boxes. Let us pause to 



summarize what we have done so far: (1) the DW amplitude was written down as a loop 
integral, (2) partial fractions was used to split the product of the two nucleon propagators 
into a sum +/2B' (^) loop integrals were evaluated and the result expressed in terms 
of TOPT diagrams, and (4) the TOPT diagrams were algebraically combined into a form 
useful for the following interpretation. 

6. 5. 3 Interpretation 



Although not obvious at first sight, convolution of the operator corresponding to Eq. (6.36) 



with wave functions as defined in Eqs. (6.30) and (6.31) results in an amplitude approx- 



imately equivalent to that which one obtains by using the operator shown in Fig. 6.9(a 
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The same is true of Eq. (6.37) with Fig. 6.9(b), and together they replace the traditional 



(one-body) impulse approximation with Eq. (6.28). Furthermore, the operator that results 



from Eq. (6.37) is expected to be small by power counting arguments. The task of this 
subsection is to verify these statements in detail. 



In Eq. (6.36) the factor (myr — 2E — 7) ^{—ojf) is interpreted as the product of the 



two-nucleon initial-state wave function with static OPE. This is the statement that 

1 1 



2^-7 -ujf 



OPE 



(6.39) 



This factor can be absorbed (after adding in the PW term) using the zero-relative-energy 
Lippmann-Schwinger equation that is employed by the phenomenological potentials we are 
using. We will continue to refer to the initial wave function as a function of p and pj, so 
absorbing this factor means that we set 1 = p. 

Likewise, in Eq. (6.37) the factor (— a;j)~^(— 2£' — 7)"^ is interpreted as the product of 
static OPE with the two-nucleon final-state wave function: U'^^^g. Absorbing this factor 
into the wave function, we set 1 = k. The remaining factors of IgB ^'^^ -^2B become the 
two-body impulse production operators. 



2B 



(Tl • (p - k)(T2 • (k - p) 



1 + 



2B 



Sk 



1 



4(k) 



w/ + 4(k) 



1 



1 + 



1 



1 



5a(p) 



+ 5a{v) 



Sp 



(Tl • (p - k)<T2 • (k - p) 

(mj2)2-(u;, + 5,(k))2' 



(6.40) 
(6.41) 



where we have now made explicit the momentum dependences of the vertices and used 
S = (cTi -|- <T2)/2. It is also important to include form factors in the OPE which match 
those of the wave functions. These form factors are present in our calculation even though 
we leave them out of this expression for the sake of generality. 



Next, note that in the evaluation of the matrix element using Eq. (6.40), the initial state 
wave function is peaked about its plane wave value p ~ Pi, and thus E ~ m^/2 — 7/2 and 
5a(p) ~ 0. On the other hand, in Eq. (6.41), we have k ~ kj and E ~ —7/2 and (^^(k) « 0. 



If we were to neglect all the 5's, we would have 

<Ti • (p - k)<T2 • (k - p) 



1 



2B 



(mj2)2-((p-k)2 + m2) 



-Sp 



o 



2B 



S -k 



1 <Ti • (p - k)(T2 • (k - p) 

(m^/2)2 - ((p - k)2 + ml) 



(6.42) 
(6.43) 
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(a) O^B of Eq. ([6l2l 



(b) O2B of Eq. (6.431 



Figure 6.9: Two-body impulse production operators (pion exchange is non-static). 



which suggests that these operators can be approximately interpreted as the diagrams in 



Fig. 6.9 Thus we have finally obtained our central result, [Eq. (6.28)], which states 
that the correct impulse approximation is a two-body operator. The contribution to pion 



production given in Eq. (6.28) is not replacing the rescattering diagram (which is also 
two-body), but rather replacing the traditional contribution which has been referred to as 
the impulse approximation (or direct production). Note that if we assign standard pion 
production power counting to these diagrams. Fig 



6.9 a 



is O 



rriM J 



while Fig. 



O ( ) . In the next section the approximate expressions given in Eqs. (|6.42|) and 



6.9(b) 



IS 



6.43) 



are numerically evaluated. Nevertheless, we acknowledge the importance of verifying that 



the 6 terms are indeed small and relegate that discussion to Appendices D.3 and D.4 



6.5.4 Evaluation of two-body operators 



Next, we calculate the threshold s-wave np — )• dn amplitudes corresponding to Eqs. (6.42) 



and (6.43). We do not present the details here as most are given in Appendix C.2 Again 



we remind the reader that for the sake of consistency we use a deuteron wave function 



calculated from a purely-OPE potential (with form factors as described in Appendix D.2). 



For the initial-state distorted waves, we use three different phenomenological potentials 
(Avl8 m], Nijmegen II [96], and Reid '93 [96]). In Table |6]l| we display the results in 
terms of the reduced matrix elements of Appendix [B.l 
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Table 6.1: Threshold reduced matrix elements calculated with an OPE deuteron and various 
phenomenological initial states. The first row shows the traditional impulse approximation 
(one-body) while the second and third show our replacement (two-body). 



Avl8 Reid '93 Nijm II 



/iDW 
^IB 


8.3 


7.1 


5.4 


xDW.a 
^2B 


17.4 


13.5 


7.8 


.DW,b 
^2B 


-1.5 


-2.2 


-6.! 



The first row of Table 6.1 gives the traditional (one-body) impulse approximation, which 
is slightly bigger than in Chapter [5] due to the use of the OPE deuteron. The next row shows 
that the new two-body operator (at leading order) is roughly twice as large as the tradi- 
tional calculation it is replacing. We mention here that the significant cancellation between 
deuteron s- and d-states remains, keeping the impulse amplitude smaller than rescattering; 
however, the cancellation is less complete when using our new two-body operator. The final 
row verifies that the diagram is smaller than the diagram, as dictated by the power 
counting. The Nijmegen II potential provides a bit of deviation from these results, and it 
will be interesting to investigate other potentials to determine the true model dependence 
of this calculation. In finding these results, it is important that the pion propagators of 



Eqs. (6.42) and (6.43) be implemented in a manner consistent with the potential used for 
the wave function of Fig. |6.6[ Namely, the cutoff procedure of the convolution integral with 



form factors needs to match that by which the potential was constructed. Appendix D.2 
contains the details of this procedure. 

Our conclusion is that the traditional impulse approximation is an underestimate. While 
it is true that several approximations were made in order to permit final expressions as 



simple as Eqs. (6.42) and (6.43), we believe this conclusion to be sound. The S terms do 



not defy their classification as sub-leading (see Appendices D.3 and D.4), and Appendix D.5 



shows that using K = OPE in the initial state is at least reasonable. In summary, we 



simply claim that Eq. (6.42) is the new impulse approximation at leading order in the 
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effective field theory. The corrections in the aforementioned appendices, in addition to Eq. 



(6.43) contribute to the next-to- leading order calculation, which needs to be systematically 
considered in a later work. 

Finally, it is important to note that although the OPE deuteron reproduces the phe- 
nomenological results for the rescattering diagram quite well, the numbers in this section 
are greatly changed if a phenomenological deuteron is used. Using Avl8 we find = 4.9, 
and by using the cutoff procedure of Avl8 for the two-body operators, we find ^4^^'^ = 33.5, 
^2B^''^ = —2.8. Thus, the ratio of our new two-body operator to the traditional impulse 
operator is ~ 7 instead of the ~ 2 presented above. At this time one is faced with a choice 
of either: (1) using a "correct" phenomenological deuteron and leaving out parts of the 
potential when calculating the two-body kernel or (2) using an inexact OPE deuteron with 
a completely self-consistent kernel. For the time being, we believe the latter to be more 
trustworthy, if not ideal. 

6. 6 Discussion 



Recall from Table 3.1 the results for the total cross section obtained by the five most recent 



experiments in terms of the parameters a and /3 of Eq. (3.22). Since the present calculation 



is performed at threshold (r/ = 0), we compute only the value of a, 



a = \A\^ . (6.44) 

1287r2sp ' ' ^ ' 



For ease of comparison, we invert Eq. (6.44), plug in the results of the mentioned experi- 



ments, and propagate the errors to find Table 6.2 



The full theoretical amplitude includes not only the impulse diagram but also the rescat- 



tering diagram, which is given in Table 6.3 along with the total amplitude using either the 
traditional one-body or the leading-order two-body impulse diagram. The uncertainty in an 
effective field theory calculation is estimated by the power counting scheme. In this work, we 



have included both the rescattering and the impulse diagrams up to O y^/niT^/mNj . There- 
fore one might assign an uncertainty of m-^/mN = 14% to the calculation but stress that 
such an estimate based solely on power counting is rough at best. Taking this uncertainty, 
we see that the theory update presented here changes the situation from under-prediction 
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Table 6.2: Threshold reduced matrix elements extracted from experiment 

Experiment ^expt 
np dn^ (TRIUMF, 1989) [65] 
TT+d pp (1991) 166J 
pp dTT+ (lUCF, 1996) [67] 
pp dTT+ (COSY, 1998) [51 
Pionic deuterium (PSI, 2010) [48] 



80.1 ± 1.1 
77.9 ±0.7 

85.2 ± 1.0 
84.6 ± 1.9 

93.8^°:^ 



Table 6.3: Rescattering and total reduced matrix elements for a variety of potentials. The 
second line shows the traditional calculation (with a one-body lA) while the third shows 
our replacement (with a two-body lA). 





Avl8 


Reid '93 


Nijm II 


RS 


69.8 


72.1 


74.0 


RS + lA (IB) 


78.1 


79.2 


79.4 


RS + lA (2B) 


87.2 


85.6 


81.8 
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of the most recent pionic deuterium experiment by ~ 1.3 cr, to under-prediction by ~ 0.7 a. 

In summary, we have developed a consistent formalism that allows one to separate effects 
of the kernel from those of the wave functions, finding a new impulse approximation kernel. 



This two-body operator, given in Eq. (6.28), replaces the traditional one-body impulse 
approximation and is the central result of the present work. We numerically investigated 
the simplest example (s-wave A^A^ — )• dir) and found the impulse amplitude to be increased 
by a factor of roughly two over the traditional amplitude. This calculation was performed 
with a regulated OPE deuteron which has advantages and disadvantages as described in the 
body of this work. Rescattering remains the dominant contribution to the cross section. We 
find that the updated total cross section is ~ 10% larger than before and is in agreement with 
experiment at leading order. We verified that corrections to the new impulse approximation 
(which together with other loops and counterterms will contribute at next-to-leading order) 
do not destroy these results. 

These findings suggest several directions for future research. Firstly, one needs to de- 
velop a power counting scheme for the "Q space" discussed in Section |6.2[ Secondly, the 
significant model dependence of the new formulation of the impulse approximation needs to 
be investigated in a renormalization group invariant way. Thirdly, it will be very interesting 
to see the impact of this increased impulse amplitude on the pp — )• ppir^ cross section which 
is suppressed due to the absence of rescattering. Finally, one could look at the energy de- 
pendence {p-wave pions) of A^A^ — t- NNtt, for which there is an abundance of experimental 
data. 
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Chapter 7 
DISCUSSION 

In this thesis, we have reported on two studies: Chapter |4] described the calculation of 
the asymmetry of np — )• di:^ resulting from CSB (most importantly the 5mi^ interaction), 
Chapters [5] and [6] eventually led to the discovery that the hybrid formalism forbids the 
traditional one-body impulse approximation, replacing it with a two-body operator. Let us 
now summarize these results and discuss the outlook for these topics. 

We performed the asymmetry calculation up to partial-NLO and found = 40.2 x 
10"^ which is to be compared with A^^^ = [17.2 it 9.7] x 10"^. To obtain our result for 
the asymmetry, we used the Cottingham sum rule to fix (5mjv; the other CSB coupling, /3i, 
was fixed by assuming it to originate from an intermediate rj which fluctuates into a vr''. 
Both of these couplings are not precisely known, and it is possible that revised values could 
improve our calculation's agreement with experiment. Nevertheless, the disagreement in 
we see when using the best available CSB parameters is significant. It is important to 
understand to what extent this disagreement is due to the CSB couplings themselves, and 
to what extent it is due to incompletely- understood theoretical issues. 

In Chapter [4| an unexpectedly large amplitude in the ^5*0 channel led to a disagreement 
with p-wave data and was fixed with an ad- hoc cutoff. The fact that Ref. [91j, with its A^A 
coupled-channel wave functions, did not find a similarly large amplitude suggests that, with 
our A^A^ wave functions, we need to perform the p-wave calculation to higher order in the 
EFT. As we discussed, the first set of loops vanish. The NNLO calculation includes various 
tree-level diagrams (involving the known q LECs) along with a vrA^^ contact interaction 
(involving the unknown d LEC). Recently, Ref. |Mj performed this NNLO calculation and 
fit d with the relevant experimental data. The determination of d requires an unambiguous 
knowledge of all other p-wave diagrams up to the order where it appears. In this thesis we 
have uncovered a problem with the way the lA diagram has been calculated. Therefore, a 
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new determination of d should now be performed. 

A second difficulty encountered in the CSB study regards the total cross section of 
np — )• dir^. In this thesis, we applied recent power counting and reducibility developments 
to calculate a total cross section roughly in agreement with collider data. At the same 
time, pionic deuterium experiments have been shown to imply a cross section well above 
this prediction. The study of Ref. [91] used the cross section obtained from the atomic 
experiment as an ingredient in their asymmetry calculation. This fact accounts for most 
of the discrepancy between our two calculations. Why do the atomic experiments disagree 
with the collider ones by such a large amount? 

As part of the CSB study of Chapter |4j we became interested in the energy transfer 
mechanism of the lA diagram. A naive attempt at a solution (wave function corrections) 
was made in that chapter. This solution was investigated in greater detail in Chapter [5] 
and the results were not satisfactory. Significant model dependence was observed resulting 
from the use of different phenomenological potentials. With a correctly implemented hybrid 
approach (including correct cutoffs/form factors), this should not be the case. Additionally, 
a concerning departure from the MCS was seen when the intermediate nucleon propagators 
were treated properly. Rather than attempt to fix these issues from within the wave func- 
tion correction paradigm, we took a step back and thought more deeply about the hybrid 
formalism. 

In Chapter [6j we investigated the algebraic reduction of the full 4D loop-integral which 
occurs in a distorted-wave calculation. We were able to clearly (even numerically) identify 
the traditional lA as an approximation that is only valid in the absence of initial-state 
interactions. By making well-controlled approximations, we were then able to identify the 
true distorted-wave lA: a two-body operator which replaces the old one-body operator. 
We are confident in our identification of a new LO expression, but the reduction from 
the relativistic formalism to the standard hybrid formalism involved a large number of 
approximations which remain to be classified in the EFT expansion. There is a significant 
task ahead for a complete calculation. 

Despite these difficulties in the lA study, the result remains that the s-wave lA amplitude 
is enhanced by a factor of approximately two. This promises a significant numerical change 
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to the calculation of p-wave pion production (recall that for p-wave pions, the lA is not 
washed out by a large RS amplitude). A recalculation of both the / = 2 Legendre coefficient 



and the analyzing power, discussed in Section 4.3.3, is the natural place to begin. We will 
then be prepared for another recalculation of the asymmetry. 

A different avenue to pursue when applying the new formalism is the theory of the 



pp — )• ppn^ reaction, which, as discussed in Section 3.5 is quite sensitive to the lA. Here one 
needs to re-evaluate the two-body s-wave amplitude of Chapter [g] between the initial ^Pq 
and the final ^Sq pp scattering wave functions, while correctly accounting for the three-body 
phase space. Because it has been known for a long time that the traditional one-body lA 
is too small to account for the total cross section by itself, it will be interesting to see how 
much bigger the new two-body operator is. 
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Appendix A 
GENERAL APPENDICES 

A.l Fundamental Constants and Parameters 



Some of the constants used in this thesis are given in Eq. (A.l). 



m^o = 134.98 MeV qa = 1-32 
rriTv = 938.92 MeV = 92.4 MeV 

ruA = 1232 MeV Ha = 2.1gA 

Eb = 2.224 MeV he = 197.327 MeV • fm 



(A.l) 



Note that different values for gA and f-,^ were used in Chapter [4| The LO values in Eq. ( A.l ) 
are more apt for Chapters [5] and [6] because they are LO calculations. 

A. 2 Lagrange Densities 

This appendix summarizes the Lagrange densities developed and used in this thesis. Most 
of these expressions were developed in Ch. [2] in terms of the pion matrix U. For the present 
purposes it is sufficient to expand up to two-pion interactions. Terms with more than two 
pions, as well as other terms not used in this thesis are omitted, and for this reason we 
write a ... at the end of each expression. In this section the superscript in , rather than 
indicating the number of derivatives or powers of y/mq, indicates the index of a Lagrange 
density, 

u = d+^-2, (A.2) 

where d is the sum of the number of derivatives and powers of m-,^, and / is the number of 
fermion fields. This represents the standard power counting for nuclear physics. 



The starting point is the general Lagrange densities in Eqs. (2.34), (2.56), (2.57) 
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and (2.63) to (2.65), to which we apply the expansions 

1 



U = 1 + —TaTTa - ^JTa ((^aft + i(^abcTc)T^bT^c + ■■■ 



1 



^abcTaT^bduT^c + • • • 



(A.3) 



Ta<9^vra + . . . . 



We also rewrite the quark mass matrix M = (mu+md) /2+{mu—md) /2t3, use B^niu+Trid) 
m^, and employ v = (1,0) in the spin-3/2 projection operators. 



3/2 v={l,0) _ 3/2 _ 3/2 
3/2 v={l,0) 1 



(A.4) 



Furthermore, employing the heavy-baryon formalism, one is free to replace 7j7j — !■ —aiCjj. 
Finally, it is common to express the Lagrange density in terms of an un-projected field 
and leave the (iso)spin transition matrices {TabjSij explicit. 



Tab = Sab - T^TaU 



Si 



1 



(A.5) 



which satisfy (TT"^)af, = Tab and {SS^)ij = Sij. 

The u = Lagrangian of HBxPT (with isovectors in bold font) with the Delta included 
as an explicit degree of freedom is 



N 



1 

4^ 



(r • (tt X tV)) N + ( r • <t • Vtt) iV 

hA r 



2// 



(A.6) 



+ At {ido -6)A+-f- N^ViTTaSijTabAl + H.c 



+ 



where r and a are the Pauli matrices acting on the isospin and spin of a single nucleon. 
The i/ = 1 Lagrangian includes propagator corrections, recoil terms, and the leading 
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s-wave CSB operator 
1 



1 



2mjv 2mN 



1 



2mN 
6m]\f 



4/, 



1 

1 2hA 



2mN 2f^ 
2 

2vr3r • TT 



TT X Vtt • ViV 



9A 



2U 



\N^T ■ -ka ■ VN + H.c. 



(A.7) 



As previously mentioned, several terms have been omitted from this expression. Two pion 
seagull terms with both derivatives on the pion fields do not get promoted in the MCS 
scheme. This same property applies to terms with the Cj low energy constants that appear 
at this order. Terms with the di low energy constants do not contribute to s-wave pion 
production. Finally, the NNNN contact terms Cs,t do not contribute because we are 
using a potential with a repulsive core [Ri{r)Rf{r) — t- as r — t- for li = 1, If = 0]. 

Although there are a host of = 2 terms, we just list the CSB term relevant for this 
calculation 

Pi 



£(2) = i^Arta^.W3Ar + .... 
A. 3 Calculation of Nucleon-Nucleon Wave Functions 



(A. 



A. 3.1 Problem definition 

In this section, we describe the numerical exercise of calculating scattering and bound state 
wave functions from phenomenological potentials. The appendices of Ref. |1061 [33] very 
helpful to this end. The starting point is the time-independent Schrodinger equation for a 
state with spin S, orbital angular momentum L, total angular momentum J, and magnetic 
quantum number mj, 



niN 



-V2 + y(cri,o-2,r; 



ULj{r) 



25+1 



pr 



Lj,mj{r)) 



m]\j pr 



2^+iL,,m,(r)), 

(A.9) 



where is the twice the reduced mass for the nucleon-nucleon system, r = (ri — is 
the relative coordinate, and p^/m^v is the kinetic energy in terms of the relative momentum 
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p = (pi — P2)/2. The spin-angle wave function is defined according to, 

I ^^'^^Lj,mj(r)) = ^ I {Sms,Lmj - ms \ Jmj)YLmj-msi^) I Snig), (A.IO) 

ms 

where the (•• | •) indicates a Clebsch-Gordan coefficient. Phenomenological potentials are 
available in the form of data tables for each individual partial wave, i.e. V{r) for r = 
0.01, 0.02, 15.0 fm, 

VsLj{r) = j dr{^^+^Lj, mj{r) \ V{<ti, <T2, r) | ''^+^Lj, mj(r)) (A.ll) 

The deuteron has the added complication of being a — ^Di coupled channel np state; 
the unbound version of this state is also relevant for CSB matrix elements. In this case the 
makers of the potential typically provide the off-diagonal term, 

K(r) = [ dr{^Di,m{r) \ V{(Ti,a2,r) \ ^Si,m{r)) 

, (A.12) 
drfSi,m{i:) \ V{ai,cT2,r) \ =^1^1, m(i-)). 



A. 3. 2 Solution for uncoupled channels 

First we will describe the method for finding the wave functions for the channels in which 
the potential is diagonal. In this case, the Schrodinger equation reduces to. 



1 /i#,_ia±i)Vno 



y "^=0. (A.13) 

rriN \r or'^ / tun \ pr 

Since we will be integrating this equation out numerically, what we need are the solutions 
to this equation in both the r rmin limit where the centrifugal barrier dominates, 

fd^ l{l + l)\n{r) ^ u^^^^l^B_^ 

and the r — >■ rmax limit where the potential vanishes, 

where ji{ni) is the spherical Bessel (Neumann) function. Note that we can discard the 
irregular solutions (the B and D terms). Also note that the situation is different for Z = 

which has no centrifugal barrier; here, we find the solution near the origin to be oc sinh(^r) 
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with = mj\fV{r^i^) —p^. The effect of the potential in distorting these solutions (at least 
for elastic scattering) is to introduce a phase shift at large r, 

Ul{r) r^po 



sin(pr - ^) ui{r) r^po ^sin(pr - !f + 5i) 



C- 



pr pr 
A few lines of algebraic manipulation then yields, 



pr 



ui[r) 
pr 



C{cos{6i)ji{pr) - sm{5i)ni{pr)) . 



(A.16) 



(A.17) 



At this point, we are able to perform the numerical integration via the following algo- 
rithm: 



For / 7^ 0: Taking ui{rrain) = u* and u'i{r 



mm; = as the boundary conditions, 



integrate the full differential equation of Eq. (A. 13) from rmin 

to r 



For / = 0: Taking UQ^r^ 



u* and nQ(rjnin) = u*S, coth(,^rmin) as the boundary 



conditions, integrate the full differential equation of Eq. (A. 13) from rmin to r 



Solve for the overall normalization and phase shift by matching onto the asymptotic 



form given by Eq. (A. 15) and its derivative. 



A. 3. 3 Solution for coupled channel 
For the coupled channels we find. 



1 1 ^2 ^ ^ p2 

TTT^ + ^o r - ^ 



u{r) _ _y^^^-^w{r) 



1 92 



6 



ruM \r or^ 



P 

ruN 



pr 
w{r) 
pr 



-Veir) 



pr 
u{r) 
pr 



(A.18) 
(A.19) 



where u{w) is the / = 0(2) solution. Equations (A.18) and (A.19) admit two solutions 



{ni,tt;i}, {^2,^2}- One may obtain independent solutions by solving the equations once 
with the boundary condition u{r^i^) = and then again with w{r^\a) = 0. These two 
solutions are combined in two different linear combinations to find the two true physical 
solutions {mq,, tOa}, {u^^wp}^ 

Ua = Aui + BU2 Uj3 = Cui + Du2 

Wa = Awi + BW2 Wp = Cwi + DW2- 



(A.20) 
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The free solution for the a channel can be expressed as a sum of incoming and outgoing 
spherical waves, 

^ o/ij^^ (pr) + bh'^^^ {pr) (A.21) 

^ c/i5^^(pr) + dhf\pr), (A.22) 

where, because the normalization will be fixed later, we can set 6=1. The incoming and 
outgoing solutions are related by the scattering matrix S according to. 



(A.23) 



In the "Stapp" convention |1U7| which yields the "nuclear bar" phase shifts, the scattering 
matrix is written, 

/ g2i<5„ cos(2e) ie^i^'^+^ii) sin(2e)\ 
S=\ , (A.24) 

\ie<^o.+Sfi) sin(2e) e^'^^P cos(2e) J 

where (5q,((5^) is the phase shift of the / = 0(2) solution and e is the mixing anglej^ 

The final step is to require that (by definition) in the a channel, the u solution must 
have an outgoing wave with amplitude e^**^" ; this means that we set a = e^^^" and solve for 
c and d, 

^ ^2i5^ cos(2e) + ie^(^"+^^) sin(2e)(i ^ d = -ie^^^--^p) tan e (A.25) 
g ^ jgi(<5a+<5/3) sin(2e) + e^^'^'' cos(2e)d c = ie'^^'^+^^ha.ne. (A.26) 

A few more lines of algebra yield the final asymptotic forms for the a channel (up to an 
overall normalization) , 

Uair) _^ sm{pr + 6a) 

P'' (A.27) 

Wair) cos(pr + TT + 6 b) 

> tane . 

pr pr 



^The older "eigen" phase shifts HQS] use a scattering matrix, 
5" = 



e^''^" cos^ e + e^'*** sin^ e (e^"" - e^'*" ) sin e cos e 
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To obtain the /3 channel solutions, we require that the outgoing wave part of have 
amphtude e^^^i^ . The algebra then proceeds analogous to the a channel we just described. 



Using Eqs. (A.20) and (A.27), their derivatives, and the corresponding /3 channel equa- 
tions, one is able to solve for the phase shifts normalizations necessary to specify the coupled- 
channel wave functions. We refrain from displaying the rest of the formulas here and mention 
one final channel, the deuteron. 

The deuteron is the bound state of the scattering channel we just described. Its wave 
functions' asymptotic forms do not involve any phase shifts and can be simply expressed, 

n(r) ^ Ae-^'' 

where 7 = \/m]sfEb is the binding momentum and r] is the asymptotic d to s ratio. 
A. 4 Spin-Angle Reduced Matrix Elements 
A.4-1 Wigner-Eckart Theorem 

When performing the calculations in this thesis, we have taken care to consistently apply 
partial wave expansions. As part of such a calculation, on encounters matrix elements of 
the form, 

{{S'L')J'mj I I {SL)Jmj), (A.29) 

where is a tensor operator of rank k and projection m which is created by coupling one 
of the spin vectors, S = (cri + <T2)/2 or S = (<ti — cr2)/2, to one of the spherical harmonics 

= {^rns,lmi \ kq)S^Ji^^{r), (A.30) 

nis ,mi 

where S represents either S or 5] and the angular momentum addition rules must be satisfied: 
k > |/ — 1|, k < / + 1, and + rrii = q. One must take care to use spherical tensors; a 
rank-one spherical tensor is related to its Cartesian components by 

^±1 = ± iVy), Vo = V,. (A.31) 
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This connection allows one to relate the position vector to the spherical harmonic Yi 




The Wigner-Eckart Theorem states that a general matrix element {j'm' \ 



(A.32) 



Ig I jm) can 

be expressed as the product of a "reduced matrix element" {j' \\ \\ j), which does not 
depend on the projections, and a Clebsch-Gordan coefficient, 

{jm, kq I j'm' 



{j'm' I I jm) = II || 



(A.33) 



In this way, one only needs to calculate the reduced matrix element once, simplifying cal- 
culations (such as the present one) that involve sums over the projections. The reduced 
matrix elements can be obtained by evaluating one particular matrix element explicitly. 
The reduced matrix elements we will be using are. 



(1 


1 s 1 


1) 


= V6 




5: 1 


0) 


= V3 


(0| 


s 1 


1) 


= -V3 


{l'\ 


Yl 


\l) 


Ll 
\/47r 



(A.34) 



{10, LO I I'O), 



where J = V2J + 1. 

When the tensor is formed by coupling two individual tensors which act on different 



spaces as in Eq. (A.30), one uses a general formula which is worked out in many textbooks 



{S'L'f II Tfc II SLJ) = J'kJ < 



see 


Ref. 


m 


's' 


L' 




S 


L 




S 

< 


I 





(A.35) 



where the 3x3 object enclosed in the braces is a 9j-symbol (a sum of products of Clebsch- 
Gordan coefficients). This expression becomes even simpler when the is a scalar produclj^ 



^In terms of spherical components a scalar product is given hy S ■ V — ^^{—l)"^SmV- 
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S-r, 



{S'L'f II 5 • f II SLJ) = 6j'j{-lf+^'+^J { 



S' V J' 
L S 1 



{S' II 5 II S) {L' II i II L), 



(A.36) 



where the 3x2 object enclosed in the braces is a 6j-symbol (another sum of products of 
Clebsch-Gordan coefficients). 



A.4-2 Evaluation of relevant reduced matrix elements 



To illustrate the methods described in Appendix A. 4.1 this section will tabulate the reduced 
matrix elements which are necessary to calculate the pion rescattering diagram amplitudes. 



This diagram is treated as an example calculation in Appendix B.2 First consider the 
operator S • f which contributes to s-wave pion production {^Pi — )• ^5*1, "^Di), 



((10)1 II S II (11)1) = V2 
((12)1 II S -til (11)1) = 1. 



(A.37) 



Next consider the operator S • ryi(r) which contributes to p-wave pion production 
(^5*0, — )• "^Di). One way to evaluate this matrix element is to insert a complete set 
of states to the right of the S • f . This set will collapse to one term, ^Pi, 

((IL')l II S -rYiir) \\ (OJ)J) = ^((IL')l || 5] • r || {SL)J) {{SL)J \\ Yi{r) \\ (OJ)J) 

SLJ 

= {{iL')i II s • i- II (01)1) ((01)1 II yi(f) II (oj)j). 

Now one simply needs to evaluate the individual reduced matrix element^ 
((10)1 II 5] • i- II (01)1) = 1 ((01)1 II n II (00)0) = 

((12)1 II 5] • f II (01)1) = -V2 ((01)1 II n II (02)2) = 



3 

Ait 

6 

47r" 



(A.38) 




(A.39) 



The uncoupled tensor operator Yi between coupled states can still be evaluated using Eq. (A. 35 1 by 
recognizing that the spin-space operator is simply the identity operator. 



131 



Appendix B 
CH. H DETAILS 

B.l DeGning Reduced Matrix Elements 

One can show that for s-wave pions, the production operator is always either a scalar 
or a rank-two tensor while for p-wave pions, it is always a rank-one tensor. This guides 
the following definition of the reduced matrix elements. Note that the Clebsch-Gordan 
coefficients which will be summed over as well as the Spherical Harmonics describing the 
angular distribution of the differential cross section are "pulled out." First we define the 
strong reduced matrix elements, 

r I . AStr I A f ^ A (1 "^f ' 2 I IrUf) ^ \ 

/ 1 MtU I .) = {^^A„ + ^ ^ ' ■^^-) ^^^j 

1* 



X (lm^,lm/ -rus \ lmf)Y^*_^^{p) 
f I -M?r=i I ^) = ^Be^''>cg,5„y,oY,\p) + ^^^^^^^^^^Ce^^^cgoy^;(p), (B.2) 

where rus = mi + 1712, egg = (1/2 mi, 1/2 m2 \ 0), and cg^ = (l/2mi, l/2m2 | 1 m^). 
Aq and A2 are the s-wave reduced matrix elements, and B and C are the p-wave reduced 
matrix elements. To clarify the notation consider A2, for example. 

A, = jirr^ (^^((10)1 II +^^((12)1 ||) Alf- (4.i^ || (11)1)) , (B.3) 

where the subscript J = 2 on the M indicates that we are using the portion of the operator 
that is a rank-two tensor in the space of total angular momentum. Note also that we have 
used the following general definition of a reduced matrix element, 

{{S'LV'm'j I T^^ I {SL)Jmj) (-^^jJp^^'j^ i^^s'LV II II {SL)J). (B.4) 
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Similarly for the CSB reduced matrix elements, 



(lm/,10 I Im/) / .^^ ^ 



X (lm,,2m/-m, | 1 m/) y^*_„^ (p) 

(B.6) 

where Aq and A2 are the s-wave reduced matrix elements and B, C, and D are the p-wave 
reduced matrix elements. Also note that the strong phase shifts have been denoted 6l for 
each of the three initial channels, the CSB phase shifts are denoted 6l for the ^Pi and 
^1)2 channels, and the coupled channel phase shifts are 6a and (5^. Since the ^Si and ^Di 
channels are coupled, B and C are split into a and /3 parts which have different phase shifts 
in the presence of initial state interactions. 

Finally, for comparison purposes we include a translation between our reduced matrix 
elements and those of Ref. |44j , 

Co = -^Aoe'^" 



Ci = -i^Be'^' (B.7) 
V Svr 



C2 = Hl—Ce'^\ 



3 

V 47r 

where the "C's" are pp — t- dn^ amplitudes and isospin symmetry has been used to determine 
the translations. 

B.2 Diagram Technique 



To establish the diagram technique consider Fig. B.l in light of Eq. (3.14). We write down 
the amplitude using momentum space Feynman rules, Fourier transform, and then convolve 
with the initial and final state wave functions. The left line is taken to be "nucleon 1" and 
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h = {k'l k - q/2) k2 = {kl -k - q/2) 



Figure B.l: Strong rescattering diagram. Solid lines represent nucleons, dashed lines repre- 
sent pions, and the pions' isospin z-components are 3 and a. 



we make the approximation that the exchanged pion carries half of the produced pion's 
energy, = ((x>g/2,q'). Momentum conservation gives q' = k — p + q/2. According to 



Eq. (A. 6) the WT vertex contributes 1/(4/^) ri^;, easb (a;g/2 + ooq) while the irNN vertex 



contributes qa/ {2fn) T2,a "'2 • (— q')- The momentum space propagator is — i/(q'^ + /i^) where 
/i^ = — LOg/4 is the effective mass of the rescattered pion. Next, as discussed in Sec. 
we Fourier transform with respect to 1 = k — p, 



3.4 



il r '^2 • q 



q'2 + ^2 



-jq-r/2 



e-*'i-'"/2cr2- (-iV) 



r "•2 q 

q'2 + /x2 



Airr 



^e-^'i-^/2^(r)^2-r, 
47r 



where h{r) = (1 + l//ir)e~^^'/r. 

The deuteron has isospin and the np wave fuction includes a T = 1, = isospinor 
|1, 0), and thus 



(0,0 I iea3bri,bT2,a I 1,0) = -2. 



(B.9) 



At this point, we have (defining M' = M / ^/2W[2E^2Ed) 

;0,0 I A^l'J 1,0) = -^m^M(r)e--r/^ 



0-2 ■ r. 



(B.IO) 



134 



To calculate the diagram with the WT vertex on nucleon 2 we consider how each part of 



the left side of Eq. (B.IO) transforms under 1 -H- 2. Since the strong part of the Lagrangian 



is invariant under isospin, Ai is invariant. The initial isospin ket, |1,0) is invariant as well, 
but |0, 0) — )• — |0,0). Also note that r — )• — r. Thus, 

. QA 3Wq/2 



0,0 I M'r. I 1,0 



/i h{r) e 



iq-r/2 



(Ti • r. 



(B.ll) 



(B.12) 



2fn 8^/2 

Defining S = (cri + a2)/2, = {ai - (T2)/2, and 

£ = exp(iq • r/2) + exp(— iq • r/2) 
O = exp(iq • r/2) — exp(— iq • r/2), 
we have the complete rescattering operator, 

0, I M'ns I 1, 0^ = -i 7i?5 h{r) (SS-f + O'E-r) 

where 

To proceed, we preform a partial wave expansion on £ and O and just keep the leading 
term. Note that we use the coordinate system defined hy q = z. Then we calculate the 
spin-angle matrix elements of the rank zero £ S ■ f and the rank one O'S ■ r operators (see 



(B.13) 



(B.14) 



Appendix A. 4. 2). The final expression for {f \ Ai \ i) simplifies to a radial integral which is 



computed numerically. Note that the first term in Eq. (B.13) corresponds to s-wave pions 
because £ carries L 



0, 



£ = 2x4njo(^)Yoo{r)Y,*o{ci) + ... 



2 jo 



qr 



(B.15) 



while r carries L = 1, 



= 2xA7Tin ( y) y^YimmUq) + • • • 



m 



(B.16) 
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and thus the operator will change the parity. Likewise, the second term corresponds to 



p-wave pions. In terms of the reduced matrix elements of Eqs. (B.l) and (B.2), we have 



A^^ = ^/2Ei 2E2 2Ed 8tt jrs V2 Ki 



A. 



RS 



B 



RS 



/2Ei 2E2 2Ed 8tt Kq 



C^^ = ^/2Ei2E2 2Ed8TT-fRsVQK2, 
where the integrals are defined 

„2 



Ki = 

K2 = 



dr r' 
dr 
dr 



Ud{r) , Wd{r 



r 


V2r ) 


Ud{r) 




r 


V2r , 


Ud{r) 


^Wd{ry 


r 


^/2r , 



Jo 



qr 



h{r] 



Ji 



h{r) 



h{r) 



ui,i{r) 
pr 

UQfl{r) 



pr 

U2,2{r) 

pr 



B.3 Observables 



(B.17) 



(B.18) 



One experimental observable is the analyzing power, Ay, defined in Eq. (4.19). In the 
strong sector, we find (neglecting A2 which is numerically small) 

V3cos{<p) sin(0)^o {V^B sm{6o - Si) + C sm{S2 - <5i)) 



(B.19) 



3Al + B^ + C^+ (C2 - 2V25Ccos(<52 - <5o)) P2(cos^) ' 
where the angular dependence is that of the nucleon relative momentum, p, with respect 
to the pion momentum, q = z. To compare with experimental results, which use p = z and 
(pT^ = 0, we need to set (pN = ^i" and On = 9-,^. 

To calculate the differential cross section as well as the asymmetry, we need to square 
the sum of all the matrix elements, sum over ruf and average over mi and m2. First we 
define 



Mo + AfiPi(cos6') + M2P2{cose) + MsP^icosi 



(B.20) 



so that 



A 



fb 



64 vr^ s |p 
Ml - \Ms 
2Mo 



AttM, 







(B.21) 
(B.22) 



136 



The results for the required quantities are 
1 



Mn 



Ml 



Mo 



M, 



48 vr 



(B.23) 



+ 2 + C^Cp) cos (5« - 5p)] 



\/2C ( 



Ac, 



1 



5V10 



^2 COs((5i - (52) 



+ Uo + 



10 



A, 



1 



An 



1 



5V10 



1 



SVlOvr 
5 



{A0A2 - 2A0A2) 



Ba + -^C!a ] cos{5a - Si) + (a ^ f3) 



A2 DC0S(J2 -(5l) 



^ {Al-2Al) 



(B.24) 



2V10 



+ 



3VW^ 2 " 2 ^ 2 



CaCp COs{5a - 5p)) 



V5 



(B.25) 



;25Ccos(52 - 5o) - BaCa - BpCp - {B^Cp + 5^C„) cos(5„ - 5p)) 



I (b^ + -\Ca ) ^ cos((52 - (5„) - (a ^ ^) 



'3 

40 7r V 5 



CA2 cos((5i 



1 

71 



A2D cos((52 — (5i) 



(B.26) 



Disregarding the small M3 term, the asymmetry is proportional to Eq. (B.24). The 



physical content of Eq. (B.24)'s first line is the interference of strong p-wave pions with 
CSB s-wave pions and the second and third lines are strong s-wave and CSB p-wave. 



Table B.l shows the strong reduced matrix elements and Table B.2 the CSB reduced 
matrix elements. The CSB rescattering numbers were calculated including the new contri- 
butions discovered by Ref. |91) . The np phase shifts (in radians) which appear in the cross 
section according to Eqs. ( B.23||bI26 ) are given in Tables B.3 and B.4 



Table B.l: Strong reduced matrix elements 



Diagram 


^0 


A2 


B 


C 


Impulse (w/ wfn corr) 








-6.59 


32.85 


Impulse Recoil 


5.62 


0.21 


1.17 


-8.56 


RS (w/ Recoil) 


81.12 





0.54 


1.66 


Delta (no cutoff) 








-33.94 


37.96 


Delta (A = 417 MeV) 








-10.48 


21.60 



Table B.2: CSB reduced matrix elements 



Diagram 


3o 


M 




Bp 






D 


Impulse (x^) 








12.23 


29.72 


-7.80 


-15.05 


-28.30 


RS(x^) 


-28.83 





-1.37 


1.79 


1.48 


1.89 


-4.92 


Delta (x-^) 








12.60 


-7.71 


-8.47 


-2.86 


22.37 



Table B.4: CSB phase shifts 
Table B.3: Strong phase shifts 



= 5i -0.44 

(5i -0.47 

5„ 0.19 

(5o -0.044 

5b -0.43 

(52 0.16 

= 52 0.44 
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Appendix C 
CH. E DETAILS 

C.l Impulse Approximation Details 

This Section begins with a presentation of the details of the calculation of the traditional 
lA pion production amplitude. Evaluating the isospin matrix element 

(00 I ri,3 I 10) = 1 (C.l) 



and using the vertex rule shown in Fig. 5.2, we obtain for Fig. |5.3Ka) at threshold 



00 I M'dp, k) I 10) = f^^^i • (P + k)(27r)^<5^(p - k), (C.2) 

where Ai' = Mj \J2m^ 2mp 2md = M/N. Since we are using position space np wave 
functions, we Fourier transform the matrix element with respect to 1 = k — p, which is 
identical to q' at threshold, 

d-^^ ^il-r/r,„\3e3 



^2^^3e^'-(2^rd^(p-k) = l (C.3) 

Note that we group the p and k with their respective wave functions prior to performing 
the Fourier transform 

ai-{p + k)^ai-{-i^np + i^d)- (C.4) 
Thus the full position space operator is 

00 I M'^ir) I lo) = -i^^ai • (^„p - ^,). (C.5) 
/ 2/^ 2mN 

To calculate the diagram with rescattering on the other nucleon, we consider how each part 



of the left side of Eq. ( C.5 ) transforms under 1^-7-2. Since the strong part of the Lagrangian 
is invariant under isospin, A4 is invariant. The initial isospin ket |1,0) is invariant as well, 
but |0,0) - |0,0) and V -V. Thus, 

00 I M'^^nir) I lo) = -2^|^^S • i^np - ^d). (C.6) 
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The final spin-angle wave function is that of the deuteron, while the initial state for 
s-wave pion production is solely ^Pi, 



/(r)) 



u{r) 3 



Si) + 



w{r) 



^(r))=4v^^"i^|3p0 
pr 



(C.7) 



where we have absorbed the unobservable (since there is only one initial channel avail- 
able) phase into the definition of the matrix element. The spin- angle matrix elements are 
calculated, 



(/(r) II S- (^-^) |h(r))=4^i 



, ,dRi{r) , ,2 , , dRf(r) , ' 



(C.8) 



where Rf{r) = ^/2u{r) jr -|- w{r') /r, Rf^2{f) = V^u{r) jr — 2w{r) /r, and Ri{r) = ui^i(r) /pr. 
Using Eqs. ( |C.6 ) and (C.8), we have the final result for the reduced matrix element. 



4"^^ = j drr^ {{00 I (g){f{T) \\)M 



(r)(||i(r))^|10)) 



2/^ 2mN 



K 



dr r 



(C.9) 
(C.IO) 



C.2 Including OPE Details 



C.2.1 Reducible OPE 



Taking just the q' terms at the OPE vertices. Fig. |5.3[ b) is given by 

9 A \ 
'2/J 

X ri,3 —p:, CTi • ( 2p , 



M'{p,k) 



-Ed - p^/m/vT 



q'2+/i(0)2 

m. 



(72 • q 



2m 



(C.ll) 
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where = — uP' . Adding to this expression emission from the right nucleon and 

approximating = rriT^mj^ as discussed at the end of Sec. 5.3, we find 



(00 I >i'(p,k) I 10) = • • q- ,., ^ V.. S • p, 



(00 I 7W'(r) I 10) 



8/^ 
87r/3 2mN 



Kof iSi2fiuj,r) + ■ a2giu,r)) 



1 



-Ed - 



S-V, 



(C.12) 

where the V acts on the initial np wave function, 5*12 = 3cri • fa2 ■ r — ai ■ a2 is the normal 



tensor operator, and 



f{u;,r) 



1 + 



+ 



(C.13) 



come from the Fourier transform [see Eq. (C.3)] of the pion propagator. Next, we evaluate 



u{r) 



where 



I +^(^1^1 I ) {S^2f{u:,r) + a^•a29{u:,r)) ^ ^^eS^ \ +^^^{^0, |, 



u{uj, r) 



w{uj, r) 



(C.14) 



u(uj,r) u(r) , , ^w(r) „, 

^ = ^ r) - 2/(a;, .)) + 2^2^/(0;, r). 
r r r 



Thus, 



(C.15) 



(/(r) II (5i2/(w,r) + ,Ti-<T2<7(^,r))S- V II i(r)) 

= 4vrzf^^/2f# + ^V^^ 
\ r \or r J r 

and we finally arrive at the full reduced matrix element, 

.3 _ ^(0)3 



d_ 
dr 



A, 



OPE,red,f 




2/3 2miv -Ed - m., 



-L/(0), 



UJ 



dr r 



uiuj^^ _^ 2\ w(aj,r ) f d 



dr r 



dr 



R,{r). 



(C.16) 



(C.17) 
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C.2.2 Irreducible OPE 



Finally, as described in Sec. 5.3, for the irreducible diagram we use [— 771^^/2 — -p^ /2m n) ^ w 
(— m7r)~^ for the intermediate nucleon propagator and take oj = mT^/2. 

,OP.«^_,|^M^,,,„,/„, ,,,,, 

C.2.3 Initial state OPE 

For OPE in the initial state, the isospin matrix element is (00 | ti^ti ■ T2 \ 10) = 1, and 
because the initial state consists of just one channel, '^Pi, 

{Si2f{u,r) + ai-a2g{oj,r)) f P,) = {2f{u,r) + g{u;,r)) f Pi). 

Evaluating the S • ^ reduced matrix elements, we find 



U{uj) = J drr^ (^V2 



C.3 Exact Wave Function Corrections Details 

Consider the nucleon propagator for reducible OPE in the initial state. Pulling out a —m^ 
and expanding this function in spherical coordinates, we have 



k2 - ^2 - ie 

e«C|r-r'| 

= -2imiv— (C.20) 

47r|r — r'l 

l,m 

where ^ = y/rr^/niN and r<(r>) is the lesser (greater) of |r|, |r'|. This spherical partial wave 
expansion was derived from the differential equation 
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where iGo = —irriNQ. First, one solves the homogeneous equation and requires both 
finiteness at the origin and outgoing wave behavior for large r. Thus, Q{r,r') = 
^ji(.^r<)/ip^(^r>). Next, the boundary condition at r = r' is obtained by integrating 
the differential equation across the boundary. In terms of g{r,r') = rr'Q{r,r'), 



d 



9>{r, r')\ 



d 



9<{r, r'% 



-1, 



(C.22) 



dr 

which yields A = i^. At this point in the diagram, the two-nucleon state is still ^Pi, so we 
preform one of the angular integrals and obtain 

iGo(r,r') ^ mNih{ir<.)h^^\ir-^). (C.23) 



Thus, 



OPE,red,i 




(m,/2) = - N 



2/3 2mN 3 
V{m^/2) = j drdr'r^r'^ (^^r 



d u{r) f d 3\ w{r) 
r \dr r I r 



(C.24) 



X ji{Cr^)h^^\^r>) {2fimj2,r') + g{mj2,r')) R,{r'). 

For the irreducible initial state OPE, the only difference is that a —2mM gets pulled out 
and the momentum becomes ^' = y^2m^j^rri/v, 



^OPE,irr,ij 



m„/2) 



-N- 



2/3 2mN 



-2imNi'L\m^/2)) 



(C.25) 



For the final state OPE, we can obtain the correct Green function from Eq. (C.23) by 
letting ^ — )• and using the correct / for the term under consideration. 

C.4 Cutoff Details 



In this section we display the exact expressions needed to implement the Gaussian cutoff of 



Sec 5.5 For the OPE diagrams, the integral of Eq. (5.9) is evaluated, 

g-/i(aj)r 



5A(a;,r) = ie'^HVA^ 



//(a;)r 



■erfc 



Ar fi{uj)\ e'^('^)'" fAr fi{uj) 



One also needs derivatives of Eq. (C.26), 



<Ti • Vcr2 • V5A(w,r) 



(5'i2/A(a;, r) + cTi • <T2^A(a;, r)) , 



(C.26) 



(C.27) 
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where, 

/A(w,r) 



A 



1 + 



+ 



+ H and A -A) 



i^lJi{uj)r + IH 



erfc 



^ /^(^) 
' 2 A 



'I 2 + A ; 



IJi(oj)r 



erfc 



' 2 A 



A 



A2 



+ {fi^ fi and A -A) . 



^(a;)r + 1 e V 2 + a ; 



g-^j(a;)r 



(C.28) 



(C.29) 
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Appendix D 
CH. H DETAILS 



D.l N ^ Nir from BxPT 



Recall that the LO NNtt interaction reads, 



£{0) c N^-A^^j^N, 



(D.l) 



where n_L^^ = Uf^ — v ■ uv^, = i{v)dfj,u — ud^u'^), and = e!"^°-'^a-l ^-^ , Expanding, we find 



i = 1,2,3 



(D.2) 



a o 



and thus the Feynman rule for an outgoing pion with momentum q and isospin a is 

= -i( :?^7i7N (^(7,:)r„. (D.3) 



2/. 



At threshold, the pion four-momentum is g = {m-n-, 0, 0, 0) making o'k^ = 0. This reflects 
the fact (well-known from current algebra) that threshold pion production proceeds via the 
off-diagonal, and therefore l/m^ suppressed, interaction g-w^'^^f^q^r. In the effective theory, 
this recoil correction shows up in the NLO Lagrangian 



where the spin vector is = — 57^(7'^ — v^). Thus the Feynman rule is 



(D.4) 



-I 

2mAr / 



la,. X 
- — [-imn) 

J TV 



-I- :^7"7i7'^(^ - ^)i-ra, 



"ivr OA 5 0/ 



(D.5) 



2mN 2/^ 

where the derivatives act on the nucleon wave functions. 
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D.2 One Pion Exchange Deuteron 

In this Appendix, we present the method by which the deuteron wave function is calculated 



for use in Sec. 6.5 This method is taken directly from the work of Friar, Gibson, and Payne 



|105j . The OPE potential is defined to have central (Y) and tensor (T) parts, 

K(r) = /'m^^^ [a, • a^Yir) + SuTir)] , (D.6) 

where = 0.079 (to be distinguished from /tt) measures the strength of the pion-nucleon 
coupling and 5i2 is the standard tensor operator. The deuteron has isospin zero and spin 
one, so we have 

K(r) = -fm^ [Yir) + Si2T{r)] . (D.7) 

The Y and T functions are expressed as derivatives of the Fourier transform of the pion 
propagator, 

Y{r) = /io(x) - hQ{x)/x 

T{r) = h'^{x) + 2h'o{x)/x (D.8) 
where x = mT^r and -FjrAfAf is the form factor for which we use, 

(A 2 2 \ 



In Ref. [105j . it is shown that 

_™ 2n— 1 



Y{r) = — - /S^-f'^ V § {5i{Px) - 2i6i.i{(3x)) 

i=0 



(D.IO) 



i=0 

32 t\/r,Q2 



where f3 = A/m-,^ and ^ = (/3 — l)/2/3 and the 6i are defined by 

= {2i - l)6i{(3x) + (/3x)25i_i(/3x) (D.ll) 

along with 6o = l/(3x and 5i = 1. One of the results of Ref. |105| is that larger values of n 
lead to better fits to experimental data. We use n = 5 and /3/\/To = 5.687805 in order to 
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Uiir) (fm"""^) 




Figure D.l: Deuteron s- and d-state wave functions (the s-state is larger). The potentials 
used to calculate the wave functions are Avl8 (black) and the cutoff OPE described in this 
section (red). 



Table D.l: Deuteron properties. 

Potential Q (fm^) (r2)V2 (fj^) 

Avl8 0.270 1.968 

OPE (n = 5) 0.282 1.939 

Experiment 0.2859(3) 1.955(5) 



precisely reproduce the binding energy Eb = 2.2246 MeV. The wave functions are calculated 
by integrating in from r^ax = 100 fm and adding together two linearly independent solutions 
such that the sum vanishes at rmin = 0.01 fm. As shown in Fig. D.l, the results are close 
to the "correct" Avl8 deuteron. 



In Table D.l we display the quadrupole moment and mean square charge radius of Avl8, 
this OPE potential, and experiment (as quoted in |105) ). It is clear that the form factors 
in the OPE potential make it difficult to distinguish this construction as less accurate than 



Avl8. Finally, in Table D.2 , we display the reduced matrix elements for the rescattering pion 
production diagram evaluated with both the phenomenological potentials and the deuteron 
of this section. Since this diagram makes the largest contribution to the cross section we 
need to verify that neglecting non-OPE parts of the potential does not dramatically change 
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Table D.2: Effect of using OPE deuteron on rescattering diagram. 



Deuteron 


Avl8 


Reid '93 


Nijm II 


Phenomenological 


67.8 


69.7 


71.1 


OPE (n = 5) 


69.8 


72.1 


74.0 



this amplitude. Indeed, we observe what should be expected: since the rescattering diagram 
is not as sensitive to the core of the deuteron, using the OPE wave function in place of the 
standard one has only a small effect. 

D.3 Effect of the 6 Terms: Diagram 

In this section we calculate the correction terms to the first two-body DW amplitude [Eq. 
(6.40)] which is shown in Fig. 6.9(a)[ Assuming that the J's truly are small compared to Eq. 



(6.42 ), we will only worry about calculating them one at a time, numbering the contribution 
of the (5's from right to left as 1, 2, and 3. Note that we will display the results as calculated 
using the OPE deuteron and the Avl8 initial state. 

D.3.1 First Og^ correction term: AOi 



Consider the rightmost 5 in Eq. (6.40) 



\n CTi • (p - k)cr2 • (k - p) 2 I, 



2miv 2 2 Q 

2 ~ ■ P' 



(D.12) 



where /i^ = 3m^/4 and Fj^^N is the form factor described in Appendix D.2 The easiest 
way to evaluate the matrix element of this operator is to let the OPE act to the left on 
the deuteron in position space. The resulting expression is then transformed to momentum 



space. We can expand the fraction in the integrand of Eq. (D.12) into spherical harmonics 
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(taking pi = z) 



2mN 2^2 / 

and note that only the / = 0, 2 terms will contribute to s-wave production. The expansion 



coefficients are shown in Fig. D.2 Clearly the / = 2 term is small and we neglect it here to 



avoid the extra algebra involved with a J = 2 operator (resulting in the A2 reduced matrix 



elements in the notation of Appendix B.l). We find 



^•^1 = -34%. (D.14) 



M 

D.3.2 Second correction term: AO2 

Next consider the term, 

An - '^i • (P - 1^)'^2 • (k - p) 2 ,,x2x 1 1 

V2mAr 2 2j ^ 
This term has a modified OPE, 

d^q jq.^ cri -qcra- (-q) ^ . 2x 1 ^C(?') 



(D.15) 



(27r)3^ q2+^2 -^^^^^"^7^2^- ^ ^ ^ ^ 4vr 



^[5i2rc(r) + ^i.(T2ycW]- 



(D.16) 
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In Fig. D.3 we compare the functions T^{r) and Y(-{r) to traditional OPE which has /x in 



place of the square root in Eq. (D.16). We use the Schrodinger equation to replace the 
p'^/2mN with V{r) and evaluate the matrix element in position space to find 



M 



+50%. 



(D.17) 



D.3. 3 Third O^^ correction term: AO3 

Calculating the effects of the 6 in the denominator of the OPE is difficult to do exactly due 
to the combination of momenta that appear, 



AO. 



<Ti • (p - k)o-2 • (k - p) 



(p - k)2 + m2 + 5(p) +m2/4 



-F^nnUp-^) 



-Sp, 



(D.18) 



(recall that 6{p) = /2mN — 771,7^/2 + 7/2). Instead we will evaluate it for fixed values of 5 
which represent the deviation of p away from pj, 



d+ = d{pi + m,r) = = 0.45 

X xf \ -2pim„ + ml 

0- = oipi — mn) = = — O.ozmTr. 



(D.19) 



The modified tensor and central functions and 1^ are shown in Fig D.4 We define the 
correction as 



AM3i,6) = M{6) - M{0). 



(D.20) 
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Figure D.4: Effect of the 6 on the OPE (a) tensor and (b) central radial functions. 



We find, 



+16% 



M 

D.3.4 Summary of corrections 

For the purposes of estimating the net result we take the average of the estimates in Sec. 
EM] and find 

« -34% + 50% - 8% = +8%. (D.22) 

M 

We have successfully shown that the corrections to the first two-body DW amplitude are 
small, and actually increase the amplitude which is already twice as large as the traditional 
impulse approximation. 

D.4 Effect of the 6 Terms: Diagram 

The second two-body DW amplitude's corrections are evaluated exactly as in the previous 
sub-sections and we just display the results here, 

M 

= -39% (D.23) 

M 

M = M = 
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\ 

\ 

\ 

\ 

\ 




Figure D.5: Impulse approximation with distorted waves: initial state heavy meson ex- 
change. Solid lines represent nucleons, dashed lines pions, and the double solid line a a 
meson. Crosses represent propagators which are absent due to the partial fractions decom- 
position. 



with the net result 



AM 



tot 



M 



-35% - 39% + 3% 



-71%. 



(D.24) 



We see that the corrections to the approximation in Eq. (6.43) are fairly large, but this has 



a negligible effect because the amplitude is already small compared to the first two-body 
DW amplitude. 

D.5 Heavy Meson Exchange 



Consider the loop on the left-hand side of Fig. D.5 which is obtained by using OPE for 



the left K in Eq. (6.31) and a exchange (the dominant intermediate-range mechanism) for 



the right K. Note that this loop only differs from Fig. 6.6 in two ways: the meson- nucleon 
vertex (here we consider only scalar-isoscalar) and the meson mass. We use a typical set of 
parameters [29], g'^/^ir = 7.1 and mo- = 550 MeV. 

The result of integrating over energy will proceed exactly as it did with the pion resulting 

To interpret the IgB^ term (again. 



D.5 



in the two diagrams shown on the RHS of Fig. 
neglecting stretched box diagrams), we absorb the sigma exchange into the initial state and 
no new term is added. However, in the term, after absorbing the pion exchange into 
the final state, we are left with a new operator. The amplitude for this operator can be 
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obtained from that of Fig. 6.9(b) with the following change, 



where iJ? = 3m^/4 and f^l = ml - (m^/2)^. We find, 



>gl^i^ , (D.25) 



= -7.24, (D.26) 

which is larger in magnitude than the pionic ^2B^'* (with the same sign) but smaller than 
^2B^'" (with the opposite sign). Since nia is relatively large, we can safely ignore the two 6 
corrections that are competing with uji and only need to evaluate 

= -27%. (D.27) 

M 

One natural question is whether the static a exchange already present in the initial-state 
wave function is a sufficient approximation for the contribution considered in this section. 
To answer this question, we can evaluate the traditional impulse approximation with 

= \pi,P2) + GV„\^)u (D.28) 

where here we employ a static a exchange that is present (at least effectively) in the wave 
function. Using this initial-state wave function, we calculate 

Ml^ = f{<P\0^\^)1, (D.29) 

and find the reduced matrix element, 

Al^ = -3.3. (D.30) 

Thus we see that the a exchange in the traditional impulse approximation is an underesti- 
mate (in magnitude) of the true non-static exchange dictated by the loop integral. 

Of course there is no a in traditional BxPT, so this section is simply telling us that to 
achieve high accuracy it is indeed important to use more than just simple pion exchange 
when forming the original box diagram. Such a calculation is beyond the scope of this work. 
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